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Abstract 


The  purpose  of  this  thesis  was  to  model  a  computerized  tomography  device  (CT) 
using  Monte  Carlo  methods  to  determine  the  scattered  radation  spectrum  inside  and  die 
dose  outside  the  suite  at  Wii^t-Patterson  Medical  Center  and  at  a  generic  suite  to  non- 
occiqiational  personnel.  This  was  driven  by  the  recent  inclusion  of  die  most  recent  NCRP 
recommended  dose  limits  into  10CFR20  of  50  mSv,  occiqiational;  1  mSv,  non- 
occupational  continuous  ejqiosure;  and  5  mSv,  infrequent  exposure.  The  rotating  fan 
beam  of  the  CT  was  modeled  for  MCNP,  a  general-purpose  Monte  Carlo  n-particle 
transport  model.  The  CT  target  was  a  standard  human  phantom  defined  in  MIRD 
Pamphlet  5.  The  source  and  phantom  were  rotated  27**  widi  respect  to  the  walls.  A  head 
scan  and  an  abdomen  scan  were  simulated.  CT  usage  factors  were  set  equal  to  average 
WPMC  values.  The  suite  walls  were  modeled  alternately  with  1/16”  and  1/8”  lead 
shielding  between  gypsum  drywall  for  the  occupational  and  non-occiqiational  dose, 
respectively.  With  CT  isocenter  as  problem  origin,  test  detectors  were  placed  30  cm 
outside  the  walls  on  two  sets  of  axes  defined  by  (a)  the  walls  and  (b)  the  source-body  axes. 
Film  badges  placed  in  die  CT  suite  at  WPMC  agreed  with  MCNP  modelir^  results  by  a 
factor  of  two,  validating  the  method.  The  analysis  showed  that  outside  bodi  the  WPMC 
CT  suite  and  the  generic  room  the  continuous-exposure  non-occupational  dose  limit  was 
exceeded  below  die  floor  and  above  die  ceiling;  the  infrequent  exposure  non-occupadonal 
dose  limit  was  exceeded  below  die  floor.  Modeled  results  exceeded  the  occiqiational  and 
non-occupational  dose  limits  outside  two  walk  of  die  generic  suite.  Modeled  dose  was 
below  the  occupational  limit  at  all  test  locations.  The  scatt^d  radiation  spectrum  is  a 
softened  source  spectrum  minus  die  60  and  68  keV  X-ray  tube  bremsstrahhing  peaks. 
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A  MONTE  CARLO  ANALYSIS  OF 


COMPUTERIZED  TOMOGRAPHY 


1.  Introduction 

Tomography,  or  sectional  radiography,  is  defined  as  “making  radiographs  of  plane 
sections  of  a  body  or  an  object,  its  purpose  is  to  show  detail  in  a  predetermined  plane  of 
the  bo^.”  (Parker,  1989)  Computerized  tomogr^^hy  (CT)  uses  computers  to  process  the 
many  image  projections.  One  means  of  obtaining  a  CT  image  is  by  passing  X-rays  tinough 
the  object  (usually  a  human  body)  to  be  imaged.  The  by-product,  however,  of  this  type  of 
imagery  is  scattered  radiation.  Because  of  the  potentially  harmful  ejSects  of  radiation  on 
the  human  body,  it  is  absolute^  necessary  to  know  the  dose  received  by  hospital  occupants 
fi-om  scattered  radiation  produced  by  X-ray  imagery.  Recent  changes  to  U.S.  Federal 
Law,  Title  10  Code  of  Federal  Regulations  Part  20  or  10CFR20  (NRC,  1995),  have 
required  additional  efiTorts  to  characterize  radiation-producing  systems.  Direct 
measurement  of  the  scattered  radiation  can  be  difficult  and  only  tests  the  radiation  field  at  a 
few,  selected  points.  Computer  modeling,  however,  is  a  non-destructive,  safe  metirod  of 
interpolating/extrapolating  the  dose  from  scattered  radiation  fi'om  a  variety  of  radiation 
producing  systems  to  other  locations. 
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1.1  Background 


The  International  Commission  on  Radiation  Units  and  Measurements  (ICRU), 
since  its  inception  in  1925,  and  the  National  Council  on  Radiation  Protection  and 
Measurements  (NCRP),  chartered  in  1964  by  Congress,  are  responsible  for  promulgating 
information  and  recommendations  based  on  leading  scientific  judgment  in  matters  of 
radiation  protection  and  measurement.  NCRP  Report  49  contains  tiie  most  recent 
comprehensive  shielding  guidance  for  diagnostic  X-ray  facilities  (NCRP,  1976).  The 
shielding  guidance  was  based  on  NCRP’s  then-recommended  annual  dose  Umits  of  SO 
mSv^  and  5  mSv  for  occi^tional  and  non-occupational  personnel,  respectively^.  As  earl^ 
as  1984,  NCRP  reduced  the  non-occiq)ational  dose  limit  recommendation  to  1  mSv  for 
continuous  e^)osure  (retaining  the  non-occupational  dose  limit  of  5  mSv  for  infirequent 
exposure).  In  Januaiy  1994  the  most  recent  NCRP  dose  limits  were  included  in  Titie  10 
Code  of  Federal  R^;ulations  Part  20.  With  the  recommendations  now  law,  there  is 
obvious  concern  whether  facilities  built  according  to  NCRP  Report  49  guidelines  and  the 
older,  greater  dose  limits  will  meet  the  new  limits. 

Robert  Met2ger,  et.  aL  used  die  Monte  Carlo:  N-Particle  transport  program 


^  Gray  (Gy)  and  sievert  (Sv)  are  the  special  units  of  dose  and  dose  equivalent,  respectively, 
adopted  by  the  ICRU  in  1 979  to  conform  to  the  SI  unit  standard.  They  replace  rad  and  rem:  1 
Gy  =  1  J/kg  =  100  rad  and  1  Sv  =  100  J/kg  =  100  rem  thus  50  mSv  =  5  rem,  etc. 

^  Occupational  personnel  are  assumed  to  spend  their  entire  work  period  in  controlled  areas 
defmed  as  areas  in  which  the  e^qjosure  of  persons  to  radiation  is  under  the  supervision  of  a 
Radiation  Protection  Supervisor  (implying  control  of  access,  occupancy  and  working  conditions 
for  radiation  protection  purposes).  Non-occipational  personal  are  assumed  to  be  restricted  to 
uncontrolled  areas  defmed  as  any  area  not  meeting  the  definition  of  a  controlled  area  (hallways, 
offices,  waiting  rooms,  etc.).  Patients  are  not  included  in  either  of  these  definitions  with  the 
assumption  that  the  dose  received  is  acceptable  given  the  benefits  of  the  treatment  or  imagery. 
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(MCNP)  to  analyze  the  effectiveness  of  NCRP  Report  49  shielding  recommendations  in 
standard  diagnostic  X-ray  rooms  (Metzger,  1993).  They  note  that  typical  shielding 
calculations  ignore  absorption  by  die  patient,  film  cassettes,  table,  and  fi-equentty,  the 
gypsum  wallboard  because  of  the  overall  complex  geometiy  involved  and  lack  of  material 
attenuation  data  in  NCRP  Report  49.  Using  MCNP  for  shielding  calculations  eliminated 
these  difficulties. 

MCNP  (1993)  is  a  general-purpose  Monte  Carlo  N-particle  code  (N:  many)  for 
calculating  the  time-dependent  energy  transport  of  neutrons,  photons,  and/or  electrons  in 
user-defined  three-dimensional  geometry.  The  code  uses  die  Monte  Carlo  methods 
derived  firom  research  done  during  die  Manhattan  Project  in  the  earfy  1940’s.  Its  strength, 
and  thus,  suitability  for  shielding  anatysis,  comes  fi'om  four  key  features:  (1)  the  ability  to 
model  complex  geometry,  (2)  the  wide  selection  of  variance-reduction  techniques,  (3)  the 
inclusion  of  extensive  cross-section  libraries  (allowing  the  use  of  materials  of  any 
composition),  and  (4)  the  abOity  to  model  sources  of  complicated  shape  or  spectra. 
(Briesmeister,  1993). 

Metzger’s  analysis  revealed  that,  in  general,  a  standard  diagnostic  X-ray  suite 
configuration^  shielded  according  to  NCRP  Report  49  guidelines  would  meet  the  new 
standards  assuming  modem  X-ray  equipment,  usage  and  films.  They  clearfy,  however, 
stated  that  die  final  result  could  not  be  extrapolated  to  CT  units  (Met:%er,  1993).  Modem 
CT  units  use  a  fan  shaped  beam  diat  rotates  around  die  patient  with  the  imaging  material 


^  A  standard  diagnostic  X-ray  suite  configuration  would  consist  of  a  sipine  patient  on  a  table  or 
an  upright  patient  with  a  wall  mounted  film  cassette  both  with  simple  conic  or  pyramid  shaped 
beams. 
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(detector)  rotating  opposite  the  source.  Given  this  radical^  different  source,  patient  and 
detector  configuration,  Metzger’s  assertion  is  understandable.  Their  method  of  anafysis, 
however,  could  be  applied  to  any  type  of  diagnostic  or  therapeutic  radiation  units. 

1.2  Problem  Statement 

Dose  outside  the  CT  suite  at  Wright-Patterscm  Medical  Center  (WPMC)  needs  to 
be  estimated  and  compared  to  dose  limit  requirements  in  10CFR20.  Captain  William 
Ruck,  WPMC’s  radiation  protection  ofilcer  (Wright-Patterson  Air  Force  Base,  OH) 
suggested  this  anafysis  for  the  General  Electric  Ifi-Speed  Advantage  Computerized 
Tomography  unit  installed  in  the  medical  center’s  CT  suite.  He  recommended,  however, 
that  the  anafysis  be  generalized  to  a  standard  suite  configuration.  This  could  be 
accomplished  using  the  manufacturer’s  recommended  room  dimensions  and  performing 
the  analysis  for  a  variety  of  shielding  conditions,  i.e.,  varying  thickness'  of  lead  shielding  in 
the  walls  or  concrete  in  the  floors  and  ceiling.  Additionally,  research  indicated  tiiat 
information  about  the  scattered  spectrum  within  a  CT  suite  was  not  available  (Simpkin, 
1995)  which  limits  analysis  of  CT  dose  by  methods  other  than  those  used  here.  MCNP 
was  able  to  provide  an  estimate  of  the  spectrum. 

1.3  Approach 

First,  the  results  of  Metzger,  et  al.  were  dtqtlicated  to  ensure  their  work  and  method 
were  correctfy  understood  and  to  confirm  their  calculations.  Next,  a  CT  unit  (source, 
detectors,  and  shielding)  was  placed  in  Metier’s  X-r:^  room  with  its  primary  axes  parallel 
to  the  walls.  This  configuration  was  used  for  test  runs.  The  analysis  runs  included 
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modeling  Wrigjit-Patterson’s  CT  suite,  containing  the  CT  unit  rotated  with  respect  to  the 
walls,  with  various  degrees  of  shielding.  In  addition,  the  manufacturer’s  recommended 
room  size  was  modeled  with  die  rotated  CT  unit.  As  a  test  of  modeling  accuracy,  radiation 
film  badges  were  placed  in  Wii^t-Patterson’s  CT  suite  and  their  reported  eTqiosure 
conqiared  to  MCNP’s  predictions. 

Since  no  modeled  scenario  can  ever  perfect^  mimic  reality  without  unmanageable 
ovefhead  (input  detail  and/or  computer  resources),  tqiproximations  are  necessaiy.  Two 
factors  affect  the  validity  of  the  approximations:  how  well  the  model  duplicates  actual 
events  and  the  level  of  significance  a  particular  input  (and  ite  operating  range)  has  on  the 
desired  output.  These  factors  were  weig)ied  against  available  computer  resources 
(memoiy,  computation  time,  etc.).  The  ^proximations  and  limitations  and  their  affect  on 
the  results  were  thorough^  investigated. 

1.4  Contents 

Chester  2  discusses  X-ray  interactions.  Cluqiter  3  presents  the  Monte  Caiio 
mediod  and  its  application  to  particle  transport.  Ch<q)ter  4  contains  a  description  of  MCNP 
(the  code  used  for  die  analysis).  Chapter  5  describes  a  standard  CT  unit,  how  it  was 
modeled  for  MCNP  and  the  assumptions/limitations  used  in  die  modeling  process. 

Chapter  6  contains  the  calculations  and  methods  used  and  analysis  of  preliminary  test  runs. 
Chapter  7  presents  die  results  and  analysis  of  the  MCNP  runs  for  inside  and  outside  the 
CT  suite  at  WPMC  and  outside  the  G£  recommended  room.  Last,  Chapter  8  contains 
conclusions,  recommendations  and  possible  follow-on  anatysis. 
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2.  X-Ray  Interactions 


2.1  Introduction 

Diagnostic  X-ray  and  computerized  tomography  produce  dieir  images  by  using 
interactions  of  X-ray  photons  with  matter.  In  both  imaging  techniques,  the  target  to  be 
imaged  is  positioned  between  a  beam  (or  source)  of  X-rays  and  a  recording  material,  film 
or  detector.  The  image  is  produced  fi*om  die  measured  or  recorded  X-ray  transmissions. 
X-ray  energy  used  for  diagnostic  purposes  is  confined  to  25  to  1000  keV  to  maintain 
reasonable  resolution  widi  acceptable  attenuation  levels  (Macovski,  1983). 

There  are  three  primary  types  of  X-ray  interactions  which  must  be  considered  in 
radiological  physics:  Compton  effect,  photoelectric  effect  and  Rayleigh  (coherent) 
scattering.  The  first  two  are  important  since  they  result  in  die  transfer  of  energy  to 
electrons  which,  in  turn,  impart  that  energy  to  matter  (targets).  Ra}dd^  scattering, 
although  prevalent  in  the  diagnostic  range,  only  removes  the  X-rays  fi-om  die  beam  with  no 
odier  implied  interactions.  In  addition,  photonuclear  interactions  and  pair  production  are 

possible  X-ray  interactions  but  they  are  not  significant  in  the  diagnostic  energy  range.  The 

consequences  of  these  interactions  depend  on  both  the  photon  quantum  energy  (£  =  hv) 
and  the  atomic  number,  Z,  of  the  absorbing  media.  (Attix,  1986;  Chilton,  1984;  and 
Profio,  1979). 
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2.2  Compton  Effect 


The  Compton  effect  can  be  split  into  two  aspects:  kinematics,  relating  the  energies 
and  angles  of  the  particles,  and  cross-sections,  which  predict  the  probability  of  an 
interaction  occimii^  An  incoming  photon  collides  with  an  electron  that  is  assumed  to  be 
initially  unbound  and  stationary.  This  is  an  accurate  approximation  for  two  reasons:  the 
energy  the  electron  receives  is  much  greater  than  its  binding  or  kinetic  energy  and,  in  the 
diagnostic  range,  any  errors  tiiat  result  from  tire  approximation  win  be  much  less  than  the 
competing  photoelectric  effect.  (Attix,  1986;  Chilton,  1984;  and  Profio,  1979). 

2.2.1  Kinematics. 

Figure  2.1  shows  a  Qrpical  photon-electron  coUisiott  The  photon  has  an  initial 


e* 


Figure  2.  1  Kinematics  of  the  Compton  effect 
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energy,  E  =  hv  and  momentum,  P  =  hv/c.  The  electron  is  stationaiy.  After  the  collision, 
the  electron  departs  at  an  angle  6,  with  kinetic  energy  T  and  momentum  p.  The  photon 
scatters  at  an  an^e  (}>  with  a  new,  lower  energy,  hv’,  and  momentum,  hv’/c.  The  collision 
kinetics  solution  is  based  upon  conservation  of  both  energy  and  momentunt  Energy 
conservation  requires 

r  =  hv-hv’.  (2.1) 

Conservation  of  momentum  in  the  original  photon  direction,  0°,  requires 


hv  hv’  ' 

- = - COS(|)+  /7C0S6 

c  c 

and  perpendicular  to  the  direction  of  incidence  requires 


(2.2) 


hv'sin<()  =  pccosG.  (2.3) 

Using  relativistic  relationships,  pc  can  be  e^qnressed  in  terms  of  T  and  rtio  (electron  rest 

mass) 


pc  =  [T(T+2in„c’)]  . 

Eqs.  (2.1),  (2.2),  and  (2.3)  (with  the  (2.4)  substitution)  are  a  set  of  three  simultaneous 
equations  in  hv,  hv’,  T,  6  and  (|).  A  complete  solution  to  die  kinematics  of  Compton 
interactions  can  thus  be  derived  given  any  two  of  the  five  parameters 

hv’= _ 5^^ _ ; 

l+(hv/m^c^Xl-w>8(t))  ’ 


T=  hv-hv'; 


cotG  = 


(2.4) 


(2.5) 

(2.6) 

(2.7) 
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where  the  rest  mass  energy  of  the  electron,  moc^,  is  0.511  MeV  and  hv,  hv’  and  T  are  also 
expressed  in  MeV. 

£q.  (2.5)  s^ows  ^t  hv,  hv’  and  an^e  ()>  of  the  scattered  photon  are  unique^ 
correlated.  Eq.  (2.6)  then  provides  &e  electron’s  kinetic  energy  while  Eq.  (2.7)  gives  its 
scattering  angje,  6.  Figure  2.2  shows  the  relationship  between  hv,  hv’  and  T  for  a 
selection  of  photon  scattering  angles,  (|),  at  energies  in  the  diagnostic  range.  Figure  2.3 
shows  the  relationship  between  the  scattering  angles,  ({>  and  6,  for  a  selection  of  photon 
energies,  hv.  (Attix,  1986;  Chilton,  1984;  and  Profio,  1979). 


0.01  0.1  1 

hv,MeV 

Figure  2.  2  Kinematic  relationship  of  hv,  hv*  and  T.  T  can  be  interpreted  as  the 
vertical  separation  of  the  (|>  curve  from  the  (|)  =  0  diagonal. 

2.2.2  Cross  Section 

J.  J.  Thompson  derived  the  earliest  theoretical  description  of  photon-electron 
scattering.  He  assumed  tire  electron  to  be  freely  oscillating  under  the  influence  of  an 
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Figure  2. 3  Scattering  angle  relationships  in  Compton  scattering. 


electric  field  and  re-emitting  photons  of  die  same  energy  as  the  incident  photon.  The 
electrcm  will  thus  retain  no  kinetic  energy  firom  the  photon.  By  1928  it  was  known  that 
Thompson’s  value  was  too  large  for  hv  >  10  keV,  thus,  Klein  and  Nishina  q>plied  Dirac’s 
relativistic  theory  of  the  electron  to  the  Compton  effect  They  also  assumed  unbound 
electrons  initially  at  rest  £q.  (2.8)  could  tiren  be  re-written  as 

2  ^hv/  vhv'  hv 


clf.o 


2  Vhv/  vhv'  hv  / 


(2.8) 


where  hv’  is  given  by  Eq.  (2.5).  Figure  2.4  shows  this  relationship  for  a  selection  of 
incident  photon  energies,  hv.  The  total  Klein-Nishina  cross  section  per  electron  (eC, 
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typically  called  the  Compton  scattering  or  incoherent  scattering  cross  section)  can  be 
obtained  by  integrating  Eq.  (2.10)  over  all  photon  scattering  angles 


<* 

1+a 

’2  (1+a) 

ln(l+2a) 

ln(l+2a) 

-v 

l+3a 

- 

l+2a 

a 

^  2a 

(l+2a)^ 

(2.9) 


where  a  =  hv  /  mo  c^.  The  Klein->nshina  cross  section  per  electron,  e  a ,  in  Eq.  (2. 1 1)  is 
independent  of  Z.  The  Klein-Nishina  cross  section  per  atom,  a<7,  is 

y^o=ZxEO  (cmVatom)  (2.10) 


Equivalent!^,  the  Klem-Nishina  cross  section  per  mass  is  obtained  from 

a  N.Z 

(cmVg)  (2.11) 

P  A 

where  Na  is  Avogadro’s  number,  A  is  atomic  weight  and  p  is  the  densi^  in  g/cm^ 

(Na  Z  /  a  is  the  number  of  electrons  per  gram  of  material),  o/p  is  called  the  partial  mass 
attenuation  coefficient,  or  less  frequently,  the  Conq)ton  mass  attenuation  coefficient 
The  cross  section  (which  can  be  thou^t  of  as  an  effective  target  area  for  the 
photon)  is  numerically  equal  to  die  probability  of  a  Thompson  scattering  event  occurring 
when  a  single  photon  passes  through  a  layer  containing  one  electron  per  cm^.  It  is  also  the 
fraction  of  photons  that  will  scatter  in  passing  througli  that  same  hQrer,  e.g.,  approximately 
665  events  for  10*’^  photons.  (Attix,  1986;  Chilton,  1984;  and  Profio,  1979). 
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Figure  2.  4  DifTerential  Klein-Nishina  cross  section  vs.  scattered  photon  angle,  <t>,  for 
hv  =  0.01, 0.05,  0.1, 0.5,  and  1.0  MeV. 

2.3  Photoelectric  Effect 

The  photoelectric  effect  is  die  most  important  interaction  of  low-eneigy  photons 
widi  matter  because  the  photoelectric  cross  sections  increase  strongj^  as  photon  energy 
decreases,  particulaily  for  high  Z  media.  Again,  there  are  kinematic  considerations  relating 
energies  and  angles  of  the  participating  particles  to  die  cross  section  predicting  the 
interaction  probability.  (Attix,  1986;  Chilton,  1984;  Profio,  1979). 

2.3.1  Kinematics. 

The  photoelectric  effect  (for  X-ray  interactions)  differs  from  die  Compton  effect  in 
that  the  photon  always  gives  up  all  of  its  energy  to  the  electrcm  it  collides  with  because, 
unlike  the  Compton  effect,  the  encounter  is  with  a  ti^dy  bound  electron,  such  as  diose  in 


2-7 


the  atom’s  inner  shell,  especially  of  hi^  atomic  number.  Figure  2.5  shows  a  typical 
photoelectric  interaction.  An  incident  photon  of  quantum  eneigy  hv  is  shown  interacting 
with  an  atomic-shell  electron  widi  a  binding  eneigy,  Et.  The  photon  will  be  totally 
absorbed  by  die  electron,  giving  it  enough  eneigy  to  overcome  its  binding  energy.  The 
electron  will  depart  at  angle  6  with  kinetic  eneigy,  T ,  and  momentum,  p.  The  atom  will 
depart  at  an^e  (|)  with  its  kinetic  eneigy,  T„  veiy  nearfy  0.  The  photon  eneigy,  hv,  must  be 
greater  than  the  electron's  binding,  Eb,  but,  the  closer  the  two  eneigies  are,  the  greater  the 
probability  of  a  photoelectric  interaction.  The  kinetic  eneigy,  independent  of  its  scattering 
angle  6,  is 

T=hv-Eb-TaShv-Eb.  (2.12) 


Figure  2.  5  Kinematics  of  the  photoelectric  effect 
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2.3.2  Cross  Section. 


Derivation  of  the  interaction  cross  section  for  the  photoelectric  effect  is  difficult 
because  of  the  complications  arising  from  the  binding  of  tiie  electron  (Attix,  1986).  There 
is  no  equation  analogous  to  tiie  Klein-Nishina  formula.  Published  tables  are  based  on 
experimental  results  supplemented  by  interpolations  for  other  energies  and  absorbing 
media. 


Photoelectrons  are  ejected  predominate^^  sideways  for  low  photon  energies.  With 
increasing  energy  the  distribution  is  pushed  more  toward  smaller,  but  still  non-zero,  angles. 
A  convenient  representation  uses  the  “bipartisan”  angje,  i.e.,  the  median  ejection  angle. 
Integrated  over  all  an^es  of  emission,  tire  cross  section  can  be  written 


hvsE^ 

otherwise 


(cm*/atom). 


(2.13) 


where  k  is  a  constant,  n  is  =  4  at  hv  =  100  keV,  gradually  rising  to  about  4.6  at  3000  keV, 
and  m  =  3  below  hv  =  100  keV,  gradually  decreasing  to  about  1  at  5000  keV.  Using  the 
same  conversions  as  in  Eqs.  (2.12)  and  (2.13)  for  hv  <  100  keV,  the  photoelectric  cross 
section  per  unit  mass  is 

(cm^/g)  (2.14) 

p  Vhv> 

This  is  called  tire  photoelectric  mass  attenuation  coefBcient. 

When  an  electron  is  removed  from  an  inner  atomic  shell,  the  resulting  vacancy  is 
pronptly  filled  by  another  electron  from  a  less  tightty  bound  shell.  This  may  be 
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accompanied  by  the  emission  of  a  fluorescence  X-ray  of  quantum  energy  equal  to  flie 
difference  between  the  two  energy  levels.  Binding  energy,  £b,  not  complete^  removed  by 
a  fluorescence  X-ray  may  be  given  to  one  or  more  ejected  electrons.  This  alternative  excess 
binding  energy  disposal  mechanism  is  called  the  Auger  effect  (Attix,  1986).  An  atom  may 
emit  a  number  of  Auger  electrons  more  or  less  simultaneously  in  a  kind  of  chain  reaction, 
resulting  eventually  in  more  than  one  vacancies  in  the  outermost  shell. 

2.4  Rayleigh  (Coherent  Scattering 

Rayleigh  scattering  is  an  elastic  collision  between  a  photon  and  particle.  The 
photon  loses  essentially  none  of  its  energy,  it  is,  however,  redirected  through  some  small 
angle.  Rayleigh  scattering  removes  photons  from  a  collimated  beam,  but  does  not 
contribute  significant  to  total  dose  since  no  energy  is  given  to  charged  particles,  and  no 
ionization  or  excitation  is  produced.  The  scattering  angle  dependents  on  both  Z  and 
photon  energy  giving  a  Rajdei^  cross  section  per  atom  of 

72 

a®R* - 2  (cmVatom),  (2.15) 

(hv)^ 

or,  in  mass  units 

—  (cm'/g)  (2.16) 

P  (hv)^ 

Ra)dei^  scattering  is  nearly  insignificant  vrdien  conqrared  to  Are  Compton  or 
photoelectric  effects  for  high  Z,  high  energy  interactions.  It  increases,  however,  for  low  Z, 
low  energy  events  such  as  diagnostic  X  rays  in  tissue,  composed  primarily  of  water,  H2O. 
(Attix,  1986;  Chilton,  1984;  Profio,  1979). 
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2.5  Other  Interactions 


There  are  two  other  photonic  interactions:  photonuclear  and  pair  producticm.  In  a 
photonuclear  interaction  an  energetic  photon  (exceeding  a  few  MeV)  enters  and  excites  a 
nucleus.  Required  photon  eneigies  are  well  above  the  diagnostic  range.  Pair  production  is 
an  absorption  process  in  which  a  photon  disappears  and  gives  rise  to  an  electron  and  a 
positron.  A  minimum  photon  energy  of  2moC^  =  1.022  MeV  is  obviously  required,  again, 
out  of  the  diagnostic  range.  (Attix,  1986;  Chilton,  1984;  Profio,  1979). 


2.6  Attenuation  Coefficients 


All  of  tiiese  effects  are  present  in  different  magnitudes  depending  on  X  ray  energy 
and  target  material.  Mass  attenuation  coefficients  can  be  combined  (neglecting 
photonuclear  and  pair  production)  to  obtain  a  total  mass  attenuation  coeffici^t 

-=-+-+^  (crnVg),  (2.17) 

P  P  P  P 


where  x/p  is  the  contribution  of  the  photoelectric  effect,  o/p  is  the  Compton  effect  and 
Or/p  is  for  Rayleigh  scattering. 

Since  all  target  media  are  not  simple  elements,  the  Bragg  rule  can  be  applied  to 
mass  attenuation  for  compounds  or  intimate  mixtures  of  eluents 


\p) 


P  f 

—  Ia+  — 


mix 


PV 


(2.18) 


where  fA,  fs, . . .,  are  die  wei^t  fractions  for  the  separate  elements  (A,  B, . . .)  present. 
(Attix,  1986;  Chilton,  1984;  Profio,  1979). 
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2.6.1  Examples  of  Mass  Attenuation  Coefllcients. 


Figure  2.6  shows  the  Compton,  photoelectric  and  Rayle^  mass  attenuation 
coefficients  for  two  common  elements,  caihon,  low  Z,  and  lead,  high  Z,  as  a  function  of 
photon  energy.  The  three  interaction  cross  sections  and  thdr  total  are  shown.  Figure  2.7 


shows  the  mass  coefficients  for  a  mixture,  gypsum,  present  in  diywaU. 


ho.MeV 


hv  (MeV) 


Figure  2.  6  Mass  attenuation  coefficients  for  carbon  and  lead. 

2.6.2  Beam  Attenuation. 

Once  die  total  cross  section  is  identified  for  a  material,  the  attenuation  of  the  beam 
through  that  material  can  be  determined.  Recall  that  the  cross  section  identifies  the 
probability  of  interaction  in  a  unit  diickness  of  material  traversed.  The  probability  that  it 
win  interact  in  an  infinitesimal  thickness  dl  is  \iJl.  If  N  particles  are  incident  upon  dl,  the 
change  in  the  number  N  due  to  absorption  is  given  by 


</N  =  -N 


ipdl). 


(2.19) 
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hf.  MeV 

Figure  2.  7  Mass  attenuation  coefficients  for  gypsum  (drywall). 

where  is  the  mass  attenuation  coefEicient.  hitegrating  over  the  tfiickness,  1,  from  0  to  L, 
and  corresponding  particle  populations  from  No  to  Nl  gives  the  law  of  exponential 
attenuation 

NL=Noe("^^\  (2.20) 

If  p.  is  a  total  attenuation  coefficient  it  is  singly  replaced  by  the  sum  of  individual 
coefficients,  a  +  t  +  Gr.  If  the  attenuation  coefficient  is  eneigy  dq)endent,  the  integration 
must  be  acconq)]ished  for  each  eneigy  interval.  The  effect  of  the  energy  dependent 
attenuation  coefficient  for  lead  (see  Fig.  2.6)  on  an  exponential  spectrum  is  shown  in 
Figure  2.8. 


2-13 


Figure  2.  8  Attenuation  of  an  exponential  energy  function  passing  through  1  cm  lead. 

2.7  Ionizing  Radiation  Quantities 

To  fiiify  understand  ionizing  radiation  a  de8ci:4)tion  and  definition  of  the  quantities 
and  units  typicalfy  used  is  necessary.  For  example,  flux  can  be  used  to  describe  the 
progress  of  the  photons  throu^out  the  problenu  Quantities  such  as  kerma  and  dose  relate 
photon  flux  to  the  effect  of  radiation  on  humans.  (Attix,  1986;  Profio,  1979). 

2.7.1  Fiuence. 

Fluence  is  the  mean  number  of  rays  striking  a  finite  sphere  surrounding  a  point,  P, 
during  an  arbitrary  time  interval.  Fluence,  dierefore,  is 


where  Ne  is  the  number  of  particles  and  a  the  great  circle  area  of  the  sphere.  It  is  expressed 
in  units  of  or  cm'^.  Fluence  expresses  the  sum  of  rays  incident  from  all  directions, 
irrespective  of  individual  quantum  or  kinetic  eneigies.  CMton’s  altemate  definition  of 
fluence  quoted  by  Atdx  is:  "fluence  at  a  point  P  is  numerically  equal  to  die  expectation 
value  of  the  sum  of  the  particle  track  lengths  that  occur  in  an  infinitesimal  volume  dV  at  P 
divided  by  dV.”  This  definition  eliminates  die  need  for  the  spherical  volume,  lending  itself 
to  dosimetiy  calculations  by  die  Monte  Carlo  method.  (Attix,  1986).  Chilton's  definition  of 


fluence,  when  multiplied  by  cross'section,  |x  (effectively,  collision  probability  per  path 


length),  ^es  collision  density: 


path  length  travelled  by  x  -  rays 
Volume 


'  collisions  ^ 

collisions 

^  cm  of  path  length; 

cm 

=  collision  density 


collisions 


L  cm' 


(2.22) 


2.7.2  Flux. 

At  any  point  during  the  above  aibitraiy  time,  a  flux,  or  fluence  rate,  is  the  change 
in  fluence  as  a  function  of  time 

do 

<p  =  — ,  (2.23) 

usually  expressed  in  units  or  m'  ^  s  or  cm  s  Fhix  expresses  die  sum  of  rays  incident 
from  all  directions,  irrespective  of  individual  quantum  or  kinetic  eneigies. 
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2.7.3  Energy  Fluence. 


The  energies  of  the  individual  rays  are  accounted  for  in  the  quantity  energy  ftuence 

dR 

=  — ,  (2.24) 

da 

where  R  is  tire  mean  total  energy  emitted,  transferred  or  received  (excluding  rest-mass 
energy)  within  the  volume.  Energy  ftuence  is  usually  expressed  in  units  of  J/m^  or  etg^cm^. 
For  the  case  of  a  sin^e  energy,  R  is  simpty  expressed  R  =  ENe  and  energy  ftuence  can  be 
related  to  particle  fhience,  Y  =  EO  .  The  energy  ftuence  rate  would  be  the  time 
differential  of  the  energy  ftuence 

d*? 

¥  =  (2.25) 

expressed  in  units  of  J  m'*  s‘^  or  erg  cm'*  s'*. 

2.7.4  Kerma. 

Kerma,  the  kinetic  mergy  released  in  material  is  mean  energy  transferred  to 
charged  particles  for  unit  mass  at  a  point  of  interest,  including  radiative-loss  energy  but 
excluding  energy  passed  from  one  charged  particle  to  another.  It  can  be  defined  in  terms  of 
energy  transferred,  8t„  and  radiant  energy,  R, 

6tr  =  (^m)u  -(^aut)r”'  +2:Q.  (2-26) 

where  (i?in)u  is  the  radiant  energy  of  fire  uncharged  particles  entering  V,  (Roiit)u"°”  is  the 

radiant  energy  of  uncharged  particles  leaving  V,  except  that  originating  from  radiative 

losses  of  kinetic  energy  by  charged  particles  while  in  V,  and  ZQ,  the  rest  mass  conversion 

term,  is  the  net  energy  derived  from  rest  mass  in  V.  Radiative  losses  are  the  conversion  of 
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charged-particle  kinetic  eneigy  to  photon  energy.  Energy  transferred  is  simpfy  the  kinetic 
eneigy  received  by  charged  particles  in  the  specified  finite  volume,  V.  Kerma,  therefore,  is 


<tetr 

dm 


(2.27) 


Average  kerma  is  the  mean  energy  divided  1^  the  mass 

Ka=^.  (2.28) 

n 

Kerma  and  average  kerma  can  be  expressed  in  units  of  eig/g,  rad,  J/kg  or  gray  (Gy)  simply 
related  by 

1  Gy  =  1  —  =  10^  rad  =  10^*  ^  .  (2.29) 

g 

The  kerma  fixnn  monoenergetic  photons  is  related  to  the  energy  fluence  by  the  mass 
energy-trrmsfer  coefficient,  (|Vp)e.z 


K='Fx 


V  P  - 


E.Z 


(2.30) 


2.7.5  Absorbed  Dose. 

The  absorbed  dose,  D,  the  mean  energy  imparted  to  matter  per  unit  mass  at  a 
point,  can  also  be  defined  in  terms  of  the  total  energy  imparted  by  ionizing  radiation  to 
matter  of  mass,  m,  in  finite  volume,  V, 


6  =  (i?in)u  -  iRcm)T  +  (^in)c  "  ”  +  2Q ,  (2-31) 
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where  the  sum  now  includes  the  uncharged  and  charged  radiant  energies  (see  Eq.  2.28) . 
The  absorbed  dose  is  therefore 

dc 

D  =  — .  (2.32) 

dm 

Dose  has  the  same  units  as  kerma,  J/lg,  eigi'g,  rads  or  gray  and,  like  the  average  kerma,  the 
average  absorbed  dose  is  the  energy  imparted  divided  by  the  mass  within  die  volume 


Da 


(2.33) 


2.7.6  Exposure. 

Exposure,  X,  is 


dm  ’ 


(2.34) 


where  dQ  is  ^e  absolute  value  of  die  total  charge  of  the  ions  of  one  sign  produced  in  air 
when  all  the  electrons  liberated  by  photons  in  air  of  mass  dm  are  completely  stopped  in  air” 
(ICRU,  1980).  Exposure  has  units  of  CVkg,  esu/g,  or  more  commonly,  roentgen  (R)  where 
1  R  =  2.58  X  lO"^  C/kg,  With  the  conversion  (in  air,  independent  of  photon  energy)  of 
33.97  J/C  for  the  mean  energy  expended  in  a  gas  per  ion  pair  formed,  1  R  =  8.78  x  10'® 

Gy  =  0.878  rad  (Attix,  1986;  Cember,  1983).  Exposure  rate  is  die  change  in  exposure 
over  time.  The  conv^ion  from  exposure  in  air  to  exposure  in  a  material  is  the  ratio  of  the 
mass  energy-absorption  coefficient  (Mw^p)nMt  for  the  material  to  that  of  air,  (Pat/p)w  which, 
for  muscle  tissue,  is  nearly  constant  over  the  diagnostic  energy  range,  i.e.  1.07  ±  3%. 
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2.7.7  Dose  Equivalent 


Dose  equivalent,  H,  (now  called  equivalent  dose  by  ICRP)  allows  the  average 
absorbed  dose  to  be  associated  with  a  risk  factor  '(^ch  is  an  estimate  of  the  relative  human 
hazard  from  different  ionizing  radiation  energies.  Effective  dose  equivalent,  £,  (now 
called  effective  dose  by  ICRP)  wei^ts  the  dose  equivalent  tissue  factors  to  associate  a 
risk  by  body  part  susceptibility: 


(2.35) 


H  =  X  D^; 

E  =  W.J.  X  H, 

where  Da  is  the  average  absorbed  dose  (section  2.8.5),  Wr  is  the  radiation  wei^ting  factor 


(also  denoted  by  Q)  and  Wt  is  the  tissue  wei^ting  factor.  The  units  of  dose  equivalent  and 
effective  dose  equivalent  are  J  kg'^  and  are  given  the  special  name  of  sievert,  Sv  (previously 
rem,  where  1  Sv  =  100  rem).  Tissue  and  radiation  weighting  factors  differ  from  source  to 
source.  10CFR20  and  NCRP  use  die  ICRP  recommendations.  Table  2. 1  contains  the 
recommended  die  tissue  wei^ting  factors  found  in  ICRP  Publication  60.  The  NRC’s 
recommended  tissue  wei^ting  factors  are  in  Table  2.2.  All  source  agree  on  a  radiation 
weighting  factor,  Wr  =  1  for  X  rays  of  all  energies.  (Scleien,  1993). 

Table  2.1  ICRP  Tissue  Weighting  Factors  (Scleien,  1993) 


1 _ Wi _ 

HHffIfiTHHi 

_ (U2 _ 

iniiMsiHHiii 

Skin 

Breast 

Limg 

^lilllllllllllH  1 

liver 

Red  bone  marrow 

Stomach 

Thyroid 

Totals 

0.02 

0.30 

0.48 

0.20 

'Remainder:  13  organs,  including  muscle,  sharing  equally  0.05,  Le.  2:0.004. 


2-19 


Table  2.2  NRC  Tissue  Weighting  Factors  (Scleien,  1993) 


Organ  or  Tissue 

Wt 

Gonads 

0.25 

Breast 

0.15 

Red  bone  marrow 

0.12 

Lung 

0.12 

Thyroid 

0.03 

Bone  surfaces 

0.03 

Remainder* 

0.30 

Whole  boch^** 

1.00 

“For  wei^ting  the  extonal  whole  body  dose  to  add  to  the  intonal  dose. 
0.30  results  from  0.06  from  each  of  5  remainder  organs  (excluding  skin 
and  lens  of  the  eye)  that  receive  the  hipest  dose. 


3.  Monte  Carlo  Method 


3.1  Introduction 

As  discussed  in  Chapter  2,  the  particle  cross  sections  can  be  considered  a 
probability  of  particle  interaction.  Given  this,  production  and  transport  of  ionizing 
radiation  can  be  treated  as  a  statistical  process  controlled  the  eneigy  of  die  impinging 

photons  and  the  atomic  number(s)  of  the  target  material.  Statistical  treatment,  and  thus, 
statistical  results,  is  well  suited  to  larger  scale  nuclear  anafysis  when  average  (statistical) 
behavior  is  the  goal  The  Monte  Carlo  method  of  simulating  individual  particles  and 
recording  some  aspects  of  thdr  average  behavior  is  one  such  statistical  tool. 

hi  the  Monte  Carlo  method,  average  behavior  of  the  particles  in  the  physical  system 
is  inferred  from  the  average  behavior  of  die  simulated  particles.  The  simulations  are 
performed  on  a  digital  conqiuter  because  the  numb^  of  trials  necessary  to  adequate^ 
describe  the  phenomenon  is  usually  quite  large.  The  statistical  sampling  process  is  based 
on  the  selection  of  random  numbers  -  analogous  to  thrown^  dice  in  a  gambling  casino  - 
hence  the  name  “Monte  Carlo”.  The  Monte  Cario  method,  in  short,  duplicates  a  process 
by  drawing  random  numbers  which  sample  a  probability  distribution  that  reflects  known 
behavior,  such  as  scattering  angle  or  track  length  distance  between  collisions.  In  particle 
transport  applications,  for  each  of  many  particle  histories,  a  particle  is  followed  from  a 
source  throughout  its  life  in  a  series  of  individual  events  to  its  death  in  some  terminal 
category  (absorption,  escape,  etc.).  Only  the  outcome  of  each  individual  step  is  determined 
by  the  random  draw. 
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3.2  Probability  Distribution  Functions 


Two  functions  essential  to  Monte  Carlo  calculations  are  the  probability^  density 
function  and  die  cumulative  probability  density  function.  The  probability  density  function, 
f(x)  is  the  probability  that  some  value  (in  diis  case,  the  random  draw,  x')  will  fall  between 
X  and  Ax.  It  is  necessary  that  f  (x)  2  0  and  normalized  over  its  restricted  interval,  [a,  b] 

f'*f(x)dx  =  l.  (3.1) 

Ja 

The  cumulative  probability  distribution  function,  F(x),  is  the  probability  that  the  random 
variable,  x’,  is  less  than  or  equal  to  x,  thus 


F(x)=  r  f(x’)dx'.  (3.2) 

J—cO 

It  is  often  convenient  to  write  £q.  (3.2)  in  the  form 


dF(x) 

dx 


=  f(x). 


(3.3) 


Figure  3.1  shows  the  relationship  between  f(x)  and  F(x). 

The  Monte  Carlo  method  relies  on  known  probability  functions.  To  use  these 
relationships,  rules  of  transformation  must  be  applied  to  the  randomly  drawn  variables. 

For  y  =  y(x),  let  g(y)dy  be  die  probability  that  y  is  between  y  and  y  +  dy,  and  f(x)dx  be  the 
probability  that  x  is  between  x  and  x  +  dx.  It  is  important  to  note  that  y'  (in  y  to  y  +  ^) 
and  x'  (in  X  to  X  +  dx)  both  describe  the  same  event  such  that  one  happens  if  and  only  if  the 
other  happens.  The  probability  density  functions  gty)  and  f(x)  must  satisfy 


|g(y)<*y|  =  |f(x)dx|, 


(3.4) 
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Figure  3.  1  Example  of  (a)  probability  density  function,  f(x),  and  (b)  cumulative 
probability  distribution  function,  F(x),  for  normalized  f(x). 


or,  since  g(y)  >  0  and  f(x)  >  0 


g(y)  =  f(x) 


dx 

dy 


(3.5) 


Now,  if  y  =  F(x)  of  the  random  variable  x,  Eq.  (3.5)  becomes 


g(F)  =  f(x) 


dx 


(3.6) 


Eq.  (3.3)  can  now  be  used  for  this  particular  transformation  to  obtain 


g(F)  =1,  0  <  F  <  1.  (3.7) 

Since  g(F)  is  the  probability  density  function  of  the  random  variable  F,  die  probability  of  F 

taking  on  a  value  between  F  and  F  +  dF  is  just  equal  to  dF  and  thus  is  uniformly  distributed 

between  zero  and  one.  The  random  numbers,  available  on  most  computers  via  library 

routines  or  built-in  language  functions  are  also  sanqiled  between  zero  and  one,  thus  F(x) 

may  be  sanqiled  uniformly 


F(x)=^ 


(3.8) 
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However,  we  need  to  sample  the  distribution  of  x  and  not  F(x),  thus,  after  each  call,  we 
must  perform  the  inversion 

x  =  F-^(^).  (3.9) 

(Lewis,  1984). 

3.3  Random  Numbers 

Random  numbers  are  generated  by  a  computer  system  with  complicated  numerical 
algoriduns.  The  use  of  the  algorithm  means  that  the  numbers  are  not  trufy  random,  but,  in 
fact,  pseudo  random,  i.e.,  their  sequence  can  be  duplicated  and  there  is  only  a  finite 
number  of  unique  “random”  numbers.  Processes  such  as  the  Monte  Carlo  method  can  use 
an  extraordinary  number  of  random  numbers  -  as  many  as  tens  of  thousands  per  particle 
history.  A  sample  Monte  Carlo  run  performed  for  this  analysis  used  300  million  random 
numbers  to  process  two  million  particles,  widi  a  maximum  of 28,000  for  a  single  history. 

It  is  important,  dierefore,  to  make  efficient  use  of  the  limited  number  of  random  numbers 
available.  This  is  accomplished  through  selection  or  rejection  schemes.  The  choice  of  the 
scheme  is  driven  by  the  nature  of  the  probability  density  ftmction.  A  simple,  intuitive 
method  of  samplir^  randomfy  jfrom  a  given  distribution  ftmction  uses  prepared  look-up 
tables.  Anafyticalfy  integrable  ftmctions  can  use  the  more  mathematically  attractive  and 
efficient  inverse  cumulative  ftmction  method  (also  called  the  selection  method).  Those  that 
are  not  anafyticalfy  integrable  must  use  the  rejection  method  (with  possible  enhancements 
to  improve  efficiency). 


3-4 


3.3.1  Selection. 


To  illustrate  selection  methods  of  random  numbers,  consider  the  case  of  choosii^  a 
random  distribution  over  a  circle  as  shown  in  Figure  3.2.  We  could  choose  random  x’s  and 


Figure  3.  2  Choosing  a  uniformly  distributed  random  point  within  a  circle, 
y’s  over  the  interval  [-a,  a].  This,  however  would  fill  the  entire  square  and  thus  waste 

those  random  numbers  falling  outside  the  circle.  On  the  odter  hand,  we  could  choose 

random  r’s  over  the  interval  [0,r]  and  0’s  over  the  interval  [0,27i].  This,  however,  woidd 

cause  “clumping”  of  points  in  the  center  of  the  circle,  i.e.  a  non-uniform  distribution. 

Instead,  let  dA,  the  differential  area  of  die  circle,  be  equalfy  likely  vriiere  dA  =  r  dr  d6 

Using  the  mediod  starting  witii  Eq.  (3.6) 


g(e)  =  f(x) 


dx 

Hb 


=  rdr 


d0 


=  Ca  (constant). 


(3.10) 


Since  we  know  6  must  be  distributed  around  the  entire  circle,  we  let  C  =  l/(2n)  to 


normalize  g(6)  to  one.  Then 


1  d© 


(3.11) 
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where  ©  is  flie  ciimulative  probability  function  over  6. 
Using  Eq.  (3.8): 

271  27t 

and  performing  the  inversion 


©-'(?)= e=274j. 

Following  the  same  steps  for  r: 


g(r)  =  f(x) 


dx 

T 


=  de 


rdr 

dr 


=  rC. 


To  normalize 


(3.12) 


(3.13) 


(3.14) 


Jg(r)dr  =  JrCj.dr  =  1  such  that  Cj.=' 


.2  ■ 


(3.15) 


then 


(3.16)  and 


(3.17) 


Thus,  for  every  point  within  the  circle,  two  random  numbers  will  be  drawn,  and 


which,  via  Eqs.  (3.13)  and  (3.17),  will  become  the  (r,  0)  pair  describing  the  location  of  the 
point.  The  selection  method  will  efficiently  produce  uniformly  distributed  points 
throughout  the  circle  without  the  need  to  reject  any  numbers.  This  method,  however,  onh^ 
works  for  anafyticalfy  integrable  functions.  Non-integrable  functions  may  be  numerically 
integrated  to  get  F(x)  followed  by  table  mteipolation  to  evaluate  F‘(^).  (Mathews,  1995). 
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3.3.2  Rejection. 


The  rejection  method  is  used  for  non-in^rable  functions  and  is  not  as  efficient  as 
the  selection  method  described  in  section  3.3.2.  In  implementing  the  rejection  method,  the 
function  is  not  integrated  (relating  explicit^  a  given  x  to  some  y  via  the  function),  so  two 
random  munbers  must  be  drawn  to  define  a  single  point.  The  position  of  the  point  relative 
to  die  function  determines  acceptance/rejection.  The  first  random  draw  determines  Xi  on 
the  interval  [xi,  x^]  (see  upper  section  of  Figure  3.3).  The  next  draw  will  choose  a  point, 
yi,  on  the  Xi  line.  If  yi  >  f  (?^),  it  will  be  rejected;  if  yt  >  f(Xi),  it  is  accepted.  (Note:  the 
ratio  of  accepted  to  rejected  numbers  is  the  volume  under  the  curve,  Le.  integration..) 

The  lower  portion  of  Figure  3.3  illustrates  an  enhancement  to  the  rejection  method 
that  may  be  used  to  improve  efficiency.  The  function  is  broken  into  a  number  piece-wise 
int^able  functions.  They  are  expanded  to  fill  y.  i^ropiiate  wei^ts  are  assigned  to  each 
interval 

wi  =  r"fi(x)dx 

JXi 

^2  =  Jx2^2(x)dx  (3.  18) 

W3  =  r''f3(x)dx 

•'*3 

If,  for  example,  in  Figure  3.  3,  <  Wi,  selection/rejection  is  determined  witii  fi(x);  if  Xi  < 

W2,  with  f2(x);  etc.  In  this  example,  the  ratio  of  the  rejection  areas  to  the  acceptance  areas 
are  seen  to  be  smaller  than  before,  thus,  fewer  numbers  wiU  be  discarded. 
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3.3.3  Discrete  Sampling. 


When,  instead  of  assigning  a  value  such  as  distance  to  collision,  the  random 
number  is  used  to  answer  a  question,  discrete  sampling  is  used.  For  example,  we  may  need 
to  know  if  the  next  event  will  be  a  capture,  elastic  scatter,  etc.  Recalling  the  cross-sections 
of  Chapter  2,  the  total  cross-section  is 

Oj  =  Oi  +  CT2  +  O3  +  O4  +  •••  (3.19) 


Xj  Xj  X3  X4 


Figure  3.  3  Rejection  method  shown  in  upper  portion.  Lower  portion  illustrates  an 
enhancement  to  the  rejection  method. 
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The  unit  line  representing  the  range  of  random  numbers  (0-1)  is  subdivided  into  intervals 
proportional  to  (oi  /  a-r) ,  (02  /  ai) ,  etc.  The  random  number  falls  into  one  of  the  intervals 
which  correspond  to  tiie  "answer"  sought.  (Lewis,  1984) 

3.4  Errors 

The  final  answer  obtained  in  Monte  Carlo  analysis  is  the  average  of  the 
contributions  fi-om  many  histories.  Equally  important  to  the  answer  is  the  statistical  error 
or  uncertainty^  associated  with  the  answer.  Not  only  the  absolute  error,  but  its  behavior  as 
you  increment  particle  histories,  must  be  anafyzed  to  ensure  confidence  in  the  results. 


3.4.1  Means,  Variances  and  Standard  Deviations. 


The  random  nature  of  a  Monte  Carlo  anatysis  typicalty  produces  a  range  of  scores 
depending  on  file  area  of  interest  and  the  variance  reduction  methods  (see  section  3.5 
Variance  Reduction)  chosen.  As  seen  above,  the  expectation  value  of  x  subject  to  f(x),  the 
score  probability  density  function  is 

E(x)  =  ^x)=  Jxf(x)dx  =  true  mean.  (3.20) 

Since  f(x)  is  seldom  explicitfy  known,  it  is  implicitty  sampled  by  the  Monte  Carlo  random 
walk  process.  The  sample  mean  is  then  an  estimate  of  the  true  mean,  x 


(3.21) 


where  x^  is  file  value  of  x  selected  firom  f(x)  for  the  i***  history  and  N  is  the  number  of 
histories  calculated  in  the  problem.  The  Strong  Law  of  Large  Numbers  states  that  if  £(x) 
is  finite,  x  tends  to  the  limit  E(x)  as  N  approaches  infinity  thus  providing  a  relationship 


between  E(x)  and  x .  The  variance  of  the  population  of  x  values  is  a  measure  of  the 
difference  in  the  two  values  and  is  given  by 

=  J(x-E(x)f  f(x)dx  =  E(x2)-(E(x))2.  (3.22) 

The  square  root  of  tiie  variance,  a,  is  called  the  standard  deviation  of  the  population.  As 
with  E(x),  a  is  seldom  known  but  can  be  estimated  by  Monte  Carlo  as  S,  v^ch  for  large 
numbers  is 


^^2 


N-1 


2  -2 
=  X  -X  . 


(3.23) 


The  quantity  S  is  the  estimated  standard  deviation  of  the  population  of  x  based  on  die 
values  of  Xi  that  were  actualty  sampled.  The  estimated  variance  of  x  is 

o2 

S|  =  —  (3.24) 

*  N 

It  is  important  to  note  that  S  x  is  proportional  to  1/N^^.  This  is  a  drawback  because 
four  times  the  number  of  histories  must  be  run  to  halve  S  x  which  can  be  computationally 
expensive.  Another  way  to  reduce  estimated  variance  is  by  reducing  S,  the  estimated 
standard  deviation,  using  variance  reduction  techniques  (see  Section  3.5). 


3.4.2  Estimated  Relative  Errors. 

A  convenient  representation  of  error  is  the  estimated  relative  error,  R 


(3.25) 
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It  represents  statistical  precision  as  a  fractional  result  with  respect  to  the  estimated  mean. 
Based  on  the  experience  of  Monte  Carlo  practitioners  and  years  of  experience  using 
MCNP  on  wide  variety  of  problems,  Table  3.1  presents  the  recommended  interpretation 
of  the  estimated  la  confidence  intervals  (Biiesmeister,  1993). 

Table  3.  1  Guidelines  for  Interpreting  the  Relative  Error,  R 


Range  of  R 

Quality  of  the  Talty 

0.5  to  1 

Useless 

0.2  to  0.5 

limited  usefiilness 

0.1  to  0.2 

Questionable 

<0.10 

Grenerally  reliable 

3.4.3  Figure  of  Merit 

By  definition  should  be  proportional  to  1/N  and  the  computer  time,  T, 
proportional  to  N,  thus,  1/(R^),  the  Monte  Carlo  figure  of  merit  (FOM),  should  be 
approximate^  constant  within  any  one  run  (with  the  exception  of  statistical  fluctuations 
early  in  the  problem.)  It  is  advantageous  to  have  a  large  FOM  since  that  would  indicate  an 
efficient  run  yielding  a  small  relative  error  from  short  computer  run  times. 


3.5  Variance  Reduction 

Variance  reduction,  as  discussed  above,  is  a  method  to  reduce  errors  efficiency. 
There  are  three  basic  types;  splitting,  Russian  roulette  and  sampling  from  nonanalog 
probability  density  functions.  Critical  to  all  three  techniques  is  the  process  of  weighting, 
which  refers  to  the  relative  contribution  of  a  particle  to  the  fiiuil  results.  Each  particle  is 
multiplied  by  a  weight,  w,  so  that  the  full  results  of  the  w  physical  particles  represented  by 
each  particle  modeled  are  exhibited  in  the  final  results.  Specific  variance  reduction 
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techniques  will  be  discussed  in  Chapter  4  with  the  MCNP  description.  (Briesmeister, 
1993). 


3.5.1  Splitting. 

Splitting  refers  to  dividing  the  particle’s  weigjht  among  two  or  more  daughter 
particles  and  following  them  independently.  Usualfy,  the  weight  is  divided  evenly  among  k 
identical  daughter  particles  whose  characteristics  are  identical  to  the  parent  except  for  the 
factor  \!km  weight.  The  weight,  however,  may  be  unevenfy  divided  to  represent  tiie 
e?q)ected  number  of  physical  particles  that  would  select  outcome  j  from  a  set  of  A:  mutually 
exclusive  outcomes. 

3.5.2  Russian  Roulette. 

Russian  roulette  takes  a  particle  of  weight  wo,  turns  it  into  a  particle  of  weight  > 
Wo  with  probability  Wq/wi  and  kills  it  (weight  =  0)  with  probability  ( 1  -  Wq/wi  ).  The 
expected  wei^t  is  Wi  ♦  wq/wi  +  ( 1  -  wq/wi  )  *  0  =  wq  so  tiie  metiiod  doesn't  bias  tiie 
results  as  will  also  be  seen  in  tiie  analog  game. 

3.5.3  Sampling  From  Nonanalog  Probability  Density  Functions. 

An  analog  Monte  Carlo  calculation  samples  the  events  according  to  their  natural 
physical  probabilities,  thus  estimating  the  number  of  particles  executing  any  given  random 
walk.  Nonanalog  techniques  are  free  to  do  anything  that  will  preserve  the  correct  particle 
density  while,  hopefully,  reducing  the  variance.  This  is  accomplished  by  adjusting  the 
wei^ts  so  the  outcome  matches  that  of  the  analog  game. 
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4.  Monte  Carlo  N-Particle  Transport  Code  System  (MCNP) 


4.1  Introduction 

The  Monte  Carlo  N-Particle  Transport  Code  System  (MCNP,  formerly  Monte 
Cario  Neutron-Photon)  is  a  general-purpose,  continuous-energy,  generalized-geometty, 
time-dependent,  coupled  neutron/photon  Monte  Carlo  transport  code.  It  was  first 
developed  by  Los  Alamos  National  Laboratory  (LANL)  using  the  Monte  Carlo  methods 
that  emerged  from  work  done  there  during  World  War  II  by  Fermi,  von  Neumann,  Ulam, 
Metropolis  and  Richtmeyer.  Version  I,  a  compilation  of  earfy  work  with  neutron  and 
photon  transport,  appeared  in  1973.  It’s  present  version,  4a,  has  over  350  man-years 
invested  in  code  production,  testing  and  optimization  by  the  Monte  Carlo  section  of  die 
radiation  transport  group  (X-6)  of  the  applied  theoretical  physics  division  (Xdivision)  of  the 
LANL.  Its  single  source  code  complies  with  the  FORTRAN  77  standard  and  thus  can  be 
run  on  many  different  types  of  computers:  Cray,  CDC,  IBM,  VAX,  PRIME,  SIEMENS, 
RIDGE,  APOLLO,  PC’s  and  SUN.  MCNP  is  able  to  analyze  neutron-only,  photon-only 
or  combined  neutron/photim  transport.  A  user  created  input  file  is  subsequently  read  by 
MCNP.  This  file  contains  information  about  the  problem  to  be  modeled:  the  evaluations 
to  use,  the  location  and  characteristics  of  the  neutron  or  photon  source,  the  geometry  of  the 
problem,  the  type  of  answers  or  tallies  desired  and  any  variance-reduction  techniques  used 
to  make  the  problem  run  more  efficiently.  (Briesmeister,  1993). 


4-1 


4.2  Method 


The  Monte  Cailo  method  solves  a  transport  problem  by  simulating  particle  histories 
that  are  transported  between  events  (for  example,  collisions)  that  are  separated  in  space  and 
time.  The  individual  probabilistic  events  that  comprise  a  process  are  simulated  sequential^ 
and  the  probability  distributions  governing  these  events  are  statistically  sampled  to  describe 
the  total  phenomenon.  The  sampling  process  is  based  on  the  selection  of  random  numbers. 
Each  particle  is  followed  from  a  source,  through  events,  to  its  death  in  some  terminal 
category  (absorption,  escape,  etc.). 

MCNP  uses  continuous-energy  nuclear  data  libraries.  There  are  over  500  data 
tables  for  neutron  interaction,  photon  interaction,  and  neutron  dosimetry  or  activation  for 
approximate^  100  different  isotopes  or  elements.  The  data  tables  are  sufficiently  dense  on 
the  energy  grid  that  linear-linear  interpolation  induces  less  than  1%  deviation  from  the 
evaluated  cross  sections.  Cross  sections  are  available  for  nearly  2000  dosimetry  or 
activation  reactions  involving  over  400  target  nuclei  in  ground  and  excited  states.  The  best 
available  evaluation  is  clearly  marked  on  the  tables  of  nuclides  included  with  the  model. 
(Briesmeister,  1993) 

4.2.1  Source. 

MCNP’s  generalized  user-input  source  capability  allows  the  irser  to  specify  a  wide 
variety  of  source  conditions  without  requiring  code  modifications.  Independent  probability 
distributions  nuiy  be  specified  for  the  source  variables  (energy,  time,  position,  and 
directions)  and  other  parameters  (such  as  starting  cell  or  surface).  Built-in  distribution 
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functions  are  available  for  various  energy  spectra  such  as  Watt,  Maxwellian  and  Gaussian; 
Gaussian  for  time;  and  isotropic,  cosine  and  mono-directional  for  direction. 

4.2.2  Tallies. 

MCMP  provides  a  measured  quantity  called  a  talfy.  A  variety  of  particle  flux  and 
particle  cuirent  tallies  are  obtainable:  (1)  current  as  a  fimction  of  direction  across  any  set 
of  surfaces,  surface  segments  or  sum  of  surfaces  and  (2)  flux  across  any  set  of  surfaces, 
surface  segments,  sum  of  surfaces  and  within  cells,  cell  segments  or  sum  of  cells.  All 
tallies  are  normalized  per  starting  particle. 

4.2.2.1  Flux  Tames. 

Standard  flux  tallies  (where  the  term  "flux”  is  used  flu-oughout  the  model 
documentation  and  output  files  regardless  of  the  actual  measured  quantity)  are  obtained 
with  point  or  ring  detectors  (used  exclusivety  in  this  analysis).  They  are  estimates  of 

Fj  =  J|<|)(f,E,t)dEdt,  (4.1) 

t  E 

where  (|),  particle  flux  density,  is  the  number  of  particles  at  distance  r  and  time  t  with 
eneigjes  between  E  and  E  +  AE.  The  units  of  the  flux  tally  are  the  units  of  the  source.  If 
the  source  has  units  of  particles  per  unit  time,  the  talty  is  also  particles  per  unit  time.  When 
the  source  has  units  of  particles,  this  tally  represents  a  fluence  tally.  Flux  can  be  obtained 
from  the  fluence  talty  for  a  time-dependent  source  by  dividing  the  talty  by  flie  time  interval 
assuming  the  source  is  constant  over  time. 
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The  flux  tally  is  obtained  at  user  defined  point  detector  locations.  A  point  detector 
is  a  stochastic  estimate  (fi'om  each  event  location)  of  die  flux  at  a  point  in  space.  It  is 
known  as  a  “next  event  estimator”  because  it  is  a  tally  of  Ifae  flux  at  a  point  if  the  next  event 
is  a  trajectory  without  further  collision  directly  to  die  point  detector.  Contributions  to  the 
point  detector  tally  are  made  at  every  event  throu^out  the  random  walk.  Each 
contribution  to  the  flux  talfy  can,  therefore,  be  thou^t  of  as  the  transport  of  a  pseudo- 
particle  to  the  detector.  The  final  quantity  reported  in  MCNP  (before  the  final 
normalization  per  starting  particle)  for  the  flux  tallies  is  the  sum  of  the  scores  fi'om  each  of 
N  individual  events.  Each  score  is  a  function  of  the  particle  weight,  w,  its  distance,  R,  to 
the  detector,  the  cosine  of  the  angle  between  surface  normal  and  particle  trajectory,  p,  and 
the  total  mean  fi-ee  path  to  the  detector,  X: 

tally  =  2;wj  X p(n,) X 2EL^.  (4.2) 

iti  2nRi* 

The  exp(-X)  term  accounts  for  the  attenuation  between  the  present  event  and  the  detector 
point.  The  p(p)  term  accounts  for  the  probability  of  scattering  toward  the  detector  (instead 
of  the  direction  selected  by  die  random  walk).  The  l/(27iR^)  term  accounts  for  the  solid 
angle,  dQ,  die  particle  must  scatter  into  to  arrive  at  the  detector.  The  R^  term  in  the 
denominator  can  cause  a  singularity  and  thus  infinite  flux  if  a  source  or  collision  event 
occurs  near  the  detector  point.  For  diat  reason,  a  fictitious  sphere  is  placed  around  the 
detector  becoming  an  average  flux  region,  i.e.  individual  contributions  within  the  sphere  are 
not  counted,  rather,  the  average  flux  is  assumed  within  die  sphere. 
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A  ring  detector  tally  is,  in  essence,  a  point  detector  tally  widi  the  point  detector 
location  not  fixed  but,  instead,  sampled  fi'om  some  location  on  a  ring.  It  can  be  thought  of 
as  a  floating  point  detector  that  moves  to  flie  location  on  the  ring  nearest  an  event.  The 
tally  is  averaged  over  the  circumference  of  the  ring.  Ring  detectors  can  be  much  more 
efficient  than  point  detectors  but  their  use  is  only  valid  when  the  problem  geometty  is 
axially  symmetric. 

Point  and  ring  detectors  are  subject  to  dramatic  Jumps  in  the  talfy  fi’om  high  weight 
particles.  It  is  therefore  vital  that  a  sufficient  number  of  particles  be  made  to  enter  a  region 
or  volume  of  interest  so  it  may  be  thorou^ify  sampled  using,  for  example,  a  combination 
of  importance  biasing  and  geometry  splitting,  typical  variance  reduction  techniques.  Tally 
fluctuation  charts  contain  a  cumulative  history  of  the  tally  as  the  run  progresses  and  must 
be  examined  to  confirm  convergence. 

4.2.2.2  Errors. 

Associated  with  each  tally  in  the  ouq)ut  file  is  the  estimated  relative  error,  R, 
corresponding  to  one  standard  deviation  of  flie  merui  divided  by  the  mean,  x,  which, 

for  a  well-behaved  tally,  will  be  proportional  to  where  N  is  the  number  of  histories. 

This  variance  refers  only  to  the  precision  of  the  Monte  Carlo  calculation  itself  and  not  to 
the  accuracy  of  the  result  compared  to  the  true  plysical  value.  MCNP  also  reports  a 
figure  of  merit,  FOM,  for  one  tally  bin  of  each  tally  as  a  function  of  the  number  of 
histories  where  FOM  =  1/(R*T)  (T,  computer  time  in  minutes).  The  more  efficient  the 
Monte  Carlo  run,  die  larger  the  FOM. 
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4.2.3  Variance  Reduction. 


The  error  in  a  Monte  Carlo  result  is  directly  related  to  the  number  of  particles 
simulated.  One  could,  of  course,  run  a  large  number  of  particles  until  acceptable  error  is 
obtained.  There  is  value,  however,  in  reducing  the  errors  by  less  costly  means,  i.e. 
variance  reduction  techniques.  The  strength  of  MCNP  lies  in  the  large  variety  of  variance 
reduction  techniques  available. 

In  the  simplest  Monte  Carlo  model,  particles  are  followed  from  event  to  event  with 
the  next  event  sampled  from  a  number  of  possibilities  according  to  natural  event 
probabilities.  This  analog  model  works  well  when  a  significant  fraction  of  the  particles 
contribute  to  each  estimate.  However,  many  cases  arise  where  the  fi:uction  of  particles 
detected  is  very  small:  the  Monte  Carlo  model  then  fails  with  unaccepable  statistical 
imcertainties. 

Although  the  analog  model  is  simple,  there  are  other  probability  models  for 
transport  that  will  estimate  the  same  average  value  witii  a  much  smaller  variance.  A  non¬ 
analog  model  attempts  to  follow  interesting  particles  more  than  uninteresting  ones,  vvdiere 
the  interesting  particle  is  one  tiiat  contributes  a  large  amount  to  the  quantities  that  need  to 
be  estimated.  To  ensure  the  average  score  remains  the  same,  the  score  is  modified  to 
remove  the  effect  of  biasing  the  natural  odds.  Thus,  if  a  particle  is  artificialfy  made  q  times 
as  likely  to  execute  a  given  random  walk,  its  score  is  weighted  by  Hq,  preserving  the 
average  score. 

Any  scheme  tiiat  increases  the  number  of  particles  available  for  scoring  will 
decrease  the  error.  Related  to  an  increase  in  particles  is  use  of  the  particles  in  an  efficient 
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manner.  An  ejfficient  scheme  would  remove  particles  not  likety  to  contribute  significantly 
to  the  ta%,  allowing  more  total  particles  to  be  run  within  the  same  time.  A  problem  with 
some  variance  reduction  techniques  is  that  diey  often  are  tailored  to  aid  one  tally  at  the 
expense  of  all  others.  This  forces  separate  runs  for  each  optimized  tally.  The  following 
non-analog  variance  reduction  schemes  operate  to  increase  the  number  of  particles, 
increase  the  efGciency  with  which  they  are  used  or  both: 

1.  Importance  Sampling:  Regions  in  the  problem  geometry  are  assigned  relative 
importances  so  particles  transported  from  a  region  of  higher  importance  to  a  region 
of  lower  importance  can  be  Russian  rouletted,  that  is,  killed  a  fraction  of  the  time 
according  to  the  ratio  of  die  r^ons’  importances.  The  survivors  will  be  counted 
more  by  increasing  their  weight.  Particles  transported  into  a  region  of  higher 
importance  will  be  split  into  two  or  more  particles  with  lower  weight. 

2.  Time  and  Energy  Cutoff:  Particles  with  a  time  or  energy  out  of  the  range  of 
interest  are  killed. 

3.  Implicit  Capture  vs.  Analog  capture:  A  particle  has  a  probability  of  capture  when 
it  collides  with  a  nucleus.  Analog  c^ture  kills  the  particle  with  that  probability.  In 
implicit  capture  (survival  biasing),  the  particle’s  wdght  is  reduced  by  the  capture 
probability,  thus  allowing  important  particles  to  survive  the  coUisional  capture. 

4.  Exponential  Transformation:  To  transport  particles  long  distances,  die  distance 
between  collisions  in  the  preferred  direction  is  artificially  increased  and  the  weight 
is  decreased.  (Used  with  wei^t  windows  because  of  the  large  weight  fluctuations.) 

5.  Forced  Collisions:  When  a  particle  enters  an  almost  transparent  cell  (low  cross- 
section),  it  may  be  forced  to  undergo  a  collision  by  splitting  it  and  its  weight  into 
collided  and  uncollided  parts. 

6.  Energy  Splitting/Russian  Roulette:  Particles  entering  important  and  less  important 
energy  ranges  can  be  split  or  Russian  rouletted,  respectively. 

7.  Correlated  Sampling:  The  i“*  history  will  always  start  at  the  same  point  in  the 
random  number  sequence  regardless  of  the  (i-l)'*'  particle’s  history  so  statistical 
fluctuations  in  the  solution  will  not  mask  small  variations  caused  by  the  problem 
specification.  In  particular,  it  allows  two  possible  slightly  differing  calculations  to 
be  compared  independent  of  differences  caused  by  random  number  generation. 

8.  Source  Biasing:  Source  particles  of  higher  importance  are  emitted  with  a  hi^er 
fiiequency  but  with  a  compensating  lower  weight  than  less  important  source 
particles. 

9.  Point  Detectors:  The  probability  of  transporting  a  particle  precisely  to  a  point 
detector  is  vanishingh^  small.  Therefore,  pseudo-particles  are  directed  to  the  point 
instead.  With  eveiy  new  particle  from  the  source  and  at  each  event  or  collision,  the 
pseudo-particle’s  probability  of  reaching  the  detector’s  location  is  used  as  its  weight 
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with  tile  weighted  pseudo-particles  summed  to  obtain  the  flux  at  the  point  detector. 
Ring  detectors  are  essential^  circular  point  detectors  over  which  tiie  flux  is 
averaged.  They  are  only  available  for  problems  with  axial  symmetry. 

10.  DXTRAN:  When  the  problem  contains  a  small  regimi  of  interest  that  otherwise 
would  be  difiQcult  to  sample,  a  user-specified  “DXTRAN  sphere”  can  be  used  that 
will  sum  particles  weighted  by  their  probabilily  of  entering  the  sphere  with  no  other 
collisions.  Any  particle  that  tries  to  enter  the  sphere  will  be  killed  since  it  has 
alrea^  been  accoimted  for  deterministically. 

11.  Weight  Window:  To  keep  tire  particle  weight  distribution  within  reasonable 
bounds,  low-weighted  particles  are  eliminated  by  Russian  roulette  and  high 
weighted  particles  are  split  as  a  function  of  energy,  geometric  location  or  both. 
MCNP  can  provide  nearfy  appropriate  weight  windows  (determined  for  onfy  one 
detector  at  a  time)  throu^  its  weight  window  generator  function. 

4.2.4  Geometry. 

The  geometry  utility  of  MCNP  treats  an  arbitrary  three-dimensional  configuration 
of  user-defined  materials  in  geometric  ceUs  bounded  by  first-  and  second-  degree  surfaces 
and  some  special  fourth-degree  surfaces  (elliptical  tori).  Cells  are  defined  by  intersections, 
unions  and  complements  of  the  surfaces  bounding  the  desired  regions.  Surfaces  are 
defined  by  supplying  coefficients  to  the  Cartesian  coordinate  analytic  surface  equations  or 
for  certain  types  of  surfaces,  by  supphmg  known  points  on  the  surfaces.  The  code  does 
extensive  internal  checking  to  find  input  errors,  hi  addition,  an  elaborate  geometiy-plotting 
capabili^  allows  the  user  to  view,  and  tiius,  define  and  check  the  geometry. 


4.2.5  Cross  Sections. 

As  discussed  in  Chapter  2,  transport  problems  require  particle  interaction  cross 
sections  for  each  material.  The  producers  of  MCNP  have  taken  the  responsibility  to  ensure 
that  the  most  current  and  accurate  nuclear  data  tables  are  available.  Tables  are  provided 
that  cover  all  the  different  interactions  modeled  by  MCNP. 
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Each  nuclide  library  contains  an  energy  grid  including  the  photoelectric  edges  and 
pair  production  thresholds.  These  are  followed  by  tables  of  coherent  and  incoherent  form 
factors  tabulated  as  a  function  of  momentum  transfer.  Next  are  incoherent  scattering, 
coherent  scattering,  photoelectric  and  pair  production  cross  sections.  The  total  cross 
section  is  a  summation  of  the  other  cross  sections. 

Angular  distributions  of  secondary  photons  are  isotropic  for  photoelectric  effect, 
fluorescence  and  pair  production  and  are  obtained  by  sampling  the  well-known  Thomson 
and  Klein-Nishina  formulas  for  coherent  wd  incoherent  scattering,  respectively.  The 
energy  of  incoherent  scattered  photons  is  calculated  ffom  die  sampled  scattering  angles, 
using  the  Compton  relation. 

4.3  Physics 

MCNP  treats  photon  interactions  the  same  as  neutrons  for  sampling  of  a  collision 
nuclide,  analog  capture,  implicit  capture  and  variance  reduction.  The  collision  physics, 
however,  are  completely  different.  (Briesmeister,  1993). 

4.3.1  Sampling  of  a  Collision  Nuclide. 

When  a  photon  interacts  within  a  material  composed  of  more  than  one  element 
(compound,  mixture,  etc.)  the  particle  ^e  it  collides  with  is  determined  by  the  relative 
proportions  of  the  individual  cross  sections  within  die  total  (see  Eqs.  2-19  and  3-19) 
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4.3.2  Capture:  Analog  and  Implicit 


Analog  capture  kills  particles  with  probability  a,  /  Ot  where  a,  and  Ot  are  the 
absorption  and  total  cross  sections  of  the  collision  nuclide  at  the  incoming  photon  eneigy. 
All  particles  killed  by  analog  capture  have  their  entire  particle  energy  and  wei^t  deposited 
in  the  collision  cell. 

hnplicit  capture  reduces  the  photon’s  weight  to  wo’ 

«-„•=(! (4.2) 

G  'j» 

If  the  new  weight  is  below  the  problem  weight  cutoft  the  photon  is  rouletted, 

resulting  in  fewer  particles  with  larger  wei^ts.  A  fraction,  a,  /  Gt,  of  the  incident 
particle’s  wei^t  and  eneigy  are  deposited  in  die  collision  cell  corresponding  to  the  portion 
of  the  particle  captured.  Implicit  capture  is  the  default  method  of  capture  in  MCNP. 

4.3.3  Collisions. 

MCNP  has  two  interaction  modes  for  modeling  collisions,  simple  and  detailed. 

The  simple  treatment  is  characterized  by  its  neglect  of  coherent  scattering  and  fluorescent 
photons  from  photoelectric  absorption  and  is  intended  for  high  eneigy  photons  (model 
default:  lOOMeV).  The  detailed  treatment  includes  coherent  scattering  and  accounts  for 
fluorescent  photons  afrer  photoelectric  absorption.  The  detailed  treatment  is  described  in 
the  following  sections. 
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4.3.3. 1  Incoherent  (Con^ton)  Scattering. 

MCNP  detennines  the  angle  0  (scattering  angle  relative  to  the  incident  line 
of  fli^t),  hv’,  the  new  energy  of  the  photon  and  hv  -  hv’,  the  recoil  kinetic  energy  of  the 
electron.  MCNP  uses  the  Klein-Nishina  cross  section  modified  by  an  appropriate 
scattering  factor  to  determine  6.  This  scattering  factor  acts  to  decrease  the  Klein-Nishina 
cross  section  (per  electron)  in  the  forward  directions  for  low  hv  and  high  Z,  independent^. 

4.3.3.2  Coherent  (Thomson  or  Ri^leigh)  Scattering. 

For  coherent  scatter  events,  on^  the  scattering  angle  0  is  computed  after 
which  transport  of  the  photon  continues.  MCNP  modifies  the  energy  independent 
Thomson  cross  section  with  a  form  factor  that  acts  to  decrease  the  cross  section  for 
backward  scattering  of  high  hv,  low  Z  interactions.  This  strongly  peaks  the  cross  section  in 
the  forward  direction  to  the  extent  that  at  high  energies,  coherent  scattering  can  be  ignored. 

4.3.3.3  Photoelectric  Effect 

The  phototelectric  effect  consists  of  the  absorption  of  the  incident  photon 
(with  energy  £)  and  the  subsequent  emission  of  fluorescent  photons  and  the  ejection  or 
excitation  of  an  orbital  electron  (binding  energy,  e  <  £).  Three  cases  of  the  photoelectric 
effect  are  treated:  the  emission  of  zero,  one  and  two  fluorescent  photons  with  energy 
greater  than  1  keV  after  the  ejection  or  excitation  of  an  orbital  electron.  In  theirs/  case, 
since  no  photons  are  emitted,  the  photon  track  is  terminated  and  scored  like  an  analog 
capture  in  the  surrunary  table  of  the  output  file.  The  cascade  of  electrons  that  fills  up  the 
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orbital  vacancy  left  by  the  ejection  (Auger  effect)  produces  electrons  and  low-eneigy 
photons  within  the  model,  hi  the  second  case,  one  photon  is  emitted  with  an  energy  equal 
to  die  incident  photon  energy  less  the  ejected  electron  kinetic  energy  less  a  residual 
excitation  energy.  The  Auger  effect  and  creation  of  electrons  and  low  energy  photons  is 
again  seen.  In  the  last  case,  two  photons  can  be  emitted  if  the  residual  excitation  of  the 
second  case  is  greater  than  1  keV.  As  before,  the  residual  excitation  is  dissipated  by 
ftirther  Auger  events  and  electron  production.  In  all  cases,  die  photons  are  assumed  to  be 
emitted  isotropically.  A  photoelectric  event  is  terminal  for  elements  Z  <  12  because  the 
possible  fluorescence  energy  is  below  1  keV.  A  single  event  is  possible  for 
12  <  Z  <31  with  double  fluorescence  possible  for  Z  >  31. 

4.4  Inputs/Outputs 

Input  to  MCNP  consists  of  a  file  containing  the  entire  problem  description: 
geometry,  source,  tallies,  etc.  Appendix  A  contains  a  sample  input  file  for  this  analysis. 

The  specific  entries  will  be  discussed  in  Chapter  S  widi  the  description  of  the  modeled  CT 
source  and  configuration.  MCNP  produces  an  output  data  file  containing  an  entire 
problem  summary.  It  includes  a  copy  of  the  input  file  (fixed  data),  explicit  source 
definitions,  sampling  values  (expected  vs  sampled),  and  a  diorough  analysis  of  each  tally 
including  the  results  of  ten  statistical  checks  for  each  tally  history.  Appendix  B  contains  die 
ouq>ut  file  that  resulted  fi-om  running  the  input  file  in  Appendix  A.  This  example  does  not 
contain  the  fixed  data  since  diat  is  contained  in  the  input  file.  It  also  eliminates  the 
sampling  results  since,  with  such  a  complicated  source,  the  amount  of  information  is 
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overwhelming  and  with  2,000,000  particles  per  source  point  run,  quite  luiiform  and  well 
sampled.  Figure  4.1  shows  the  virtually  indistinguishable  normalized  input  CT  spectrum 
and  die  sampled  spectrum  (for  onfy  66,000  particles)  which  supports  the  decision  not  to 
present  sampling  results. 

In  addition  to  the  regular  output  file,  an  optional  MCTAL  file  can  be  requested.  It 
contains  strict^  formated  tally  and  special  data.  The  file  can  be  can  be  reformatted  with  a 
relatively  simple  FORTRAN  program.  The  reformatted  data  is  better  suited  to  data 
processing.  That  option  was  utilized  in  tiiis  work.  MCNP  produces  the  RUNTPE  file  that 
is  used  to  continue  runs  that  were  interrupted  or  imejqiectedfy  killed.  Frequent  and 
periodic  dumps  to  the  RUNTPE  file  ensure  minimal  loss  of  data  in  case  a  continuation  is 
necessaiy.  The  number  of  saved  dumps  can  be  controlled,  dius  limiting  the  size  of  die 
RUNTPE  file. 
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Figure  4.  1  Normalized  sampled  vs  expected  source  spectra. 
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5.  Computerized  Tomography 


5.1  Introduction 

Tomography  refers  to  the  cross-sectional  imaging  of  an  object  from  either 
transmission  or  reflection  data  collected  by  illuminating  the  object  from  many  different 
directions  (Kak,  1987).  Although  tomography  is  primarily  an  imaging  technique,  that 
particular  aspect  win  not  e}q)]icidy  be  covered  since  diis  diesis  is  concerned  with  the  source 
and  system  configuration  and  their  impact  on  the  non-occupational  dose  outside  the  suite. 
Widi  diat  focus,  die  foUowing  sections  wiH  cover  a  description  of  the  source,  the  detector, 
the  suite  configuration,  and  how  all  of  these  were  modeled  within  MCNP. 

5.2  Source 

A  CT  source  is  unique^  defined  by  its  method  of  scanning,  its  spectrum  and  its 
fan-shaped  beam.  Detector  choice  is  determined  by  source  since,  with  CT,  unlike  other 
radiographic  imaging  processes,  the  characteristics  of  the  source  force  a  definition  of  the 
detector.  (Kak,  191987;  Mackovski,  1983). 

5.2.1  Scan  Modalities. 

In  1973  a  revolutionary  concept  in  tomography,  known  as  computerized  axial 
tomography,  was  introduced  that  provided  an  isolated  image  of  a  plane  within  a  volume 
completely  eliminating  all  other  planes.  Figure  5. 1  (a)  shows  die  operation  of  the  first 
generation  of  scaimers.  They  consisted  of  a  single  source  and  single  detector  that  were 
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(d)  4th  Generation 


Figure  5. 1  Four  scanner  modalities. 

translated  simultaneously  along  strai^t  parallel  paths  (1)  widi  the  taiget  between  them  thus 
producing  a  single  projection.  The  source>detector  set  were  then  rotated  about  the  object 
(2)  and  the  translation  repeated.  The  problem  with  this  method  was  the  relatively  long  total 
scan  periods,  on  the  order  of  several  minutes.  This  was  acceptable  for  relatively  stationary 
targets,  but  insufticient  for  restless  patients  (like  children)  or  abdomenal  images,  where  it 
was  difficult  for  patients  to  hold  their  breath  for  the  necessary  lengtii  of  time. 

The  second  generation  of  scanners  consisted  of  a  thin  fan-shaped  source  beam  witii 
multiple  detectors  tiiat  were  again  moved  simultaneously  along  parallel  axes.  Figure  5.1  (b). 
The  multiple  projections  (to  each  of  the  multiple  detectors)  for  each  lateral  scan  reduced 
the  number  of  rotations  required,  and  thus  the  scan  times,  to  a  fraction  of  a  minute. 
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The  scanner  in  WPMC’s  CT  suite  is  a  third  generation  scanner  with  a  single  fan- 
beam  source  and  multiple  detectors  simultaneous^  rotating  around  the  target.  Figure  5.1 
(c).  This  geometry  eliminates  the  need  for  die  translation  motion.  Scan  times  for  third 
generation  scanners  are  as  low  as  one  second  widi  more  than  1000  projections  taken  in  a 
given  360  d^ee  rotation.  Five  to  seven  hundred  elements  are  contained  in  die  detector 
array,  allowing  a  very  high  resolution  image.  Some  rotational  smearing  may  be  present 
(from  the  continuous  rotation),  but  is  minimal  in  die  final  image  (Kak,  1987). 

The  newest  scanners,  fourth  generation,  have  a  sin^e  fan-beam  source  and  a 
stationary  detector  ring  conqiletely  surrounding  the  target.  The  source  is  rotated  around 
the  target,  within  the  detector  rir^  Figure  5.1  (d). 

The  CT  unit  at  WPMC  is  a  GE  Hi-Speed  Advantage.  Its  fan  beam  has  an  angular 
extent  adjustable  from  25  to  48  degrees  (25  to  48  cm  at  isocenter).  A  35  degree  beam  was 
used  in  die  analysis  as  a  typical  choice  (Poteat,  1995).  The  beam’s  thickness  is  adjustable 
from  one  to  ten  mm.  A  ten  mm  beam  thickness  was  used  for  the  analysis.  The  source 
rotates  63  cm  from  isocenter.  Scan  times  for  WPMC’s  GE  CT  are  adjustable  from  one  to 
four  seconds,  with  the  time  determined  by  the  type  of  slice  and  operating  current.  During 
March,  April  and  M<ty  of  1995  all  scans  were  one  second  long  with  die  exception  of  head 
scans  which  are  done  for  two  seconds  at  half  the  current.  (GE,  1991). 

5.2.2  Detectors. 

The  scan  modalities  of  later  generation  CPs  presented  unique  detector 
requirements.  The  fast  scan  times  required  rapid-response  detector  materials.  Since  the  X 
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rays  impinge  on  the  detector  perpendicular^,  relative^  deep,  and  thus  highly  collimating, 
detectors  can  be  used.  Thus,  this  allows  the  use  of  a  gaseous  ionization  detector  material 
like  xenon.  Solid  scintillation  detector  materials  are  also  common.  The  crystal  traps  most 
of  the  X-ray  photons  (absorption  is  dependent  on  photon  energy  and  crystal  size).  The 
photons  most  likely  undergo  photoelectric  absorption  which  results  in  the  production  of 
secondary  electrons.  Within  the  crystal,  the  kinetic  energy  of  the  secondary  electrons 
transforms  into  flashes  of  li^t  which  are  detected  by  a  photomultiplier  tube  or  solid  state 
photodiode.  (Kak,  1987) 

The  CT  at  WPMC  uses  a  solid  scintillation  material,  possibly  cadmium  tungstate 
(CdWO^.  The  exact  material  composition,  although  confumed  as  solid  state,  is 
proprietary  to  GE.  The  detector  was  modeled  as  CdW04  wWi  a  density  of  7  g'cm^.  The 
material  was  required  to  absorb  99%  of  the  photons  impinging  on  it  to  match  the  machine 
specifications.  As  defined  by  third  generation  scanning,  the  detector  is  an  arc  with  a  radius 
defined  by  tire  source-to-detector  distance,  110  cm.  The  detector  array  of  852  elements  is 
one  meter  long.  (GE,  1991). 

5.2.3  Spectrum. 

Fewell  and  Shuping  published  CT  spectra  for  a  variety  of  manufacturer’s  X-ray 
tubes  (Fewell,  1981).  Their  spectrum  for  the  GE  tube  was  used  for  preliminary  runs,  but 
the  tube  spectrum  for  the  GE  HiSpeed  Advantage  was  obtained  fi'om  GE  and  used  for  the 
final  runs.  Figure  5.2  shows  normalized  spectra  firom  Fewell  and  Shuping  and  fi'om  GE. 
The  GE  tube  consists  of  a  source  of  electrons  impinging  on  a  tungsten  target  which 
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Figure  5.  2  Normalized  CT  spectra  from  Fewell  and  Shaping  (Fewell,  1981)  and  GE 
(Fox,  1995a). 

produces  a  Bremsstrahlung  radiation  spectrum  from  the  collisional  interactions  between  die 
electrons  and  the  tungsten.  The  photon  beam  is  collimated,  narrowing  it  to  one  to  ten  mm 
thickness  as  describedabove.  This  beam  is  passed  throu^  a  Teflon  bowtie-shaped  filter: 
flat  on  one  side  widi  an  IS"*  order  polynomial  surface  on  die  other  side  (see  Figure  5.3). 
The  puipose  of  the  bowtie  filter  is  to  reduce  die  total  energy  of  die  fan  beam  at  the  wider 
angles.  This  is  desirable  since  the  thickness  of  the  human  body  is  less  at  the  hi^er  fan 
beam  angles.  The  filter  reduces  unnecessaiy  radiation  in  the  patient  (Fox,  199Sa).  The 
exact  filter  shape  and  eneigy  dependent  attenuation  coefitients  (propriety:  available  for  the 
anafysis  but  not  release)  were  used  to  modify  die  35°  fan  beam  at  one  degree  intervals. 
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Figure  5.  3  GE's  CT  X-ray  tube  bowtie  filter  shown  with  a  representative  fan  beam. 

5.3  Configuration 

CT  units  may  be  placed  parallel  to  the  walls  or  rotated.  The  rotated  configuration 
is  often  desirable  since  it  gives  the  patient  easier  access  to  die  table  and  also  allows  the 
technician  to  observe  the  patient  during  the  scan.  The  rotated  configuration  is  found  in 
approximately  70  -  80%  of  CT  installations  (Livingston,  1995).  The  predominance  of  the 
rotated  configuration  dictated  its  use  for  the  majority  of  the  runs  (see  Figure  5.4).  As 
described  in  the  introduction,  runs  were  accomplished  with  die  CT  unit  in  the  X-ray  room 
described  by  Metzger.  Those  results  can  be  generalized  to  the  non-rotated  conjuration. 
The  Hi-Speed  Advantage’s  product  data  booklet  supplied  by  G£  contained  a 
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Figure  5.  4  Room,  patient  and  source  configuration  for  sample  input  file.  This  is  the 
configuration  that  was  used  for  the  analyses. 
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recommended  minimum  suite  size  which  was  used  for  the  general  case.  The  exact 
dimensions  and  shielding  in  die  CT  suite  at  WPMC  were  used  for  WPMC  analysis. 

The  patient  lies  on  a  translatable  table  that  is  positioned  within  the  source-detector 
ling.  The  table  was  assumed  to  be  effectively  transparent  to  X  rays  since  it  supports  the 
entire  patient  during  scans  and  a  non-tansparent  table  would  interfere  with  imaging. 

5.4  CT  Model  for  MCNP 

Appendix  A  contains  a  sample  MCNP  input  file.  It  is  for  one  source  point  (beam 
position)  in  the  analysis  of  a  CT  head  slice  in  the  WPMC’s  CT  suite.  It  contains  a  stylized 
human  phantom  consisting  of  a  homogenous  flesh  material  (Snyder,  1978)  in  a  supine 
position  rotated  about  the  vertical  axis.  The  walls  are  sheets  of  gypsum  diywall  on  either 
side  of  a  1/8”  layer  of  lead.  The  ceiling  and  floor  are  each  15  cm  of  concrete.  The  CT 
detector  and  aluminum  structural  supports  are  modeled.  The  typical  CT  fan  beam 
(modeled  as  35  lateral  pencil  beams)  emanating  fi-om  a  single  point  source  is  used. 

5.4.1  Cell  Geometry. 

Section  I  of  Appendix  A  shows  a  list  of  cell  cards:  a  line  in  the  input  file  defining 
some  specific  aspect  of  the  problem.  Cells  are  volumes  of  interest  containing  a  single 
material.  are  built  fi'om  surfaces  (defined  later).  The  first  number  on  a  card  is  its 
unique  numerical  identifier.  The  next  two  numbers  refer  to  tiie  material  contained  in  the 
cell  and  its  density.  A  cell  can  be  defined  as  a  ‘Void”  meaning  all  particles  that  enter  it  are 
killed  (see  cell  20).  The  next  list  of  positive  and  negative  numbers  on  the  card  are  the 
surfaces  defining  the  cell.  The  sign  of  the  surface  indicates  its  sense^  i.e.  up  fi'om,  down 
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from,  to  the  left  of,  to  the  right  of,  etc.  Associating  a  sign  with  a  surface  actually,  tiien, 
defines  the  volume  in  some  direction  relative  to  the  surface.  A  space  between  two  surface 
numbers  indicates  tiie  intersection  of  the  two  volumes.  A  spherical  shell  could  be 
described  as  the  intersection  of  the  region  inside  a  larger  sphere  with  tiie  region  outside  a 
smaller,  concentric  sphere.  A  colon  between  two  numbers  indicates  tiie  union  of  two  (or 
more)  volumes.  A  cell  entry  shown  with  a  pound  sign  (#4),  where  die  number  is  another 
cell  number,  indicates  that  die  cell  contains  everything  but  that  cell.  In  the  previous 
spherical  shell  example,  the  cell  could  have  been  described  as  everything  inside  the  larger 
sphere  but  not  in  the  smaller  sphere  (defined  elsev^o'e  as  a  cell).  The  last  number  is  the 
cell’s  in^ortance  used  for  the  variance  reduction  techniques. 

Note  cell  105,  the  torso,  in  die  appendix.  It  is  the  intersection  of  the  interior  of  an 
elliptical  cylinder  parallel  to  the  y  axis  (suifrice  10)  bound  top  and  bottom  by  two  x-z 
planes  (surfaces  6  and  8).  The  torso  cell  does  not  include  the  lungs,  hence,  #101  and 
#103.  It  is  composed  of  the  homogenous  flesh  material,  m2,  with  a  density  of  0.9869 
g/cm^.  It  has  a  photon  importance  of  1 . 

5.4.2  Surface  Geometry 

Section  n  of  Appendix  A  shows  a  list  of  surface  cards.  The  first  number  in  a 
surface  card  is  again  the  unique  numerical  identifier.  This  is  followed  by  a  letter(s) 
describing  die  type  of  surface  and  a  list  of  numbers  that  satisfy  the  description  of  the 
surface.  For  example,  planar  surfaces  are  indicated  by  a  p,  and  are  described  by  their 
coefficients  in  die  equation  of  a  plane:  Ax  +  By  +  Cz  +  D  =  0  (see  surface  70  or  71). 
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Volumetric  surfaces  (sphere,  ellipsoid,  cone)  are  described  by  their  size  and  the  location  of 
their  center.  Siuface  2  is  an  ellipsoid  whose  upper  half  delBnes  the  top  of  the  head.  The 
coeflScients  (A,  B,  C,  D,  E,  F,  G,  x’,  y’,  z’)  satisfy  the  equation:  A(  x  -  x’  )*  +  B(  y  -  y’  f 
+  C( z- z’ )^  +  2D(x-x’ )  +  2E(y-y’ )  +  2F( z- z’ )  +  G  =  0. 

5.4.3  Materials. 

Section  HI  of  Appendix  A  contains  descriptions  of  the  materials  used  in  the 
problem.  Each  material  is  uniquely  numbered  and  described  as  a  list  of  the  elements  and 
their  atomic  or  wei^t  fraction.  (For  other  than  photon  anafyses,  the  desired  cross  section 
library  would  also  be  specified.)  Material  11  (mil),  cadmium  tungstate,  an  example  of  a 
material  specified  by  atomic  fiiactions,  has  the  chemical  formula  CdW04.  Material  1  (ml), 
lung  tissue,  an  exanqrle  of  a  material  specified  by  weight  fractions,  is  composed  of  10.21% 
hydrogen,  10.01%  carbon,  2.8%nitrogeri,  etc. 

5.4.4  Source. 

MCNP’s  general  source  utility  was  used  to  describe  the  CT  source.  It  requires  the 
energy  spectrum  and  the  beam  positions,  directions  and  shape.  Since  MCNP  does  not 
have  the  capability  to  directty  model  the  continuous,  rotating,  fan-shaped  beam,  the  CT 
source  was  modeled  as  32  discrete  point  sources,  each  consisting  of  35  individual  pencil 
beams  which  simulated  the  fan  sh^d  beam  (see  Figure  5. 5).  Each  source  point  was  run 
individually  because  of  the  detector  configuration.  Recall  that  the  detector  rotates  opposite 
to,  but  syncronous  with  the  source.  MCNP  could  not  model  a  detector  that  moved  as  the 
source  was  sampled  from  different  beam  positions.  The  source  is  on  a  radius  larger  than 
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Figure  5.  5  View  of  the  modeled  CT  source  in  the  x-z  plane.  Note  the  labels  at  the 
source  positions. 


the  detector  radius.  Attempting  to  place  a  continuous  detector  in  the  configuration  would 
block  the  source,  thus  the  32  individual  runs  necessary  to  model  a  single  source  rotation. 
Section  IV  of  Appendix  A  contains  the  complete  modeled  source. 

The  first  uncommented  card  is  the  general  description  of  the  source.  Each  variable 
(eig,  pos,  etc.)  points  to  an  associated  distribution,  dxxx.  The  position  (pos  =  dl)  is 
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described  as  a  distribution  of  positions,  each  with  equal  probability  (spl).  Since  the 
appendix  is  an  input  file  for  only  one  of  the  discrete  beam  positions,  die  position,  that  is, 
the  point  of  origin  for  the  source  point,  for  each  of  the  35  pencil  beams  is  die  same.  The 
beam  shape  (dir  =  d2),  for  each  of  the  35  individual  beams  is  a  veiy  narrow  cone,  or 
pencil,  beam.  The  si2  card  is  the  cosine  of  the  angular  extent  of  the  beam,  one  degree.  It 
represents  the  distribution  values  of  die  particles  for  any  selected  pencil  beam.  Its 
probability,  on  the  sp2  card,  is  always  one. 

The  beam  direction  (vec  =  Qios  =  d3)  is  defined  by  a  vector  passing  through  die 
source  point  with  a  direction  determined  by  its  position  within  the  fan  beam.  Envision  a 
line  fi'om  the  source  point  to  isocenter.  The  fan  extends  17.5  degrees  to  either  side  of  this 
line.  The  beam  vector  is  determined  by  the  beam's  position  in  the  fan  beam  but,  defined 
relative  to  the  problem  geometry  axes.  Each  of  the  pencil  beams  (initially  produced  with 
the  pos  variable)  is  associated  with  one  vector  in  the  list  of  distribution  labels,  ds3.  The 
labels  point  to  the  individual  vector  descriptors:  si201,  si202,  etc.  All  of  the  vectors,  again, 
have  equal  probability  (sp201  =  sp202  =  . . .  =  sp235  =  1). 

The  eneigy,  or  distribution  of  energy,  is  input  with  the  erg  variable,  erg  =  Q)os  = 
d4.  As  described  above,  the  primary  source  passes  through  die  bowtie  filter  (uniform 
about  its  center)  to  reduce  its  eneigy  at  the  outer  fan  beam  angles.  The  spectral  differences 
produced  in  each  pencil  beam  by  the  filter  are  seen  in  the  18  different  spectra  pointed  to 
with  ds4,  the  distribution  labels,  which  then  point  to  the  individual  sets  of  spectal  values 
(sill7,  sill6,  etc.). 
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The  last  hvo  cards  in  the  source  description,  dfO  and  deO,  are  the  flux-to-e?^sure 
conversion  factors.  These  cards  are  used  to  scale  the  talfy  as  a  function  of  particle  energy. 
They  are  typically  used  to  convert  the  tally  to  an  eneigy-dependent  dose  equivalent  rate  in 
rem  per  hour.  The  conversion  factors  shown  here  (used  by  Metzger  in  his  diagnostic  X- 
ray  system  analysis)  are  the  ICRP-2  standard  (MCNP,  1993)  provided  with  MCNP,  scaled 
by  rem-to-Roentgen  factors  and  converted  from  hours  to  seconds.  This  results  in  a  tally 
reported  as  an  exposure  rate  in  Roentgen  per  sec  which  was  converted  to  gray  per  sec  (see 
Section  2.8.6). 

5.4.5  Tallies. 

Section  V  contains  the  tallies  (see  section  4.2.2. 1)  that  were  computed  for  this 
example.  Flux  tallies  are  identified  by  a  numerical  designator  ending  in  “5”.  The  “p” 
indicates  a  photon  (X-ray)  tally.  In  fius  example,  with  the  phantom  rotated  about  the 
vertical  axis,  tallies  were  computed  along  hvo  different  axes  sets,  the  first  set  defined  by  the 
bo(fy  axis  and  source  plane  and  the  second  set  by  the  room  walls  (see  Figures  5.4  and  5.5) 
30  cm  outside  the  walls.  In  addition,  floor  and  ceiling  tallies  were  done,  again,  30  cm 
outside  the  enclosure.  Not  shown  in  the  example  file  is  a  ring  detector  that  was  placed  30 
cm  within  the  source  ring.  This  talfy  was  used  a  baseline  for  kerma  rate  computation 
(section  6.1.2). 

The  “dd”  card  describes  the  Russian  roulette  game:  the  first  number  is  the  fraction 
of  the  average  weight  below  which  particles  will  be  rouletted.  Example:  a  value  of  .01 
would  kill  an  particles  with  weights  below  1%  of  the  average.  The  “0”  entry  here  indicates 
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that  low  weight  particles  will  not  be  killed.  The  “le25”  entiy  is  used  only  to  discriminate 
which  score  histories  to  write  to  the  ou^ut  file.  A  very  large  value,  as  seen  here,  will 
prevent  all  entries  fi-om  being  written  to  the  output  file.  This  was  done  solely  to  minimize 
ou^ut  file  size.  Since  obtaining  the  spectrum  of  the  escaped  dose  was  a  goal  in  this 
analysis,  the  “eO”  card  was  used  which  lists  the  energy  bins  desired.  The  “prdmp”  card 
indicates  the  fi'equency  file  continuation  file  and  talfy  dumps  are  written.  As  with  the  “dd” 
card,  the  high  value  limits  the  number  of  dumps,  again,  to  minimize  file  size.  This  card 
also  specifies  that  a  MCTAL  file  (see  section  4.4)  be  written.  The  last  card,  “nps”, 
specifies  the  number  of  particles  to  run. 

5.5  MCNP  Output 

Appendix  B  contains  the  output  file  produced  by  a  run  using  the  input  file  in 
Appendix  A.  As  stated  previously,  it  does  not  contain  the  fixed  data  or  the  sampled  values. 
The  fixed  data  and  sampled  values  would  be  located  before  the  output  seen  in  the 
appendix. 


5.5.1  Photon  Creation. 

Section  I  of  Appendix  B  is  a  summary  of  photon  creation  information.  Note  that 
many  of  the  listed  values  are  those  produced  by  the  different  variance-reduction  techniques 
(weight  cutoff,  DXTRAN,  etc.).  This  section  also  contains  run-time  figures:  computer 
time,  source  particles  per  minute  and  random  number  generation  figures.  This  particular 
run  used  almost  300  million  random  numbers  for  the  2  million  particles. 
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5.5.2  Cell  Activity. 


The  information  contained  in  Section  n  is  useiiil  for  improving  cell  importances. 
Idealfy,  the  number  of  particle  tracks  entering  a  cell  (tracks  coliunn)  would  be  equal  for  all 
cells  to  equalize  error  contributions  per  cell.  Large  differences  could  indicate  a  need  for 
further  geometric  splitting  of  either  of  two  adjacent  cells.  The  above  guidelines,  however, 
must  be  tempered  with  the  knowledge  of  the  problem  geometry.  For  example,  cell  9 
shows  a  much  larger  number  of  entering  tracks  dian  cell  8.  Cell  9,  however,  is  the  cell 
containing  all  of  the  space  within  the  room  diat  is  not  explicit^  another  object  (torso, 
detector,  etc.).  It  would  naturally  have  a  greater  number  of  tracks.  The  tracks  column 
increments  every  time  a  particle  enters  a  cell,  even  if  it  is  re-entering  the  cell.  The 
population  column  lists  a  track  only  the  first  time  it  enters  a  cell.  Comparison  of  the 
population  column  with  the  tracks  column  can  provide  a  rough  estimate  of  die  amount  of 
backscatter  into  a  cell. 

5.5.3  Tally 

All  of  the  information  on  a  given  talfy  is  grouped  together.  Section  HI  shows  one 
example  of  a  tally  summary.  It  starts  with  a  description  of  die  tally  -  name,  location,  and 
type.  If  the  tally  has  been  modified  by  a  dose  fimcdon  (deO  and  dfO,  section  S.4.4)  a 
comment  to  indicate  that  will  appear.  Recall  that  the  energy  spectrum  of  the  tally  was 
desired.  This  information  is  listed  (with  its  error)  for  the  collided  and  uncoUided  photon 
flux.  The  total  tally  and  error  are  provided  at  die  bottom  of  the  spectrum.  The  score 
contributions  by  cell  can  ^o  assist  importance  function  determination.  Cells  that 
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contribute  little  to  the  score  do  not  require  large  importances.  The  score  misses 
information  summarizes  the  zero  score  tallies  and  their  inqMrtance  to  the  tally  value.  The 
next  table  shows  die  impact  the  largest  scoring  particle  histoiy  in  the  run  would  have  on 
tally  and  error  if  it  were  to  occur  in  the  next  history. 

MCMP  performs  ten  statistical  checks  of  each  talfy  to  help  the  user  test  for 
reliability.  The  desired  and  observed  behavior  of  the  tally,  error,  variance  of  the  variance, 
figure  of  merit  and  probability  density  fimction  tail  are  compared  and  given  a  pass/fail 
score. 

The  last  information  provided  for  the  individual  tallies  is  the  tally  density  chart.  It  is 
a  graphic  representation  of  the  unnormed  probability  density.  It  should  be  relatively 
smooth  with  no  holes  (which  could  indicate  undersampled  regions  of  the  probability 
density  fimction).  The  probability  density  chart  is  not  provided  if  the  tally  passes  all  ten 
statistical  checks. 

Eveiything  described  above  is  repeated  for  each  tally.  The  last  entiy  in  Section  in  is 
a  summary  of  the  statistical  checks  of  all  tallies. 

5.5.4  Tally  Fluctuation  Charts. 

The  tally  fluctuation  chart  in  Section  IV  is  a  cumulative  history  of  the  tally  at 
regular  intervals  in  the  run.  The  mean,  error,  variance  of  the  variance  and  %ure  of  merit 
should  be  examined  for  convergence.  Erratic  results  could  indicate  a  poorly  defined 
problem,  the  need  for  additional  particle  histories  or  inappropriate  variance  reduction 
parameters. 


5-16 


5.5.5  Final. 


This  last  section  contains  run  tennination  information.  This  run  ended  when  two 
million  particle  histories  had  been  completed.  It  could  also  show  that  a  run  had  been 
terminated  because  of  geomehy  problems,  missing  inputs,  etc.  Section  V  again 
summarizes  computer  run-time  information.  This  last  section  is  particular^  useful  during  a 
run.  The  UNIX  ‘tail”  command  prints  the  last  lines  of  a  file,  i.e.  run  status,  to  the  screen 
eliminating  the  need  to  open  a  file  and  scroll  to  the  bottom,  preventing  possible  file 
corruption. 
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6.  Calculations  and  Metiiod 


Metier  provided  a  method  to  determine  dose  outside  a  diagnostic  X-ray  facility 
using  Monte  Carlo  anatysis  techniques  (Metier,  1993).  His  method  is  the  basis  for  the 
method  used  in  this  study  to  determine  dose  outside  a  conqmtciized  tomography  (CT) 
suite.  The  first  section  below  describes  the  calculations  necessary  for  both  an  X-ray  room 
and  a  CT  suite.  The  following  sections  describe  how  Metier’s  results  were  reproduced 
for  a  standard  X-ray  room.  The  last  section  describes  the  CT  modeled  widun  die  standard 
X-ray  room  with  the  body  and  source  axes  parallel  to  die  wads.  Within  that  configuration, 
the  importance  of  variations  in  the  phantom  and  CT  detector  were  investigated.  Widi  those 
investigations  completed,  the  set  of  32  beam  positions  (representing  discrete  points  around 
a  full  scan  rotation)  for  the  CT  in  the  X-ray  room  was  run. 

6.1  Calculations 

The  following  sections  present  the  calculations  and  data  manipulation  done  to 
obtain  the  dose  outside  the  CT  suite.  In  general,  the  dose  calculation  is  determined  fi'om 
the  MCNP  tallies,  X-ray  tube  spatial  beam  width  and  potential  usage  factors,  CT  tube 
current  protocols,  and  average  air  kerma.  The  specific  calculations  are  described  below. 

6.1.1  Work  Load. 

Bodi  X-ray  and  CT  units  can  be  operated  in  different  tube  potential,  current  and 
beam  size  modes  determined  by  the  image  required  (Nfichael,  1995;  GE,  1995a).  Work 
load  refers  to  the  degree  of  use  of  a  particular  source  mode  or  the  proportion  of  time  it  is 
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operated  in  each  mode  (NCRP,  1976).  It  is  important  to  include  die  workload  in  the 
anafyses  because  it  direct^  impacts  die  total  dose.  The  area  to  be  imaged  determines  the 
beam  width,  operating  tube  potential,  number  of  rotations  per  scan  and  beam  current.  The 
total  woridoad  for  a  given  time  int^al  is,  therefore,  a  summation  over  the  beam  width  and 
tube  potential  usage  factors  (Ubw  and  Utp)  and  the  the  area  to  be  imaged: 


Work  Load 


mA-min 
time  interval. 


=  Yu  Yu  yf  ScanSets  'j 

Z^[a„eteMvalJ, 

1  J  K  ^ 

f  Rotations^  ^  ^  ^  ^ 

X  -  X  Current.  (mA) 

I  k''  ' 


ScanSet  Jy. 


X  RotationTimej^(sec)  x 


l(niin) 
60  (sec) 


(6.1) 


Unique  woildoads  were  determined  for  the  head  and  abdomen  scans. 


6. 1. 1. 1  Tube  Potenial  Usage  Factor. 

Typically,  X-ray  tube  potential  is  be  varied  from  45  kVp  to  140  kVp, 
determined  by  the  size  or  area  of  the  body  to  be  imaged.  Tube  potential  usage  factors  are 
determined  by  the  length  of  time  each  potential  is  used.  The  times  are  normalized  to  one 
and  jqiplied  to  the  workload  and  kerma  (air  kerma  varies  witii  tube  potential;  NCRP,  1989). 
The  tube  potential  usage  factors  are  different  for  tiie  chest  (standing  patient  with  wall- 
mounted  recording  film)  and  radiographic  (siqiine  patient  on  table)  rooms  which  Met2ger 
obtained  from  published  reports.  The  radic^aphic  room  usage  factors  used  by  Metzger 
are  seen  in  Figure  6.1.  Since  tube  potential  usage  factors  were  not  available  for  CT  units, 
those  used  for  titis  analysis  were  based  on  Qrpical  usage  protocols  at  WPMC  (Michael, 
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Figure  6. 1  Normalized  radiographic  room  tube  potential  worldoad  used  by  Metzger. 


1995).  As  can  be  seen  in  Table  6.1,  Iheir  CT  tube  was  operated  exclusive^  at  120  kVp 
corresponding  to  a  usage  factor,  Utp  ,  of  one. 


Table  6.  1  GE  HiSpeed  Advantage  CT  Protocols 


Scan  Type 

1  Potential 

1  (kVp) 

Current 

(mA) 

# 

Rotations 

Head 

■BB 

140 

2 

Sinus 

280 

1 

Body 

- 

- 

- 

- 

Chest 

120 

280 

1 

35 

Abdom  en 

120 

280 

1 

40 

Pelvis 

120 

280 

1 

20 

Chest/Abd/Pel 

120 

280 

1 

80 

Abd/Pel 

120 

280 

1 

50 

Liver 

120 

280 

1 

40 

Renal 

120 

280 

1 

30 

Pan 

120 

280 

1 

30 

6. 1.1.2  Current  Workload. 


X-r<or  tube  current  is  again  determined  by  the  patient  area  to  be  imaged. 
NCRP  49  contains  average  current  work  loads  for  diagnostic  X-ray  rooms  which  Met2ger 
modified  fi*om  the  suggested  1000  mA-min  per  week  down  to  400-500  mA-min  per  week 
based  on  the  use  of  modem,  faster  expostue  films.  These  figures,  however  cannot  be 
applied  to  CT  units.  The  CT  usage  summaiy  fi-om  WPMC,  averaged  over  March,  April 
and  May  1994  and  die  CT  protocols  (Tables  6.2  and  6.1,  respective^),  was  used  for  this 
analysis.  (The  number  of  individual  bo^  scans  are  a  percentage  of  die  total  number  of 
reported  bo^  scans.)  Current  workload  was  determined  separate^  for  the  head  and 
abdomen  scans  and  qiplied  to  the  head  and  abdomen  MCNP  results.  The  total  current 
woiMoad  for  the  CT  unit,  about  16,000  mA-min  per  week,  is  much  larger  than  the  X-ray 
current  workload  of 400  to  500  mA-min  per  week. 

6. 1.1.3  Beam  Size. 

Metzger  used  a  conic  X-ray  beam  codimated  to  1536  cm^  at  the  patient. 
The  G£  HiSpeed  Advantage’s  fan-shaped  beams  an^e  is  adjustable  fi'om  25  to  48  degrees 
and  beam  width  fi'om  one  to  ten  mm.  The  beam  ang|le  is  reduced  by  blocking  the  source 
which  thus  reduces  die  total  energy  out.  Since  data  was  not  available  to  quantify  the 
reduction,  the  original  source  strength  was  used  with  a  35°  fan  angle  (approximatefy  39  cm 
at  isocenter),  10  mm  wide  as  a  worst-case  assumption.  Since  the  source  utility  in  MCNP  is 
only  able  to  model  point  sources  uniform  about  a  central  axis,  the  narrow  beam  width 
required  modeling  the  fan  beam  as  35  individual  beams. 
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Table  6.  2  Summary  of  WPMC's  CT  Usage 


1  Scanl^ 

Number  of  Scan  Sets 

Total 

week 

WorkLoad  I 
mA-min  wk* 

March 

Aprfl 

hfoy 

Head 

101 

110 

yi9 

2863 

1  Sinus 

■B 

44 

38 

133 

1671 

1  Head  Total 

4534 

Body 

243 

204 

266 

713 

- 

- 

Qiest  (3%) 

7.3 

6.1 

8.0 

21.4 

1.6 

269 

Abd(40%) 

97.2 

81.6 

106.4 

285.2 

21.9 

4095 

Pelvis  (3%) 

7.3 

6.1 

8.0 

21.4 

1.6 

154 

Chest/Abcl/Pel(4%) 

9.7 

8.2 

10.6 

28.5 

2.2 

819 

Abcl/Pel(40%) 

sn.i 

81.6 

106.4 

285.2 

21.9 

5119 

liver 

11 

18 

7 

36 

2.8 

517 

Renal 

15 

7 

12 

34 

2.6 

366 

Pancreas 

5 

6 

1 

12 

0.9 

129 

1  Abdomen  Total 

11468 

6.1.2  AirKerma. 

Average  air  kerma  is  a  function  of  tube  potential  and  distance  from  the  source. 
NCRP  Report  102  provides  a  list  of  average  air  kerma  rates  by  X-ray  tubes  potential  and 
distance  from  the  source.  Given  the  MCNP  tallies  for  a  known  source  for  two  locations,  A 
and  B,  air  kerma  can  be  detennined  at  detector  location  B  if  the  air  kerma  at  detector 
location  A  is  known  simply  by  scaling  the  average  air  kerma  from  location  A  by  a  ratio  of 
the  detector  tallies 

TaIIy„ 

Air  Kerma  o  =  Air  Kerma  .  x  - £-  (6,2) 

Tally^ 


6-5 


This  method  of  relating  keimas  via  a  detector  tally  ratio  eliminates  constants  produced  from 
converting  from  d^ees  to  radians,  etc.  It  is  only  valid  for  exposure  tallies:  flux  and 
fluence  tallies  cannot  be  related  to  kermas  because  die  particles  in  the  flux/fluence  have 
different  eneigies.  This  method  was  applied  to  die  anatysis  by  placing  a  33  cm  ting 
detector  on  the  same  plane  as  and  concentric  within  the  63  cm  CT  source  ring  (see  Figure 
5.5).  The  eiqiosure  tally  at  that  rii^  detector  could  be  associated  with  average  air  kerma  30 
cm  from  a  120  kVp  source  (17.8  centi-Gray  per  100  mA  sec  =  0.1068  Gy  per  mA-min; 
NCRP,  1989),  and  thus  could  produce  a  conversion  from  the  exposure  tallies  to  air  kermas 
at  those  locations. 

Air  kerma  is  related  to  a  material  kerma  by  a  ratio  of  mass  energy-absoiption 
coefficients.  For  example,  consider  the  conversion  from  air  to  muscle  tissue  kerma: 

(l^en  P)tia8ae  ^  j ^  30/^  3^ 

(M-en  P)air 

(valid  from  4  keV  to  10  MeV,  Attix,  1986).  This  is  can  be  used  to  convert  dose  (grey)  to 
equivalent  dose  (sievert).  To  convert  the  exposure  tallies  to  effective  dose  for  the  dim 
badge  comparison,  the  method  described  in  £q.  (6.2)  was  used.  The  resulting  dose  in  air 
was  multiplied  by  the  radiation  and  tissue  westing  factors  (see  section  2.8.7)  used  in 
10CFR20  (tiiose  advocated  by  ICRP)  to  obtain  an  effective  dose  in  skin. 


6-6 


6.1.3  Final  Calculation. 


A  complete  CT  scan  is  modeled  with  32  equally  spaced  beam  positions  around  the 
source  ring.  Taking  die  average  of  the  32  tallies  at  each  detector  gives  the  equivalent  tally 
of  a  single  scan.  Initial  results  indicated  that  some  areas  of  the  scan  required  additional 
beam  positions  and  others  could  be  eliminated  (see  Chapter  7  for  a  list).  These  were 
incoiporated  by  appropriately  weighting  the  new  points  and  the  original  points  surrounding 
them.  For  example,  a  wei^t  of  1/32  is  assigned  to  each  of  the  32  equally  spaced  beam 
positions.  If  an  additional  point  was  placed  between  beam  positions  11.25  and  22.5, 

Figure  6.3  shows  the  adjusted  weights,  assuming  each  beam  position  contains  the  area 
from  halftvay  between  itself  and  its  ne^bor  on  each  side.  The  errors  are  given  the  same 
weights  as  die  exposure  tallies.  The  composite  error  is  thus 
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1/32  1/32  1/32  1/32 

Old  weights 


Figure  6.  2  Adjusted  weights  with  the  addition  of  points. 

where  Si  is  die  total  error  for  a  given  tally,  Sjd  is  the  is  the  error  reported  for  a  single  tally 
per  particle,  Wi  is  the  tally  weight  and  Sxtot  is  th®  composite  weighted  exposure  talfy  error 
from  all  beam  positions  at  a  single  detector  location. 

6.2  X-Rav  Source  in  X-Rav  Room 

6.2.1  Reproduction  of  Metzger’s  Results. 

Metzger’s  Monte  Cario  analysis  of  a  diagnostic  X>ray  room  was  duplicated 
primarily  as  a  training  exercise,  but  also  however,  to  ensure  the  his  method  was  property 
interpreted.  Metzger  performed  the  dose  calculations  for  the  chest  and  die  radiographic 
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rooms  each  with  1/16”  and  1/32”  lead,  but  it  was  sufBcient  here  to  duplicate  only  the  1/16” 
lead  thickness.  These  MCNP  runs  were  accomplished  using  input  files  provided  by 
Metier.  Any  differences  noted  would  be  a  result  onfy  of  the  computer  used  and  its 
random  number  generator.  Tables  6.3.  A  and  6.3.B  present  a  comparison  of  the  results. 


Table  6.  3  Comparison  with  Metzger's  Results 


|a:  Chest  Room,  1/16"  Lead 

Air  Kama,  Gy 

DifT 

Error,  % 

WMiin 

Pos 

Danis 

Metiger 

% 

Datvis 

Metager 

Error?* 

Posx 

1.46E-09 

1.5E-09 

2.7 

5.2 

8 

Y 

Negx 

1.68E-09 

1.5E-09 

12.2 

9.6 

8 

Y 

Point 

Negy 

1.37E-07 

1.4E-07 

2.1 

1.7 

2 

Y 

Del 

Posy 

3.62E-10 

3.6E-10 

0.6 

13.1 

10 

Y 

CeOiiig 

1.76E-09 

1.9E-09 

7.2 

35.9 

18 

Y 

Floor 

3.12E-09 

3.1E-09 

0.8 

15.0 

15 

Y 

Ring 

Negy 

1.45E-07 

1.4E-07 

3.8 

1.6 

2 

Y 

Det 

Posy 

3.49E-10 

3.6E-10 

3.2 

10.1 

10 

Y 

Source 

0.0085 

0.0085 

0.0 

0.1 

0.9 

Y 

B:  Radiographic  Room,  1/16" 

Lead 

(Using  adjiKted  source  tal^  source 

taB3^4) 

Afa*Keniia,Gy 

DifT 

Error,  % 

WitMn 

Pos 

Dands 

Metzger 

% 

Dands 

Metzger 

Error?* 

Posx 

5.98E-11 

5.5E-11 

8.6 

15.7 

9 

Y 

Point 

Negx 

5.17E-11 

5.5E-11 

6.0 

13.4 

9 

Y 

Det 

Ceilmg 

2.07E-10 

1.8E-09 

88.5 

16.7 

11 

N 

Floor 

6.96E-09 

6.8E-09 

2.3 

9.8 

8 

Y 

Source 

4.69E-03 

4.7E-03 

0.0 

0.1 

0.1 

Y 

Ring 

Posy 

3.86E-12 

3.7E-12 

4.4 

12.6 

23 

Y 

Det 

Negy 

l.llE-10 

1.6E-11 

596.2 

91.0 

22 

N'’ 

aEnorCritena;  Can  etrorbss  overly?  ^^ative  and  positive  y  had  different  ling  detector  ladii. 


The  chest  room's  source  detector  (see  section  6.1.2)  was  placed  one  meter  jfrom 
the  source  in  the  center  of  the  beam.  This  presented  a  problem  widi  the  supine  phantom  in 
the  radiographic  room,  however,  since  it  would  placed  the  detector  in  the  middle  of  the 


6-9 


torso.  Metier  placed  it,  instead,  50  cm  away  and  modified  the  results  using  the  inverse 
sqared  law,  i.e.  the  distance  was  halved,  therefore,  to  correct  to  a  one  meter  position,  the 
exposure  is  divided  by  four.  Most  agreed  wifiun  statistical  error.  There  are  some  notable 
exceptions,  however.  The  negative-y  tally  ring  detector  for  the  radiographic  room  was 
discovered  to  have  its  radius  much  smaller  than  the  positive-y  talfy  ring  detector,  producing 
an  talfy  with  high  error  that  differed  significant  fi-om  Metzger’s  result.  The  ceiling  results 
for  the  radiographic  room,  for  an  unknown  reason,  exhibited  hi^  error  and  veiy  different 
results. 

6.2.2  MCNP  Versions. 

The  above  runs  were  completed  with  Version  3a  of  MCNP.  When  Version  4a 
came  available,  the  runs  were  repeated.  The  exposure  and  errors  for  all  detectors  were 
identical  to  all  reported  decimal  places.  Some  of  the  tests  that  follow  were  accomplished 
with  Version  3a  and  some  with  Version  4a.  Given  the  exact  agreement  of  results,  those 
runs  done  with  Version  3a  were  not  repeated  with  Version  4a. 

6.3  CT  Source  in  X-Rav  Room 

The  following  sections  describe  the  tests  that  were  done  to  define  and  improve  the 
inputs.  Unless  otherwise  noted,  they  were  ah  done  using  the  CT  source  and  configuration 
modeled  wifitin  the  X-ray  room  described  by  Metzger. 
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6.3.1  Phantom  Considerations. 


The  phantom  used  for  this  analysis  is  the  one  defined  in  the  Medical  Internal 
Radiation  Dosimetiy  (MIRD)  Pamphlet  No.  5  (Snyder,  1978).  The  pamphlet  includes 
madiematical  descriptions  of  die  entire  hrnnan  body,  including  a  stylized  total  body,  organs, 
and  bones.  The  pamphlet  also  contains  the  elemental  composition  of  flesh,  bone  and  lung 
tissue  used  in  the  anafysis.  With  such  complete  phantom  information,  the  questions  were: 
“How  much  detail  is  necessaiy?  Model  the  entire  bo<fy  or  is  the  head  sufBcient?  How 
much  of  die  skeleton,  if  any,  is  necessaiy?  hiclude  the  oigans?”  To  satisfactorily  answer 
these  questions,  test  runs  were  done. 

The  elemental  composition  and  density  of  internal  oigans  reported  in  die  pamphlet 
are  nearly  the  same  as  the  flesh  tissue,  with  the  exception  of  the  lungs  which  have  a  much 
lower  density.  Therefore,  the  only  internal  oigans  retained  were  die  lungs.  The  MIRD 
head  consists  of  an  elliptical  cylinder  topped  by  half  an  ellipsoid.  The  remainder  of  the 
body  consists  of  another  elliptical  cylinder  for  the  torso  and  two  truncated  cones  for  the 
legs.  Problems  with  MCNP  geometiy  were  encountered  when  the  leg  cones  were  used. 
Th^  were  changed  to  a  single  half  ellipsoid  of  the  same  total  length,  width  and  thickness 
(80  cm  X  40  cm  X  20  cm)  resulting  in  a  volume  2.7%  greater  diat  die  cones. 

6.3.I.I  Head  Only  vs  Whole  Bo^. 

To  determine  if  using  only  the  head  would  produce  significantly  different 
results  fi-om  using  the  entire  bo(fy,  an  MCNP  run  using  (a)  only  the  top  of  the  head  and 
face  and  (b)  the  entire  bo^  were  compared.  Table  6.4  lists  those  results.  The  floor  and 
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positive  and  negative-y  exposures  showed  a  significant  difference  indicating  that  inclusion 
of  the  entire  body  was  necessary,  hi  addition  to  the  above  comparison,  exposure  fi'om 
detectors  positioned  off  of  die  head  (positive-y)  and  off  of  die  feet  (negative-y)  within  the 


Table  6.  4  Head-Only  to  Whole  Body  Comparison 


=========== 

PositiMl 

HeadOki^ 

WhdeBody 

Diff;% 

WUUn 

Erroi? 

Detector  Taify, 
Rs'^  paitkle'^ 

Flrror, 

% 

Detector  Taify, 
Rs'^  parikle'^ 

Error, 

% 

Posx 

3.06&18 

7.6 

3.11E-18 

7.5 

1.6 

Y 

Negx 

4.66&18 

8.4 

4.82E-18 

8.7 

3.4 

Y 

CdHng 

4.59E-20 

15.7 

5.48E-20 

22.7 

19.4 

Y 

Floor 

5.48E-19 

9.8 

4.25E-19 

7.2 

22.4 

N 

Source 

l.OOB-12 

8.2 

9.98E-13 

8.3 

0.2 

Y 

Posy 

3.40E-20 

10.1 

5.96E-20 

38.0 

75.3 

N 

Negy 

1.46E-20 

18.0 

1.04E-20 

27.78 

28.8 

Y 

room  were  compared.  Both  exhibited  the  same  behavior  as  the  source  was  rotated,  i.e. 
relatively  luiiform  about  a  mean.  The  ratio  of  positive-y  exposure  to  negative-y  exposure, 
however,  changed  firom  about  two,  head-onfy,  to  about  six,  entire  body.  This  reinforced 
the  decision  to  include  the  entire  bo^,  particularly  for  the  analysis  inside  the  CT  suite. 

6.3. 1.2  Skeleton. 

The  importance  of  including  the  skeleton  in  the  phantom  was  anafyzed 
using  onfy  the  head  phantom  since  it  had  the  hugest  jn'oportion  of  bone  to  tissue  (by 
volume)  and  thus  would  be  most  affected  by  neglecting  the  bone.  The  skull  was  defined 
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Top  View 


Side  View 


Figure  6. 3  Top  and  side  view  of  head  phantom  showing  location  and  size  of  skull. 

as  two  slighlfy  non-concentiic  ellipsoids  (Figure  6.3).  The  homogeneous  flesh  tissue  was 
used  inside  and  outside  the  skull.  Ring  detectors  were  placed  outside  the  walls,  ceiling  and 
floor  and  inside  the  room,  20  cm  on  either  side  of  isocenter  on  die  body  axis  (gantiy 
detectors) 

Table  6.5  shows  the  differences  in  exposure.  Note,  in  particular,  the  positive  and 
negative-y  results;  their  difference  is  much  smaller  than  the  error.  Appendix  C,  Figures 
C.l  duough  C.8  show  the  normalized  spectra  (without  the  roentgen  conversion)  obtained 
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Table  6.  5  Comparison  of  Bone-In  to  No-Bone  Detector  Tallies 


Poaikn 

Bait 

NoBdie 

% 

vmi 

Enot?* 

DetedorTd^ 

Rs^potide^ 

Emr, 

% 

DetednrTdfyy 
Rs'^  pariide^ 

lira; 

% 

Bdbx 

&171E-18 

10.3 

1.0285-17 

11.1 

20.54 

N 

Negx 

6i09(E-18 

121 

739(E-18 

129 

17.58 

Y 

52155-18 

10.3 

537^18 

81 

295 

Y 

Eoor 

3.18(E-17 

7.0 

3.5935-17 

7.9 

11.58 

Y 

Posy 

5.146&20 

25.3 

4.83®-20 

186 

641 

Y 

Negy 

21095-20 

B.9 

2140B-20 

17.7 

1.43 

Y 

Neg^nliy 

a7Q55-16 

0.8 

1.3885-15 

0.7 

3728 

N 

l23UErl5 

0.5 

1.8225-15 

0.6 

3250 

N 

*  Wtlinentrcaiteiia:  I>>1faeenxTbarsovedap? 


for  each  of  the  detectors  comparing  the  results  for  bone-in  and  no-bone.  The  error  bars  in 
Figure  C.  1  shows  that  the  difference  falls  within  error.  The  only  spectra  that  show 
noticeable  differences  are  the  gantry  detector  tallies,  Figures  C.7  and  C.8,  but  they  are 
inside  the  room  and  the  interest  in  this  thesis  is  the  dose  outside  the  walls.  Since  the 
difference  in  the  area  of  interest  (inside  the  room)  is  insignificant  the  skeleton  was  not 
included  in  the  phantom. 

6.3.2  CT  Detector. 

Proper  modeling  of  the  CT  detector  was  critical.  Any  direct  X  rays  penetrating 
through  or  not  impinging  on  the  phantom  had  to  hit  the  detector.  According  to 
manufacturers  specifications,  the  detector  is  a  99%  absorber  (G£,  1991).  As  described 
earlier,  third  and  fourth  generation  scanners  use  either  ionizing  gas  (xenon)  or  solid 
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scintillation  detectors.  Dr.  Stanley  Fox,  CT  Marketing,  General  Electric,  verified  that  the 
Hi-Speed  Advantage  CT  scanner  at  WPMC  uses  a  solid  state  detector,  but  he  could  not 
reveal  its  exact  composition  since  tiie  information  is  proprietaiy  (Fox,  1995).  Susan 
Poteat,  in  WPMC’s  radiation  safety  office,  confirmed  that  cadmium  tungstate  (CdW04) 
was  a  possible  candidate.  The  CRC  Handbook  of  Material  Science  (Lynch,  1975), 
however,  did  not  list  a  density  for  this  material.  Based  on  the  densities  of  carbon  tungstate 
and  several  cadmium  compounds,  an  approximate  density  of  7  g/cm^  was  tested. 
Measurement  of  the  detector  array  unit  in  WPMC’s  CT  indicated  the  array  was  no  thicker 
than  five  cm.  An  MCNP  run  was  prepared,  therefore,  with  a  five  cm  thick,  one  meter 
long  1 10  cm  radius  detector  composed  of  cadmium  tunstate  with  a  density  of  7  g/cm^.  An 
air  filled  cell  was  placed  immediately  behind  the  detector  cell  conforming  to  its  shape.  The 
number  of  tracks  entering  the  detector  cell  and  the  adjacent  cell  were  compared  (see 
Appendix  B,  Section  n,  “photon  activity  in  each  cefl”).  With  tiie  exception  of  some 
minimal  backscatter,  the  tracks  entering  the  air  filled  cell  could  onty  come  fi-om  the  detector 
cell,  tiius,  the  absorption  percent  could  be  determined  for  the  composite  source.  The 
absorption  of  the  detector  was  99.7%. 

6.3.3  Exposure  Tallies  and  Spectra. 

With  the  information  fi'om  the  above  tests,  a  complete  set,  consisting  of  32  beam 
positions,  was  run  with  MCNP  for  the  CT  unit  in  the  X-ray  room  (with  1/16”  lead 
shielding)  as  described  by  Metzger  (Metzger,  1993).  Appendix  D,  contains  the  exposure 
tallies  by  beam  position  and  normalized  average  special  plots  (with  the  roen^en  conversion 
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removed).  The  beam  positions  are  labeled  “Pt  XX”  where  XX  is  the  truncated  beam 
position  ai^e,  i.e.  11.25  degree  beam  angle  =  Pt  11,  33.75  degree  beam  angle  =  Pt  33. 

The  exposure  tally  at  each  detector  from  each  beam  position  (with  error  bars)  and  tiie 
mean  (with  its  error)  are  shown  in  Figures  D.l  throu^  D.7,  Figures  D.8  throu^  D.IO 
show  average  spectral  plots  for  each  detector  location.  The  source  detector  specta  shows 
that  tile  source  spectrum  was  maintained  at  30  cm.  These  exposure  and  spectral  results  for 
the  X-ray  room  could  be  generalized  to  a  CT  positioned  parallel  to  the  walls  assuming  the 
results  were  appropriate^  wei^ted  to  account  for  room  size  and  variations  in  shielding. 
From  these  results,  some  general  trends  were  observed  showing  the  combined  effects  of 
detector  and  beam  position. 

6.3.3. 1  Tidfy  Detector  Location. 

Figure  6.4  shows  a  representative  scatter  survQ'^  for  the  GE  ffiSpeed 
Advantage  CT  scanner  (GE,  1995).  As  expected,  tiie  dose  is  not  isotropically  distributed. 
The  same  general  behavior  was  expected  in  tiie  MCNP  runs.  The  MCNP  results  show 
similar  behavior  about  the  on-axis  direct  be<un  position  or  all  positive  and  negative  sets  of 
results.  Next,  althou^  the  negative-x  and  positive-x  axis  detector  mean  exposures  are 
approximately  equal  (Appendix  D,  Figures  D.l  and  D.2),  the  mean  exposure  of  the 
positive-y  axis  detector  is  greater  than  the  negative-y  axis  detector  (Figures  D.3  and  D.4). 
This  is  due  to  the  positive-y  detector  being  located  closer  to  the  source  than  the  negative-y 
detector.  The  positive  and  negative-y  axis  results  display,  however,  the  same  behavior,  i.e. 
no  real  dependence  on  source  position.  The  floor  mean  exposure  is  greater  than  the  ceiling 
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mean  exposure,  again,  because  it  is  much  closer  to  the  source.  The  floor  and  ceiling 
results,  like  the  x-axis  results,  display  a  very  marked  dependence  on  source  position,  but 
again,  they  are  similar  about  the  on-axis  direct  beam  position  (Figures  D.5  and  D.6). 

Figure  D.7  shows  die  unifomiity  in  die  source  detector  exposure  as  a  function  of  beam 
position. 

Although  an  attempt  was  made  to  make  the  geometty  identical  for  all  beam 
positions,  problems  were  encountered  with  die  90**  and  270**  source  beam  positions.  The 
geometiy  utility  within  MCNP  for  these  two  beam  positions  was  detecting  a  surface  within 
the  CT  detector  that  caused  it  to  consider  a  portion  of  the  detector  to  be  undefined,  and 
thus  ‘‘lose”  particles  in  this  undefined  area.  To  clear  this  fault,  the  detector  for  only  these 
two  points  was  split  into  two  adjacent  cells  with  the  same  total  volume,  position  and  shape. 
This,  in  turn,  produced  another  undefined  region  where  the  gantry  shielding  and  detector 
met  To  solve  this  problem,  die  gantiy  shield  for  both  beam  positions  was  changed  fi:‘om  a 
65  cm  inner-radius  cylindrical  shell  to  a  68  cm  inner-radius  cylindrical  shell.  It  appears  that 
increasing  the  radius  of  the  gantiy  shielding  allowed  an  opening  through  which  particles 
leaked.  Since  MCNP  forced  these  changes,  diey  could  not  be  immediateh^  corrected.  In 
addition,  note  that  die  ceiling  and  floor  exibit  very  high  central  peaks.  It  was  suspected  that 
a  problem  in  MCNP  caused  only  these  irr^^ilar  peaks.  Source  positions  were  anatyzed 
O.S**  and  0.01**  fi'om  the  90**  and  270**  beam  positions.  The  results  firom  diese  additional 
runs  more  closed  matched  the  beam  position  dependance  seen  for  the  positive-x  and 
negative-x  talfy  results.  It  is  likefy,  therefore,  that  die  exposure  results  for  90** 
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Figure  6.  4  Typical  scatter  survey  of  the  GE  HiSpeed  Advantage  scanning  a  body 
phantom,  units  are  mrem/scan,  but  should  be  used  for  comparison  as  relative.  The 
values  would  change  as  the  phantom  is  changed,  but  characteristic  contours  remain 
the  same. 
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and  270°  beam  positions  for  the  floor  and  ceiling  talfy  detectors  are  unreliable,  thus,  they 
were  removed  where  suspect.  Where  possible,  the  0.5  and  0.01  results  were  inserted. 


6.3.3.2  Summary. 

In  general,  the  x-axis,  floor  and  ceiling  exposure  tallies  in  the  remaining 
anafysis  configurations  are  expected  to  be  strongil^  dependent  on  beam  position, 
particulaify  diose  on  and  near  the  on-axis  direct  beams.  The  y-axis  e}q)osure  tallies  will  be 
much  more  uniform.  Since  anomalous  results  appear  at  exac%  the  90°  and  270°  beam 
positions,  they  were  deleted  where  suspect.  More  in-depth  analysis  of  exposure  and 
spectral  behavior  follows  in  Chapter  7  with  die  discussions  of  the  other  configurations. 
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7.  Results 


Three  different  room  configurations  were  modeled:  inside  and  outside  the  CT  suite 
at  WPMC,  and  outside  the  GE  minimum  recommended  CT  suite.  Both  a  head  and 
abdomen  scan  were  modeled  inside  the  suite  at  WPMC.  The  results  fi'om  inside  were 
compared  to  the  results  fi’om  film  badges  that  were  placed  in  the  suite  for  a  two  week 
period.  Head  and  abdomen  scans  were  modeled  for  the  dose  outside  the  CT  suite  at 
WPMC  with  1/16”  and  1/8”  lead  wall  shielding.  Only  a  head  scan  in  the  GE  minimum 
recommended  room  with  1/16”  lead  wall  shielding  was  model,  with  die  WPMC  results  for 
abdomen  and  1/8”  shielding  extended  to  the  GE  results. 

7.1  Validation;  Film  Badges  Inside  WPMC  CT  Suite 

To  validate  the  modeling  effort  and  MCNP’s  results,  radiation  film  badges  were 
exposed  during  normal  operation  within  WPMC’s  CT  suite  fi'om  14  to  28  August  1995. 
The  reported  dose  equivalent  (in  mrem,  for  deep,  eye  and  shallow  tissue)  was  compared  to 
the  results  of  an  MCNP  run  modeling  the  location  of  the  film  badges  for  bodi  a  head  and 
abdomen  slice.  Twelve  film  badges  were  placed  on  the  source,  bo^  and  x  and  y-room 
axes  (one  at  each  location,  8  total)  and  on  the  ceiling  and  floor  (two  at  each  location,  4 
total),  see  Figure  5.4.  The  duplicate  badges  on  the  floor  and  celling  were  used  to  verify  the 
precision  of  the  reported  badge  measurements.  To  speed  computation,  the  suite  was 
modeled  with  no  walls.  This  was  valid  since  test  runs  had  shown  that  backscatter  fi’om  the 
walls  into  the  room  was  less  than  1%.  The  results  fi’om  inside  the  suite  were  used  as  a 
standard  for  later  analyses. 
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7.1.1  Comparison. 


Before  con:q>ariiig  the  film  badge  results  to  the  MCNP  results,  note  some  general 
behavior  that  relies  solefy  on  distance  fi'om  isocenter.  The  floor  badges  were  about  one 
meter  below  isocenter  and  the  ceiling  badges  about  two  meters  above.  The  floor  dose, 
therefore  is  expected  to  be  about  four  times  huger  than  the  ceiling  dose  (using  the  inverse 
squared  law  for  distance).  Since  the  positive  and  negative  sides  of  die  rooms  were  about 
equidistant  from  isocenter,  the  badge  results  for  die  positive  and  negative  source  should  be 
about  equal.  The  same  should  be  seen  for  the  x-room  axis  badges.  The  y-room  and  y- 
body  axes  results  cannot  be  compared  simply  by  distance  since  the  presence  of  the  patient 
affects  the  results.  l^Tth  the  exception  of  the  negative  x-room  axis  badge  and  the  higher 
floor  badge  results,  these  relationships  were,  in  general,  as  expected. 

Table  7.1  shows  the  results  fi’om  the  MCNP  runs  for  the  head  and  abdomen  scans 
using  Eq.  (6.4)  to  determine  effective  dose.  The  average  weekly  current  workloads  were 
3980  mA-min  per  week  and  10913  mA-min  per  week  for  the  head  and  abdomen  scans, 
respectively,  which  differed  fi’om  the  4534  mA-min  per  week  and  11468  mA-min  averages 
seen  in  Table  6.2.  Table  7.2  presents  a  comparison  of  the  badge  to  MCNP  results.  Note, 
first,  that  although  tiie  two  measured  ceiling  doses  were  equal,  the  floor  doses  differed 
significantly.  This  difference  was  used  to  assign  a  15%  uncertainty  in  the  precision  of  the 
badge  results.  The  wall  badges  were  located  about  four  feet  above  the  floor  with, 
unfortunatcty,  the  potential  of  being  blocked  sometime  during  tiie  15  day  cTqiosurc  period, 
by  the  carts,  tables  and  miscellaneous  accessories  located  around  tiie  room.  It  was  reported 
by  tiie  medical  technician  that  the  badges  on  the  negative  and  positive  x  walls  were  likety 


7-2 


Table  7.  1  MCNP  Results  from  Inside  the  CT  Suite  at  WPMC 


Eiposive  by  Scan  Rs'*  particle^ 

Efifectivs  dose,  nrenf 

Head 

Flradional 

Ehtr 

Abdomen 

Fl^ona 

IFhtr 

Abdonni 

Totii 

Saurcie^4x 

4.77B-18 

0.0223 

3.64E-18 

0.0105 

167 

34.9 

Souroe^  -x 

5.9®-18 

0.0264 

4.mi8 

0.0111 

20.8 

45.3 

Bodjy,-ty 

224B-17 

0.0076 

7.54B.18 

0.0128 

783 

722 

15015 

Bo(fy,-y 

27ffi-18 

0.0545 

208E-18 

0.0072 

9.6 

20.0 

29.6 

Boon\+x 

3.85B-17 

0.0102 

3.63B-17 

0.0049 

134.6 

347.5 

482.1 

Room,  -X 

274B-17 

0.0096 

1.18B-17 

0.0058 

95.6 

1128 

2084 

Roon^-ty 

258E-17 

0.0091 

1.33B-17 

0.0103 

90.1 

127.7 

217.8 

Room,-y 

3.37B-18 

0.0455 

4.8^18 

0.0153 

11.8 

46.3 

581 

Geiliiig 

25SE-18 

0.0618 

234B-18 

0.0319 

9.0 

224 

3L4 

Floor 

21^17 

0.0360 

1.2S-17 

0.0266 

74.3 

119.3 

193.6 

dose = Kdatne  x  0. 1068(GyinA4niny0.0086S(R/Gy)  x  anetit(ni/Vniin/vd0  x  2(vil^ 


X  1000(innsnyien]^  x  0.01(iBtnl^ 

Table  7.  2  Comparison  of  Film  Badge  and  MCNP  Results  Inside  WPMC  CT  Suite 


Bad^ 

Nundier 

Locatkm 

Badge 

Effective  Dose,  MCNP  Effective 
nvnn  Dose,inrmi 

Badge/ 

MCNP 

11 

Source  +x 

30 

51.5 

MSM 

8 

Source,  -x 

30 

66.1 

■B 

5 

Body, +y 

150 

150.5 

1.00 

4 

Body,-y 

20 

29.6 

0.68 

1 

Room,  +x 

90 

482.1 

0.19 

2 

Room,  -X 

0 

208.4 

- 

6 

Room,  +y 

230 

217.8 

1.06 

3 

Room,  -y 

60 

58.1 

1.03 

9 

Ceiling 

60 

31.4 

1.91 

12 

Ceiling 

60 

31.4 

1.91 

7 

Floor 

320 

193.6 

1.65 

Floor 

250 

193.6 

1.29 

Effective  Dose=Relative  x  Aver  Air  K!aina(0. 1068  Gy/mA  mmy0.00868(R/Gy) 
xCurr€nt(3980  or  10920  niAmin)x2vksxl000(iiirestn/raii)x  0.0  l(reni/R) 


7-3 


blocked  for  some  or  all  of  the  IS  days.  Note  that  the  negative  x-room  axis  badge  dose  was 
zero  and  that  the  positive  x-room  axis  badge  and  MCNP  doses  differ  significandy.  With 
these  exceptions,  the  reported  badge  doses  were  within  a  factor  of  two  of  die  MCNP 
results,  some  greater,  some  lesser. 

7.1.2  Baseline  Results. 

The  MCNP  results  from  inside  the  CT  suite  in  WPMC  suite  were  used  as  a 
standard  against  vdiich  the  remaining  anafyses  were  compared  because  they  show  the 
scattered  spectrum  and  dose  as  a  function  of  onfy  die  CT  unit.  Appendices  E  and  F 
contain  the  exposure  (by  beam  position)  and  spectral  plots  obtained  for  the  head  and 
abdomen  scans,  respectively.  The  head  scan  results  are  used  to  describe  the  exposure  and 
spectral  differences.  They  are  then  compared  to  the  abdomen  exposure  and  spectral 
results. 


7. 1.2.1  Exposure. 

The  e;q)osure  tallies  as  a  function  of  beam  position  for  each  tally  detector 
location  illustrate  two  effects  on  the  exposure  behavior  as  the  source  is  rotated:  (a)  the 
position/direction  of  die  beam  relative  to  the  talfy  detector,  i.e.  the  beam  pointing  direct^  at 
or  away  from  the  talfy  detector  and  (b)  the  cross  sectional  area  of  the  target  (head)  being 
scanned.  Figure  7.1  shows  a  sample  e^qiosure  plot  with  error  bars,  mean  and  mean  error. 

The  x-z  plane  exposures  (+/-  x-source  axes.  Figs.  E.l  and  E.2,  ceiling,  Fig.  E.9, 
and  floor.  Fig.  E.  10)  illustrate  effect  (a).  They  show  a  dramatic  increase  in  contribution 
from  beam  positions  directed  at  the  talfy  detector.  The  twin-peaked  response  seen  in  the 
source  axes  and  floor  exposure  are  caused  by  photons  scattering  around  the  CT  detector 
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Figure  7. 1  Example  of  exposure  tally  (by  beam  position)  plot  Direct  beam  from 
180"  beam  position. 


with  minimal  absoiption  by  die  gantiy  shielding  (see  Figure  7.2).  The  ceiling  exposure 
behavior  differs  from  the  source  axis  and  floor  exposures  with  a  central  peak.  Since  the 
ceiling  is  much  fardier  away  from  the  source  than  the  floor,  the  scattered  photons  have  a 
greater  opportunity  to  scatter  back  toward  the  talfy  detector  (see  Figure  7.3). 

The  effect  of  target  size  is  seen  in  the  positive  and  negative  y  exposures  (Figures 
D.5  and  D.6).  The  Pt  90  and  270  beam  positions,  directly  below  and  above  the  head, 
respectively,  strike  the  narrowest  target  cross-section,  thus  less  head  scatter  and  a  lower 
exposure  (see  Figure  S.5).  Recall  tiiat  99.7%  of  tiie  beam  missing  the  head  is  absorbed  by 
the  CT  detector.  The  negative  y-bo^  axis  exposure  does  not  have  as  large  of  a  variation  as 
the  beam  rotates  since  the  majority  of  the  body  is  along  that  axis  which  is  an  additional 
source  of  absorption. 


7-5 


Figure  7. 2  Possible  scattering  mechanisms  producing  peaks  from  off  axis  beam 
positions  :  (1)  shows  photons  stopped  by  the  CT  detector.  (2)  shows  a  direct 
contribution  to  tally  detector  A  from  the  axis  beam.  (3)  shows  a  possible  scattering 
path  that  would  result  in  a  greater  exposure  at  tally  detector  location  B  (its  position 
shown  relative  to  the  source  location)  from  off-axis  beam  positions. 


Ceiling  Detector 


Figure  7.  3  As  photons  travel  away  from  the  source,  they  are  scattered  across  a  wider 
area  with  the  CT  detector  blocking  some  of  the  scattered  photons.  Since  the  ceiling  is  a 
greater  distance  from  the  source,  the  photons  are  more  abie  to  scatter  back  towards  a 
tally  detector  located  opposite  the  source. 
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The  exposure  at  talfy  detectors  located  on  the  room  axes  (Figs.  £.3,  E.4,  E.7  and 
E.8)  inside  the  CT  suite,  like  the  y-room  and  y-source  axis  results,  shows  only  tai^et  size 
effects.  They  show  means  and  variations  in  amplitude  similar  to  the  y-room  and  y-source 
axis  results.  The  relative  mean  exposures  conq)are  well  with  the  typical  scatter  survey 
(Figure  6.4). 


7. 1.2.2  Spectral  Results. 

D.  Simpkin  (Simpkin,  1995),  identified  the  lack  of  patient  scatter  spectra 
fi'om  CT  units  since  available  test  eququnent  was  unable  to  measure  the  hig|h  flux  inside  the 
room  (Simpkin,  1995).  MCNP  bins  die  talty  by  user  specified  energy  intervals  thus 
providing  the  scattered  spectrum.  The  conversion  to  Roentgens  was  removed  fi'om  the 
average  spectra  in  Appendix  £  and  the  resulting  spectra  normalized  to  one.  First,  observe 
that  aldiou^  the  individual  tally  detector  locations'  exposure  tally  means  and  behavior  are 
quite  different  inside  the  CT  suite,  they  all  have  very  similar  spectra  (Figs.  E.12  through 
E.  17)  fi'om  10  to  70  keV.  Above  70  keV,  however,  the  x-z  plane  talty  detectors  spectra 
exhibit  a  larger  value  at  the  higher  energies.  In  addition  to  the  mean  spectrum.  Figs.  £.11, 
E.12,  £.16  and  £.17  contain  additional  spectra  fi’om  beam  positions  directed  at  or  nearly  at 
the  talty  detector.  They  exhibit  a  noticeable  increase  in  the  h^er  energy  ranges, 
particularty  those  adjacent  to,  but  not  on,  the  axis.  (All  of  the  other  beam  positions 
mimicked  the  average.)  Recalling  the  scattering  mechanisms  described  above  (Figures  7.1 
and  7.2),  the  higiher  energy  humps  are  cleaity  from  photons  passing  throu^  the  CT 
detector  and  gantry  with  litde  attenuation. 
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Next,  observe  the  sharp  spikes  at  24  and  28  keV  in  the  average  and  individual 
specta  as  seen  in  Fig.  £.  1 1  or  E.  12.  These  are  fluorescent  peaks  from  the  cadmiiim 
tungstate  CT  detector.  The  24  keV  spike  is  from  the  Li  to  K  shell  transition  of  cadmium 
and  the  28  keV  spike  is  probably  from  summed  X  rt^rs  resulting  from  the  de-excitation 
cascade.  Sharp  spikes  are  also  seen  in  the  lowest  energy  bin  containing  0-10  keV.  They 
are  the  fluorescent  peaks  from  the  intermediate  shell  transitions  and  from  low  Z  elements  in 
the  problem,  ptimaiify  aluminum  (used  for  the  gantry  shields).  Remnants  of  the  source 
spectrum’s  peaks  at  60  and  68  keV  are  also  seen. 

Last,  a  spectral  shift  is  referred  to  as  either  hardened  or  scfiened.  Removing  the 
lower  energy  components  results  in  a  hardened  spectrum ,  fliat  is,  one  with  a  greater 
average  energy,  and  removing  the  higher  energy  components  results  in  a  softer  spectrum, 
one  with  a  lower  average  energy.  Figure  7.4  shows  that  the  source  spectrum,  upon 
scattering  off  the  phantom  was  softened  slightly,  with  the  maximum  shifted  approximately 
lOkeV. 


7. 1.2.3  Head  vs  Abdomen. 

Appendix  F  contains  the  exposure  and  spectral  plots  from  the  MCNP  run 
inside  the  CT  suite  at  WPMC  for  an  abdomen  scan.  Different  exposure  behaviors  were 
expected  from  the  abdomen  because  the  target  is  considerably  different. 

First,  althou^  the  MERD  phantom’s  abdomen  and  head  x-z  plane  cross  sections 
are  both  elliptical,  the  major  and  minor  axes  are  opposite.  The  variation  with  source  angle 
seen  with  the  head  slice  from  relative  target  fliickness  should  be  opposite  that  of  the  head 
scan.  This,  however  was  not  seen  (see  Figs.  G.3  and  H.3).  Second,  the  abdomen  is  much 
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Figure  7.  4  Scattered  spectrum  vs  source  spectrum. 

thicker  than  the  head  in  the  x  direction  (14  cm  vs  40  cm).  The  contribution  to  the  x-axis 
taify  detectors  opposite  the  0**  and  180°  beam  positions  should,  therefore,  be  reduced. 
Figure  7.5  shows  a  representation  of  the  x-source  axis  tally  detectors  for  the  abdomen  scan. 
The  reduction  can  be  seen  in  Figs.  F.  1,  F.2,  F.9,  and  F.  10  when  compared  to  the 
analogous  head  scan  figures  in  Appendix  £.  Tissue,  however,  a  primarily  low-Z  material, 
is  not  as  strong  absorber  as,  say,  lead  or  cadmium  tungstate,  so  there  is  not  a  large 
reduction. 

Third,  the  surface  of  the  target  being  irradiated  is  considerabfy  closer  to  the  0°  and 
1 80°  beams  and  thus  closer  to  the  tally  detectors  behind  fiiese  beams.  Recalling  that  a 
fluence  tally  computed  in  MCNP  is  a  function  of  distance  and  weight,  diis  surface  radiation 
may  contribute  more  to  the  exposure  tally  \^en  the  beam  is  directed  awcy  from  tire  tally 
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detector  (see  Figure  7.S).  Associated  with  this  is  a  possible  increased  effect  from  die 
rotation  of  the  phantom  and  CT  unit  about  the  z  axis.  The  position  where  the  phantom 
receives  maximum  radiation  will  be  closer  to  die  negative  y  talfy  detector  when  the  beam 
position  is  0”  and  converse^,  closer  to  the  positive  y  tally  detector  when  the  beam  position 
is  180°.  The  results  illustrated  in  Figure  1.5  for  the  positive  x-  and  y-room  axes  exposures 
are  seen  in  Figs.  F.3  and  F.7:  die  maximum  occurred  when  the  beam  was  directed  away 
from  die  tally  detector.  The  second  part  of  this  e?q>ectation  is  seen  in  Figs.  F.7  and  F.8. 
The  maximum  is  seen  at  0°  for  the  n^tive-y  tally  detector  and  at  180°  for  die  positive-'y 
talfy  detector. 


Detector 

r’ 
r 


Figure  7.  5  Effect  on  exposure  as  a  function  of  distance  to  radiated  area. 
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Fourth,  the  beam  positions  adjacent  to  the  on-axis  positions  (  Le.  Pts  11  and  22 
adjacent  to  Pt  0)  should  still  exhibit  die  sharp  rise  in  exposure.  The  maximum,  however, 
may  come  from  a  position  fiirdier  from  the  on-axis  position  than  noted  with  the  head  scan 
because,  referencing  Figure  7.2,  to  maintam  the  neaify  straight  line  sequence  of  interactions 
from  source  to  tally  detector,  the  source  would  need  to  be  further  from  the  on-axis  position 
to  get  around  the  wider  abdomen.  Figs.  F.l,  F.2,  F.9  and  F.  10  show  this  change.  The 
exposure  from  the  off-axis  beam  positions  closest  to  die  on-axis  positions  (Pts  11,  168,  78 
and  258,  respective^)  are  reduced. 

Fifth,  the  phantom  is  approximately  centered  longitudinally  on  isocenter  during  an 
abdomen  scan.  Differences  noted  before  between  the  positive-  and  negative-y  axes  (room 
and  bo^)  should  be  reduced.  The  negative-y  tally  detectors,  however,  are  still  about  twice 
as  far  from  isocenter  as  the  positive  talfy  detectors,  so,  greater  average  exposures  from  the 
positive  direction  are  still  expected.  The  positive  and  negative  mean  exposures  are  closer 
than  seen  with  the  head  scan.  In  addition,  the  change  in  exposure  as  die  beam  rotates  is 
more  obvious  in  the  negative-y-body  and  room  talfy  detector  tallies  (Figs.  F.6  and  F.8). 

All  of  the  above  differences  were  present  simultaneously  for  all  talfy  detectors. 

Some  dominated  leaving  the  others  virtually  unnoticeable. 

A  problem  in  the  CT  model  was  found  as  the  data  was  being  analyzed.  A  fixed 
fan-beam  an^e  of  35®  was  used  throughout  the  anafysis.  However,  along  its  major  axis, 
the  abdomen  subtends  a  slightly  wider  angle  (37°)  than  the  fan-beam  angle.  This  means 
that  at  die  90°  and  270°  beam  positions,  no  portion  of  the  beam  will  go  past  the  phantom 
directfy  to  the  CT  detector.  The  most  noticeable  effect  expected  is  a  reduction  in  the 
contribution  to  e;q>osure  from  the  beam  positions  adjacent  to  the  on-axis  beam  positions 
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(Pts  1 1  and  348  for  die  R  0  on-axis  beam  position).  Another  effect  may  be  a  flattening 
near  the  90"  and  270"  points  in  exposure  variation  for  the  ceiling  and  floor  tally  detectors. 
This  mistake  is  not  e^qpected  to  contribute  significant  error  to  die  results. 

The  differences  in  the  spectra  between  the  head  and  abdomen  scans  for  all  talty 
detector  locations  was  insignificant  as  can  be  seen  by  comparing  Figs.  F.  1 1  througjh  F.  1 5  to 
their  counterparts  in  Figs.  E.ll  throu^  E.IS.  Therefore,  the  spectra  outside  the  walls 
should  also  be  conqiarable  for  the  head  and  abdomen  scans.  The  similarity  in  the  spectra 
leads  to  the  supposition  that,  when  the  shielded  walls  are  added  to  the  scenmo,  similar 
changes  to  the  head  and  abdomen  scans  should  be  noted.  For  example,  the  difference 
between  doses  at  a  given  talfy  detector  fi'om  a  head  scan  to  an  abdomen  scan  will  decrease 
by  tile  same  factor  for  the  1/16”  lead  shielded  room  and  the  1/8”  lead  shielded  room  at  the 
same  tally  detector.  That  is,  the  head-to-abdomen  dose  ratio  will  be  the  same  outside  either 
room. 


7.1.3  Modifications. 

The  results  fi’om  inside  the  CT  suite  at  WPMC  drove  modifications  in  subsequent 
runs.  The  most  significant  change  was  eliminating  the  anafysis  of  some  tally  detector 
locations  for  certain  beam  positions.  The  exposure  was  expected  to  decrease  across  the 
board  as  the  shielded  walls  were  added  because  of  the  similarity  of  their  spectra.  Note,  fin* 
example,  the  exposwe  by  beam  position  plot  for  the  positive  x-source  axis  tally  detector 
(Fig.  E.  1).  The  onfy  significant  contributions  to  the  total  exposure  were  fi*om  the  beam 
positions  fi’om  R  0  to  R  101  and  fi’om  R  258  to  R  348.  Thus,  tiie  positive  x-source  axis 
exposure  was  not  computed  for  the  beam  positions  fi-om  R 1 12  to  R  247.  When  the 
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exposure  displayed  uniform  behavior  at  all  beam  positions,  like  the  negative  y-body  axis 
exposure  (Fig.  E.6),  only  a  few  beam  positions  were  anafyzed.  Conservation  of  computer 
resources  was  the  primaiy  motivation  for  this  reduction  in  e?^sure  computations.  Table 
7.3  lists  the  exposures  that  were  computed  for  each  beam  position.  Comparison  wi&  the 
exposure  plots  in  Appendix  E  explain  the  choices  that  were  made.  Prior  to  completion  of 
the  runs  inside  the  suite,  runs  were  started  for  the  analysis  outside  the  suite,  thus,  extra 
exposure  data  points  are  included  in  some  plots.  As  described  in  Section  6.1.3,  the 
exposure  tallies  and  errors  that  remained  were  weighted  appropriate^. 

In  addition,  a  need  was  determined  for  additional  beam  positions  near  the  e?q)osure 
maximums  to  better  define  e^qiosure  behavior  and  thus  e?qK)sure  mean.  However,  time 
constraints  caused  by  limited  computer  resources  prevented  tiieir  analysis  concurrent  with 
the  results  presented  here. 

The  results  also  raised  questions  about  exposure  behavior  between  and  around  the 
selected  ta%  detector  locations.  The  concern  whether  the  true  maximums  were  indeed  at 
the  source  and  room  axes  locations  was  investigated  by  placing  tally  detectors  at  50  cm 
intervals  in  the  positive-  and  negative  y  direction  in  line  with  the  existing  tally  detectors. 

This  was  accomplished  for  the  anafysis  of  the  CT  suite  in  WPMC  for  the  positive  and 
negative  x  walls.  Those  results,  shown  in  Figure  7.6,  indicate  that  there  is  indeed  a  “hot 
spot”  not  sampled  exacdy  at  the  source  or  room  axis  tally  detectors.  The  negative  x  wall 
displayed  a  maximum  20  cm  from  the  source  axis  talfy  detector  and  the  positive  x  wall 
displayed  a  maximum  10  cm  from  the  source  axis  tally  detector.  The  existence  of  these  hot 
spots  may  indicated  the  existence  of  others  titat  were  not  sought  out. 
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Table  7.  3  Modified  Exposure  Tallies  by  Beam  Position 
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7.2  WPMCCT  Suite 

The  main  objective  of  this  thesis  was  to  determine  the  dose  outside  the  walls  of  the 
CT  suite  at  WPMC  to  determine  if  its  shielding  met  die  recently  revised  requirements  for 
non-occupational  dose  in  10CFR20.  The  following  results  deal  directly  with  that  issue. 
The  suite  was  modeled  alternately  with  1/8”  lead  shieldii^  and  1/16”  lead  shielding  in  die 
walls.  The  Monte  Carlo  simulation  was  conducted  for  both  a  head  scan  and  an  abdomen 
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Figure  7.  6  Exposure  outside  of  and  along  the  (a)  positive  x  wall  and  (b)  negative  x  'wall.  They  show  maxima  at  (a)  +150 
cm  from  the  wall  axis  and  (b)  -100  cm  from  the  wall  axis.  The  lower  plot  is  an  expanded  view  of  the  upper  plot 


scan  at  10  locations  outside  the  walls  for  selected  beam  positions.  The  weighted  average 
exposure  for  the  head  and  abdomen  scans  were  totaled  to  obtain  the  dose  at  each  talfy 
detector  location. 

7.2.1  1/16’’ Lead  Shielding. 

The  CT  suite  at  WPMC  was  modeled  with  1/16”  lead  shieldii^  for  both  head  and 
abdomen  scans.  The  exposure  plots  (by  beam  position)  and  normalized  spectral  results 
(with  Roen^en  conversion  removed)  are  in  Appendices  G  and  H.  The  exposure  behavior 
of  the  head  and  abdomen  scans  are  as  expected:  diminished  considerabfy  in  magnitude  and 
starting  to  lose  some  structure  as  a  function  of  source  positions  when  conq)ared  to  die 
exposure  tallies  from  inside  the  suite.  Recall  that  the  largest  exposures  (on-  and  near-axis 
beam  positions  directed  at  the  tally  detector)  had  a  hi^  energy  component  in  dieir 
spectrum.  The  lead  shielding  absorbed  most  of  that,  and,  since  the  total  exposure  tally  is 
effectively  an  integration  of  the  spectrum,  the  peaks  with  the  hi^  energy  component  were 
reduced  relatively  mcnre  than  the  other  beam  positions  in  the  scan,  hi  particular,  compare 
Pt  22  for  the  positive  x-source  axis  in  Fig.  E.l  with  the  same  beam  position  in  Fig.  G.l. 
More  than  half  of  its  spectrum  was  above  80  keV  for  this  particular  point  (Fig.  E.  1 1), 
which  has  all  been  removed  by  die  room  shielding.  A  similar  comparison  of  Pt  157  in 
Figs.  E.2  and  G.2  reveal  the  same  effect  as  do  Pts  78  and  258  for  the  ceiling  and  floor 
exposures,  respectively,  in  Figs.  E.9,  E.IO,  G.9,  and  G.IO.  The  same  effects  are  noted  with 
the  abdomen  exposure  plots.  The  smooth  variation  seen  in  the  y-axes  tally  detector 
exposure  plots  is  starting  to  disappear  with  die  addition  of  the  lead  shielding  (compare  Fig. 
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E.3  to  Fig  G.3  and  E.5  to  G.5).  The  exposure  plots  showing  a  uniformity  for  all  beam 
positions  have  maintained  their  uniformity  (compare  Fig.  E.8  to  Fig  G.8).  In  addition,  note 
that  as  the  mean  went  down,  so  did  die  spread  in  the  individual  exposure  values.  Compare 
Fig.  E.5  with  a  mean  of  2.3  x  10'^^  and  spread  of  about  1  x  10’*^  to  Fig.  G.5  widi  a  mean 
of  1.5  X  10*“  and  a  spread  of  about  1  x  10*“. 

Dramatic  spectral  differences  are  noted  between  the  spectra  outside  and  inside  the 
CT  suite.  This  was  expected  since  lead  shielding  is  chosen  to  reduce  die  higher  energy 
photons,  i.e.  modify  the  spectrum.  First,  note  that  the  spectra  for  the  head  and  abdomen 
scans  were  indistinguishable  (Figs.  G.ll  through  G.15  and  H.11  through  H.15).  Next, 
observe  that  the  wall  spectra  (G.  1 1  diorough  G.  14)  are  abruptly  cut  off  near  90  keV. 
Comparison  with  the  attenuation  cross-section  for  lead  (Figure  2.7)  shows  die  K-edge  at 
88  keV  which  effectively  attenuates  the  entire  CT  spectrum  above  that  point. 

7.2.2  1/8”  Lead  Shielding. 

The  CT  suite  at  WPMC  was  next  modeled  with  1/8”  lead  shielding  for  head  and 
abdomen  scans.  The  exposure  and  spectral  results  are  in  Appendices  I  and  J.  The 
exposure  plots  for  the  head  and  abdomen  are  as  expected;  further  diminished  in  magnitude 
and  structure.  Once  again,  the  laigest  individual  exposures  (from  beam  positions  directed 
at  or  nearly  at  die  tally  detector)  which  had  a  hi^  eneigy  component  m  their  spectrum 
inside  die  room  are  reduced  relatively  more  than  the  other  beam  positions  in  the  scan. 
Compare  now  Pt  11  for  the  positive  x-source  axis  in  Fig.  G.l  with  die  same  beam  position 
in  Fig.  1.1.  A  similar  comparison  of  Pt  168  in  Figs.  G.2  and  1.2  reveal  the  same  effect. 
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The  slight  variation  left  in  the  1/8”  lead  shielded  room  seen  in  the  y-axes  talfy  detector 
exposures  has  almost  completely  disappeared  (compare  Fig.  G.3  to  Fig  1.3  and  G.S  to  1.5). 
The  exposures  that  were  uniform  for  all  beam  positions  have  again  maintained  their 
uniformity  (compare  Fig.  E.8  to  Fig  G.S). 

The  spectral  differences  between  die  1/16"  and  1/8"  shielded  rooms  are  not  as 
pronounced  as  those  noted  between  the  spectra  outside  with  1/8”  lead  shielding  and  inside 
the  CT  suite,  however,  there  are  some  differences.  First,  however,  note  that  the  spectra  for 
the  head  and  abdomen  scans  were  again  in^stinguishable  (Figs.  1. 1 1  through  1 1 5  and  J.  1 1 
through  J.  IS).  Each  spectrum  now  displays  a  softer  spectrum  compared  to  the  1/16”  lead 
shielding  spectra.  The  energies  near  88  keV  have  been  reduced  more  extremely,  allowing 
the  energies  between  40  keV  and  70  keV  to  appear  huger  by  comparison. 

7.2.3  Dose. 

Table  7.4  contains  the  calculated  dose  outside  the  CT  suite  at  WPMC  for  1/16” 
and  1/8”  lead  shielding  for  the  head  and  abdomen  scans  and  Oieir  total.  The  dose  was 
calculated  using  Eqs.  (6.1),  (6.2)  and  (6.4).  Recall  that  the  annual  dose  limits  required  by 
10CFR20  were  1  mSv  (non-occupational,  continuous  exposure),  5  mSv  (non- 
occupational,  infrequent  exposure)  and  50  mSv  (occipational).  The  ratio  of  tissue  to  air 
mass  energy  absorption  coefficients  was  used  to  convert  of  these  results  in  units  of  Gray  to 
units  of  Sieverts  as  discussed  in  section  2.8.6.  Both  the  1/8”  lead  shielded  and  1/16”  lead 
shielded  rooms  were  modeled:  first,  for  corrrpleteness  and  second,  because  fire  CT  suite  at 
WPMC  has  1/8”  lead  shielding  in  the  walls  adjacent  to  public  areas  (negative  x  and  positive 
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y  walls)  and  1/16”  lead  shielding  in  the  walls  adjacent  to  the  occupational  areas  (positive  x 
and  negative  y  walls).  Therefore,  non-occupational  dose  is  indicated  in  the  positive  y  and 
negative  x-room  and  x-source  talfy  detector  locations.  The  occupational  dose  is  indicated 
in  the  negative  y  and  positive  x-room  and  x-source  tally  detector  locations. 

With  these  guidelines,  the  only  ta%  detector  locations  that  exceed  die  non- 
occupational  dose  limits  are  the  ceiling  and  floor.  Since  the  badge  results  showed  an 
underestimation  by  die  Monte  Carlo  simulation  inside  the  room,  one  of  two  conclusions 
must  be  drawn;  (a)  the  dose  above  the  ceiling  and  below  the  floor  is  unacceptable  by 
10CFR20  standards  or  (b)  the  concrete  shielding  was  insufficient^  modeled.  Recall  that 
six  inch  slabs  of  concrete  were  placed  in  the  ceiling  and  floor.  The  floor  and  ceiling  of 
WPMC,  however,  have,  in  addition,  two-inch  thick  waffle-weave  concrete  above  and 
below  the  slab.  As  a  worst  case  scenario  (and  to  conserve  computer  resources)  this 
additional  concrete  was  not  modeled.  It  would  have  contributed  to  die  floor  and  ceiling 
attenuation,  but,  the  exact  amount  is  not  presen%  known. 

Table  7.4  shows  that  the  occiqiational  dose  limit  of  SO  mSv  is  met  at  aU  talfy 
detector  locations.  The  floor  and  ceiling  exposures  are  significantly  higher  than  the  rest, 
but  still,  well  below  the  10CFR20  dose  limits. 

7.3  General  Electric  Minimum  Recommended  Room  Size 

The  product  data  sheet  accompanying  the  GE  HiSpeed  Advantage  Scanner  System 
recommended  a  13’  x  19.5’  (397  cm  x  596  cm)  exam  room  with  a  9’  (274  cm)  ceiling.  A 
drawing  of  the  configuration  showed  isocenter  183  on  fi'om  die  negative  x  wall  and  213 
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Table  7.  4  Results  for  the  CT  Suite  in  WPMC  for  1/16”  and  1/8”  Lead  Shielding 


WPMC, 

1/16"  Lead  Shielding 

Head 

Abdomen 

1  Total  i 

Detector 

Annual  Tissue 

Fractional 

Error 

Annual  Tissue 

Fractional 

Error 

Annual  Tissue 

Fractional 

Location 

Equivalent 
Dose,  mSv 

Equivalent 
Dose,  mSv 

Equivalent 
Dose,  mSv 

Error*^ 

Source,  -x 

0.453 

0.079 

0.797 

0.098 

1.251 

MEM 

Source,  +x* 

0.391 

0.078 

0.731 

0.086 

1.122 

gg 

Body,  +y 

1.481 

0.071 

1.324 

0.091 

2.804 

0.066 

Body,  -y^ 

0.125 

0.176 

0.142 

0.102 

0.267 

0.152 

Room,  -X 

2.787 

0.059 

4.281 

0.075 

7.069 

0.054 

Room,  +x^ 

2.571 

0.053 

5.250 

0.052 

7.821 

0.047 

Ceiling* 

0.787 

0.057 

1.043 

0.046 

1.830 

0.050 

Floor^ 

10.817 

0.035 

7.671 

0.032 

18.488 

0.031 

Room,  +y 

2.872 

0.059 

3.969 

0.071 

6.841 

0.054 

Room,  -y^ 

0.429 

0.117 

1.200 

0.086 

1.630 

0.102 

WPMC 

,  1/8"  Lead  Shielding 

Head 

Abdomen 

Total  1 

Detector 

Annual  Tissue 

Fractional 

Error 

Annual  Tissue 

Fractional 

Error 

Annual  Tissue 

Fractional 

Error 

Location 

Equivalent 
Dose,  mSv 

Equivalent 
Dose,  mSv 

Equivalent 
Dose,  mSv 

Source,  -x® 

mgam 

0.123 

0.123 

0.185 

0.096 

Source,  +x 

Mg 

0.122 

0.146 

0.182 

0.110 

Body,  +y® 

gg 

0.212 

0.121 

0.321 

0.085 

Body,  -y 

0.018 

0.124 

0.046 

0.128 

0.064 

0.096 

Room,  -X® 

0.160 

0.098 

0.317 

0.120 

0.477 

0.078 

Room,  +x 

0.171 

0.105 

0.369 

0.100 

0.541 

0.080 

Ceiling® 

0.740 

0.077 

1.285 

0.068 

2.024 

0.059 

Flooi* 

4.292 

0.033 

7.753 

0.029 

12.045 

0.025 

Room,  +y® 

0.178 

0.095 

0.355 

0.111 

0.533 

0.075 

Room,  -y 

0.033 

0.199 

0.085 

0.182 

0.118 

0.151 

^Detector  locations  applicable  to  occupational  dose. 

^Detector  locations  applicable  to  non-occupational  dose, 
c 

Error  is  wetted  by  usage  factors. 


cm  from  the  positive  y  wall  and  with  the  unit  rotated  21.86  degrees  from  tiie  positive  y- 
room  axis  towards  die  negative  x-room  axis.  Except  for  the  dimensions  and  angle  of 
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rotation,  the  room  is  identical  to  the  CT  suite  at  WPMC  with  1.4  cm  diywall  sheets  on 
either  side  of  the  1/16”  and  1/8”  lead  and  six  inch  concrete  slabs  at  the  floor  and  ceiling. 


7.3.1  1/16”  Lead  Shielding  with  Head  Scan. 

Appendices  J  and  K  contain  the  exposure  and  spectral  plots  for  the  minimum  room 
size  recommended  by  G£  for  die  IhSpeed  Advantage  CT.  Since  the  exposure  behaviors 
and  spectra  are  very  similar  to  those  seen  in  the  CT  suite  with  1/16”  lead  shielding,  a 
detailed  anafysis  win  not  be  presented  here. 

7.3.2  Extension  of  WPMC  Results. 

The  differences  in  average  dose  from  the  head  to  the  abdomen  scan  and  from 
1/16”  to  1/8”  lead  shielding  for  the  CT  suite  at  WPMC  were  extended  to  the  results  from 
the  GE  room  with  1/16”  lead  shielding  for  the  head  scan.  Figure  7.5  shows  two  paths  for 
obtaining  the  GE  room  results  for  1/16”  abdomen,  1/8”  head  scan  and  1/8”  abdomen  scan. 
Since  the  conversion  factors  are  constants,  bodi  paths  produce  the  same  results.  Table  7.5 
shows  the  ratios  seen  for  the  WPMC  suite  exposures  diat  were  applied  by  source  location 
and  conversion  type  to  the  GE  room  results.  Table  7.6  contains  the  derived  tissue 
exposure  results  for  flie  GE  room.  Table  7.7  contains  the  annual  tissue  dose  equivalents 
for  the  GE  room  using  the  usage  factors  described  earlier. 

Again,  the  positive  x  and  negative  y  tally  detector  locations  u-e  pertinent  to  die 
occupational  dose  and  the  negative  x  and  positive  y  tally  detector  locations  for  the  non- 
occupational  dose.  The  results  from  the  GE  room  show  unexpected^  large  doses  at  the 
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Table  7.  5  Ratios  Between  Exposures  Outside  the  CT  Suite  at  WPMC 


WtoW’ 


Ddectcr 

location 

Headto^^tdcnm  Giqurisoi 

]/W»adt  Fhicfiani 
1/lff'i^txlDnai  ihor 

]/S"Ifead:  Haotiond 
]/8";Uxkmn  Ehor 

positive  and  negative  waD-axes  tally  detectors.  They  are  well  above  the  limits  set  in 
10CFR20  for  the  occupational  dose  and,  even  with  the  additional  shielding,  above  the  non- 
occupational  dose  limits.  The  floor,  again,  is  above  the  non-occupational  limit.  It  is, 
however,  below  the  occupational  dose  limit.  The  positive  and  negative  x-room  and  x- 
source  exposures,  particular^  in  the  room  with  1/16”  shielding,  were  much  greater, 
relatively,  than  the  WPMC  CT  suite  exposures.  Recall  that  the  differences  between  the 
rooms  were  size  and  angle  of  rotation  of  the  CT  unit.  As  seen  before,  room  size,  i.e. 
distance  from  isocenter,  can  have  a  significant  impact  on  dose.  It  is  also  possible  that  the 
change  in  rotation  picked  up  a  “hot  spot”  not  sampled  in  the  WPMC  scenario.  Refer,  for 
example,  to  Figure  7.6  showing  the  exposure  variation  at  SO  cm  intervals  aloi^  the  wall. 
The  talfy  detector  located  -100  cm  fi'om  the  negative  x  wall  axis  shows  one  such  hot  spot. 
Similarly,  the  tally  detector  located  +150  cm  fi'om  the  positive  x-room  axis  has  the  same 
abnormally  large  exposme.  These  could  indeed  be  points  of  excess  radiation  or  they  could 
show  a  weakness  in  how  the  CT  unit  (CT  detector,  shields,  etc.)  was  modeled.  That  could 
only  be  dcteimincd  by  sampling  additional  tally  detector  locations  and/or  increasing  the 
fidelity  of  tiie  model. 
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Table  7.  6  Converted  GE  Room  Exposure  (R  s‘^  particle'*) 


G£  Room,  1/16"  Lead  Shiekiiiig,  Exposure 

Detector 

Head 

Fractional 

Abdomm 

Fractional 

Location 

Error 

Error 

Source,  -x 

1.870E>19 

0.032 

1.300E-19 

0.117 

Source,  +x 

1.396^19 

0.034 

1.033E-19 

0.119 

Bod^,+y 

2.785E^20 

0.077 

9.845E-21 

0.139 

Body,-y 

4.349F>21 

0.112 

1.945&-21 

0.227 

Room,  -X 

6.666E-19 

0.022 

4.048E-19 

0.098 

Room,  +x 

9.131F>19 

0.019 

7.373F>19 

0.077 

Ceiling 

7.762E-21 

0.062 

4.065E-21 

0.087 

Floor 

1.178E>19 

0.029 

3.301E-20 

0.057 

Room,  +y 

3.854&-20 

0.079 

2.106E-20 

0.122 

Room,  -y 

7.815E-21 

0.153 

8.634E-21 

0.210 

G£  Room,  1/8"  Lead  Shielding,  Exposure 

Detector 

Head 

Fractional 

Abdomen 

Fractional 

Location 

Eiror 

Error 

Source,  -x 

2.549E-20 

0.149 

2.008E-20 

0.190 

Source,  +x 

2.132B-20 

0.167 

1.723E-20 

0205 

Body,+y 

2.037E-21 

0.151 

1.580E>21 

0.205 

Body,-y 

6.311E-22 

0.243 

6.297&-22 

0.275 

Room,  -X 

3.819E-20 

0.117 

3.001Er20 

0.172 

Room,  +x 

6.085E-20 

0.119 

5.185E-20 

0.136 

Ceiling 

7.294E>-21 

0.114 

5.197E-21 

0.113 

Floor 

4.674E-20 

0.061 

3.348E-20 

0.069 

Room,  +y 

2.391E^21 

0.137 

1.885E-21 

0.179 

Room,  -y 

5.941F>22 

0.276 

6.130E-22 

0.290 
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Table  7.  7  Annual  Tissue  Dose  Equivalent,  GE  Room 


GE  Room,  1/16” 

Lead  Shielding 

Detector 

Annual  Tissue  Dose  Eqivalant,  mSv  | 

Location 

Head 

Ab 

Total  1 

Source,  -x 

18.233 

32.067 

50.299 

Source,  +x‘^ 

13.616 

25.472 

39.088 

Body,  +y 

2.716 

2.428 

5.145 

Body,  -y^ 

0.424 

0.480 

0.904 

Room,  -X 

65.013 

99.867 

164.880 

Room,  +x'^ 

89.051 

181.868 

270.918 

Ceiling"^ 

0.757 

1.003 

1.760 

Floor^ 

11.489 

8.136 

19.363 

Room,  +y 

3.759 

5.195 

8.954 

Room,  -y^ 

0.762 

2.130 

2.892 

GE  Room,  1/8”  Lead  Shielding 

Detector 

Annual  Tissue  Dose  Eqivalant,  mSv 

Location 

Head 

Ab 

Total  1 

Source,  -x® 

2.486 

4.954 

7.439 

Source,  +x 

2.079 

4.250 

6.779 

Body,  +y® 

0.199 

0.390 

0.593 

Body,  -y 

0.062 

0.155 

0.275 

Room,  -X® 

3.724 

7.402 

11.620 

Room,  +x 

5.934 

12.790 

18.761 

Ceiling® 

0.711 

1.282 

1.993 

Floor® 

4.559 

8.260 

12.818 

Room,  +y® 

0.233 

0.465 

0.692 

Room,  -y 

0.058 

0.151 

0.266 

^Detector  locations  applicable  to  occupational  dose. 


^Detector  locations  applicable  to  non-occupational  dose. 


8.  Summary  and  Discussion 

The  purpose  of  this  thesis  was  to  model  a  computerized  tomography  (CT)  unit  and 
suite  using  Monte  Carlo  mediods  to  determine  the  non-occupational  dose  outside  the  suite. 
This  was  driven  by  the  recent  inclusion  of  die  most  recent  NCRP  dose  limit 
recommendations  into  10CFR20,  which  contain  die  national  standards  for  protection 
against  radiation.  The  dose  limits  above  background  are  1  mSv,  non-occupational 
continuous  exposure,  5  mSv,  non-occupational  infrequent  ^posure,  and  50  mSv, 
occupational.  The  occupational  limit  would  be  applied  to  die  medical  technicians  in  the 
suite,  the  non-occiqiational  continuous  exposure  limit  to  allied  medical  personnel  widi,  for 
instance,  offices  adjacent  to  the  suite,  and  the  non-occupational  infrequent  exposure  limits 
to  other  medical  center  patients  traversing  die  corridors  adjacent  to  the  suite.  Prior  to  this 
anafysis,  CT  dose  was  difficult  to  determine  because  die  scattered  radiation  spectrum  from 
a  CT  was  not  available.  MCNP,  however,  bins  its  result  by  particle  energy,  thus,  providing 
the  scattered  spectra. 

These  goals  were  accomplished  using  a  method  outlined  by  Robert  Metzger 
(Metzger,  1993)  using  MCNP,  a  general-purpose  Monte  Carlo  n-particle  transport  model 
developed  by  Los  Alamos  National  Laboratory  (MCNP,  1993).  The  CT,  a  continuous^ 
rotating  source,  was  modeled  as  32  discrete  beam  positions,  each  comprised  of  35  narrow 
conic  beams  to  form  the  characteristic  CT  fan  beam.  The  target  was  the  standard  human 
phantom  composed  entirely  of  a  homogeneous  tissue  defined  in  MERD  Pamphlet  5  (with 
some  minor  changes;  Snyder,  1978).  Separate  runs  were  done  simulating  a  head  scan  and 
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an  abdomen  scan.  Average  usage  of  the  CT  at  WPMC  was  used  to  establish  usage  factors 
and  workload  (Michael,  1995).  The  suite  walls  were  modeled  alternately  with  1/16”  and 
1/8”  lead  shielding  sandwiched  between  1.4  cm  duck  gypsum  dtywall.  The  ceiling  and 
floor  were  modeled  as  six  inch  slabs  of  concrete. 

Validation  of  the  modeling  effort  was  accon:q)lished  by  comparing  the  exposure  of 
film  badges  placed  widiin  the  CT  suite  to  the  results  obtained  from  modeling  the  badge 
locations  with  MCNP.  The  badge/MCNP  results  agreed  by  approximately  a  factor  of  two 
(some  more,  some  less)  with  the  exception  of  badge  locations  that  were  likefy  blocked 
during  the  exposure  period.  In  general,  the  MCNP  results  ctunpared  favorably  with  typical 
scatter  surveys  of  the  G£  HiSpeed  Advantage  CT  scanner  unit.  The  results  obtained  in  the 
analysis  are  presented  below.  They  are  followed  by  recommended  improvements  to  the 
methodology  and  possible  follow-on  research. 

8.1  Analysis  Results 

In  addition  to  the  primary  goal  of  obtaining  the  non-occupational  dose  outside  the 
CT  suite,  the  analysis  for  1/16”  lead  shielding  (typically  used  in  walls  adjacent  to 
occupational  areas)  and  1/8”  lead  shielding  (typically  used  in  walls  adjacent  non- 
occupational  areas)  provided  the  occupational  dose  as  a  by-product.  MCNP,  the  model 
used  for  this  analysis,  can  bin  the  total  number  of  particles  at  a  detector  location  by  particle 
energy.  Thus,  a  scattered  energy  spectrum  was  obtained  at  each  location.  The  scattered 
spectra  from  a  CT  unit,  previously  unavailable  (Simpkin,  1993),  was  also  determined. 
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8. 1. 1  Non-Occupational  Dose. 


The  dose  from  a  CT  unit  at  a  given  location  is  a  function  of  many  variables:  source 
strengdi,  spectrum  and  shape;  target  size,  location  and  composition;  CT  shielding, 
frequency  of  use  and  mode  of  use;  and  wall/floor  shielding  between  the  source  and 
detector  location.  AH  of  these  factors  were  included  in  the  dose  analysis.  10CFR20 
contains  the  maximum  non-occupational  annual  dose  limits  for  continuous  and  infrequent 
exposure.  These  would  correspond  to  allied  medical  persoimel  and  non-related  medical 
center  patients,  respective^.  The  negative  x-room  axis,  negative  x  source  axis ,  positive  y- 
room  axis,  positive  y-bo^  axis,  floor  and  ceiling  results  for  1/8”  lead  shielding 
(corresponding  to  the  shielding  and  layout  at  WPMC)  were  used  to  determine  the  non- 
occupational  dose.  The  analysis  (Table  7.4,  1/8”  shielding)  determined  diat  die  dose  limit 
for  continuous  exposure  (allied  medical  personnel)  of  1  mSv  required  by  10CFR20  was 
exceeded  below  the  floor  and  above  the  ceiling  of  the  CT  suite.  The  dose  limit  of  5  mSv 
for  infrequent  exposure  (non-related  patients)  was  exceed  onfy  below  the  floor.  The 
calculated  dose  values  above  the  ceiling  and  below  the  floor,  however,  may  be  more  a 
function  of  modeling  than  actual  dose  since  the  floor  and  ceiling  concrete  was  modeled  as 
six  inch  slabs  instead  of  the  sbc  inch  slab  plus  two  inches  of  waffled  concrete  on  each  side 
that  are  in  die  CT  suite. 

8.1.2  Occupational  Dose. 

The  positive  x-room  axis,  positive  x  source  axis,  negative  y-room  axis,  negative  y- 
bo^  axis,  floor  and  ceiling  results  for  1/16”  lead  shielding  were  used  to  determine  the 
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occupational  dose.  The  analysis  (Table  7.4,  1/16”  shielding)  determined  that  at  all 
locations,  dose  to  occupational  workers  was  below  the  dose  limit  of  SO  mSv  required  by 
10CFR20. 

8.1.3  Scattered  Spectrum. 

In  general,  die  spectrum  of  scattered  radiation  is  a  sli^tly  softened  source  spectra 
and  is  best  represented  by  the  bone-in  test  results  seen  in  F^.  C.7  and  C.8.  The  spectra 
obtained  inside  the  WPMC  CT  suite  show  a  spectral  contribution  from  the  cadmium 
tungstate  CT  detector  which  may  not  be  relevant  to  all  CT  units,  hi  addition,  the  phantom 
used  in  the  analyses  inside  the  suite  did  not  include  the  skeleton,  which  was  shown  to  affect 
the  spectra. 


8.1.4  Application  of  WPMC  Results  to  a  Generic  CT  Suite. 

The  talfy  behavior  and  spectral  results  for  the  WPMC  CT  suite  with  1/16”  lead 
shielding  for  a  head  scan  were  very  similar  to  the  analogous  G£  minimum  recommended 
room  results.  The  agreement  suggested  that  the  differences  noted  for  the  WPMC  suite 
from  1/16”  to  1/8”  lead  shielding,  and,  from  a  head  to  an  abdomen  scan  could  be  applied 
to  the  single  run  done  for  the  GE  room  with  1/16”  lead  shielding  for  the  head  scan  via 
conversion  factors.  The  conversion  factors  were  determined  by  comparing  1/16”  to  1/8” 
lead  shielded  room  results  at  each  detector  location  for  the  head  and  abdomen  scans  and  by 
comparing  the  head  to  abdomen  scan  results  for  the  1/16”  and  1/8”  shielded  rooms.  The 
l/16”-to-l/8”  conversion  factors  (for  each  detector  location)  were  applied  to  the  results  for 
GE  1/16”  shielded  room  (head  scan)  to  obtain  GE  1/8”  shielded  room  (head  scan)  values. 
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Similarly,  Ihe  head-to-abdomen  conversion  factors  were  applied  to  tiie  1/16”  and  1/8” 
shielded  room  results  for  the  head  scan  to  obtain  results  for  an  abdomen  scan  in  those 
rooms. 

In  general,  comparing  the  total  annual  dose  equivalents  by  location  for  tiie  WPMC  suite 
in  Tables  7.4  and  for  tfie  GE  room  in  Table  7.8,  similar  results  are  seen.  The  dose 
equivalents  outside  the  GE  room  are  slightfy  higher  than  outside  the  suite  at  WPMC.  This 
was  expected  because  the  GE  room  is  smaller  than  tiie  WPMC  room.  The  detectors 
located,  however,  on  the  positive  and  negative  x  walls  outside  the  GE  room  were 
unexpected^  high.  Thorough  review  of  the  input  files  revealed  no  cause  for  this,  i.e.  no 
mistakes  were  found  in  the  files.  Recalling  fiiat  the  CT  unit  is  rotated  at  slightly  diJSerent 
angles  for  the  two  rooms,  a  possible  explanation  is  that  the  GE  room  x-wall  detectors  were 
placed  directfy  on  or  much  closer  to  a  “hot  q)ot”  than  the  detectors  outside  the  WPMC 
room. 

With  die  notable  exception  of  the  positive  x-room  axis  detector  location,  the  same 
conclusions  concerning  conformity  to  the  10CFR20  dose  limits  are  reached.  That  is,  the 
annual  dose  limit  for  continuous  exposure  (allied  medical  personnel)  of  1  mSv  required  by 
10CFR20  was  exceeded  below  the  floor,  above  the  ceiling,  and  outside  the  negative  x  wall 
of  the  GE  room.  The  annual  dose  limit  of  5  mSv  for  infi'equent  exposure  (non-related 
patients)  was  exceeded  only  below  the  floor  and  outside  the  negative-x  wall.  At  all  the 
locations  but  the  positive  x-room  axis,  dose  to  occupational  workers  was  below  the  annual 
dose  limit  of  50  mSv  required  by  10CFR20. 
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8.2  Recommended  Improvements 


As  with  any  analysis,  there  is  room  for  improvement  in  this  thesis.  The  following 
sections  describe  problem  areas  and  suggested  fixes. 

8.2.1  Geometry. 

Two  needed  improvements  to  problem  geometry  were  made  obvious  in  the 
analysis:  (1)  better  definition  of  the  concrete  floor  and  ceiling,  and  (2)  defining  a  phantom 
that  fits  within  the  CT  fan  beam.  The  first  could  be  investigated  via  test  runs  using  much 
simpler  scatter  radiation  sources  than  die  CT  unit  with  phantom.  The  test  would  determine 
primarily  if  the  addition  of  die  waffled  concrete  impacted  the  external  dose.  The  second 
problem  could  be  corrected  either  by  wideniig  the  fan  beam  (via  additional  conic  beams 
and  )  or,  more  simpfy,  by  reducing  the  dimensions  of  the  phantom’s  abdomen. 

8.2.2  Source. 

The  source  spectra  described  in  the  input  files  was  unwiel^  since  every  pencil 
beam  in  the  fan  had  a  different  spectrum  resulting  fi'om  passing  beam  through  the  Teflon 
bowtie  filter.  Examination  of  the  different  spectra  (not  included  in  diis  thesis  because  they 
were  proprietary  to  GE)  revealed  that  their  magnitude  across  the  spectrum  differed  fi'om 
one  another  by  a  neaify  constant  factor.  Streamlining  the  source  definition  could  be  done 
by  replacing  the  different  spectra  with  probabilities  determined  by  a  normalization  of  the 
integrated  individual  spectra.  For  example,  ifthe  spectra  at  IT.S^’ had  one  half  the 

integrated  source  strength  of  the  central  spectra,  the  beam  at  17.  (with  the  same  spectra 
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as  the  central  beam)  would  be  assigned  a  relative  probability  of  one  half  that  of  the  central 
beam.  Although  it  would  not  reduce  run  time  per  number  of  particles,  it  would  certainly 
reduce  input  file  size. 

8.2.3  AirKerma. 

Since  the  air  kerma  value  for  the  particular  tube  used  in  this  analysis  was  not 
available,  die  air  kerma  value  used  was  obtained  from  a  chart  of  average  air  kerma  for 
diagnostic  X-ray  tubes  of  varying  potentials  (NCRP,  1989).  This  average  air  kerma, 
however,  did  not  account  for  the  changing  spectrum  as  a  function  of  fan  beam  angle  as 
described  above.  An  improvement  to  the  air  kerma  figure  would  require  file  source 
strength  for  each  energy  bin  in  particles  per  some  time  interval,  which,  when  given  file 
values  for  energy  deposited  in  air  as  a  function  of  particle  energy  (the  definition  of  kerma, 
section  2.8.4),  could  be  used  to  determine  the  actual  air  kerma  of  fliis  particular  X-ray 
tube. 


8.2.4  Sampling  Interval. 

Recall  that  some  of  the  talty  plots  displayed  large  peaks.  The  beam  positions  around  these 
peaks  show  a  need  for  sampling  additional  beam  positions  to  better  define  their  shape  and 
thus  the  average.  For  example.  Fig.  G.l,  the  positive  source  axis  detector  talty  from  file 
WPMC  room  with  1/16”  shielding,  shows  a  twin-peaked  behavior  on  either  side  of  the  Pt  0 
beam  position.  Clcarty,  these  peaks  contribute  heavity  to  the  mean  exposure.  If  the  peaks 
were,  in  fact,  wider,  the  talty  was  underestimated;  if  the  peaks  are  very  narrow,  the  talty 
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was  overestimated.  All  detector  locations  would  not  need  the  additional  sampling  points,  it 


could  be  done  on  as  as-needed  basis. 


8.2.5  Computer  Resources. 

This  analysis,  despite  all  efforts  to  conserve  computer  resources  by  reducing  file 
size,  run  time  and  number  of  points  sampled,  still  required  imposing  computer  resources. 
Each  individual  beam  position  run  was  done  on  a  separate  Sun  Sparc20  (AFTT  computer 
lab,  accessible  to  all  AFTT  students  and  faculty).  Run  times  of  one  to  two  weeks  were 
logged  for  two  million  particles  (sharing  resources  with  other  users)  for  ten  detector 
locations  plus  the  source  detector.  Elimination  of  detector  locations  resulted  in  linear 
reductions  in  run  times.  Each  of  the  32  beam  positions  produced  a  set  of  computer  files 
containing  the  input  file,  ouq>ut  file,  continuation  file  (RUNTPE)  and  tally  file  (MCTAL) 
with  a  total  approximately  three  megabytes.  The  ei^t  runs  (CT  in  X-ray  room;  CT  inside 
WPMC,  head  and  abdomen;  CT  at  WPMC  with  1/16”  and  1/8”  shielding  for  head  and 
abdomen;  and  CT  in  GE  room,  head  and  abdomen)  dius  used  almost  800  megabytes  of 
disk  space.  It  is  safe  to  say  that  the  average  radiation  safety  officer  would  not  have 
computer  resources  of  this  magnitude  readily  available. 

A  possible  answer  to  this  problem  requires  viewing  the  problem  a  different  way.  ffi, 
instead  of  performing  this  anafysis  for  a  given  room,  it  focused  instead  on  the  CT  unit  and 
phantom,  a  data  set  could  be  produced  that  would  have  greater  usefulness.  Imagine  a 
spherical  surface  surrounding  the  CT  unit  and  phantom.  If  the  dose  were  determined  on 
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this  surface,  it  could  be  extrapolated  out  to  and  through  whatever  surfaces  conJSned  it. 
This,  in  effect,  would  be  equivalent  to  defining  a  non-isotropically  radiating  point  source 
representative  of  the  scattered  radiation.  The  initial  effort  would  require  comparable 
computer  resources  described  above,  but  they  would  not  need  to  be  repeated  for  every 
diagnostic  facility. 

8.3  Follow-On  Research. 

At  die  conclusion  of  any  research  defining  or  using  a  new  method  of  analysis, 
foUoW’On  studies  immediate^  come  to  mind.  This  thesis  is  no  exception.  The  results 
raised  new  questions  about  methods,  materials  or  improved  sources.  In  particular,  two 
areas  of  further  investigation  present  themselves:  complete  definition  of  die  scattered 
spectrum  and  dose  at  all  points  surrounding  die  source  and  target. 

8.3.1  Scattered  Spectrum. 

Earfy  in  thesis  research,  D.  Simpkin  (Simpkin,  1995)  questioned  die  ability  to 
complete  this  thesis  because  of  the  lack  of  a  scattered  spectrum  fi*om  a  CT  unit.  Previous 
methods  of  determining  dose  would  have  required  it.  The  use  of  MCNP,  however, 
precluded  this.  MCNP  could  produce  a  scattered  spectrum  given  a  well-defined  source 
and  problem  geometry.  Thus,  one  product  of  this  research  was  the  elusive  scattered 
spectrum.  It,  however,  was  only  defined  at  the  particular  detector  locations.  To  be  truly 
useful,  the  spectrum  should  be  defined  at  all  points  surrounding  die  source  and  target  for  a 
variety  of  taigets,  sources,  geometries  and  compositions.  MCNP  has  a  surface  current 
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capability  that  would  provide  this  type  of  information.  This  is  an  area  ripe  for  further 
investigation. 

8.3.2  Angle  and  Position  Dependent  Dose. 

Another  area  for  further  investigation  was  presented  in  section  8.2.4,  that  is, 
determining  die  dose  on  a  spherical  surface  surrounding  the  CT  unit  and  target.  Similar 
techniques  used  to  define  the  spectrum  at  all  points  could  be  used.  In  fact,  the  use  of 
MCNP  would  natural^  tie  diese  two  together  since  the  reported  e?q)osure  talfy  is  a  function 
of  the  calculated  spectrum.  As  discussed  above,  such  a  product  would  be  very  useful  to 
radiation  protection  professionals  to  determine  the  most  efficient  and  effective  shielding 
needed  aroimd  a  CT  suite. 

8.4  Conclusions 

The  primaiy  goal  of  determining  the  equivalent  dose  outside  a  computerized 
tomography  (CT)  suite  was  accomplished.  The  dose  was  compared  to  the  recently  revised 
dose  limits  required  by  10CFR20  for  occupational  and  non-occupational  areas.  At  most 
test  locations,  the  non-occupational  and  occupational  dose  limits  were  not  exceeded.  The 
limits  are  exceeded  below  floor  and  above  ceiling  of  the  WPMC  CT  suite  and  the  GE 
room,  subject,  however,  to  further  anafysis.  The  GE  room  also  exceeded  the  non- 
occupational  limits  outside  the  negative-x  wall.  In  addition,  a  scattered  radiation  spectrum 
firom  a  CT  unit,  previously  unavailable,  was  produced.  The  method  of  using  a  Monte 
Carlo  model  described  by  Met2ger,  et.  al.  (Metzger,  1993)  is  considered  valid  to  determine 
dose  outside  diagnostic  X-r^  facilities. 
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Appendix  A  Sample  Input  File 

Section  I;  Cells 

WP  CT  room  w/MIRD  phantom  on 
c  table  in  center  of  room, 
c  Uses  geometiy  splitting;  floor,  ceiling 
c  and  lead  lining  are  divided 
c 


c 

Lungs: 

101 

1  -.2958  7  -12  13 

imp5)=4 

103 

1  -.2958  7  -14  15 

imp:p=4 

c 

Torso 

105 

2  -.9869  8  -6  -10  #101  #103 

imp:p=l 

c 

Legs 

106 

2 -.9869  -8-16 

imp5)=l 

c 

Face 

7 

2  -.9869  -5  6  -1 

imp:p=l 

c 

Top  of  head: 

8 

2  -.9869  5-2 

imp:p=l 

c 

hiterior  of  room 

9 

3  -.0012  (-8:6:10)  #8  #7  #30  #106  #40  #50 

401-411  -201  211301-311 

imp:p=l 

c 

Inner  gypsum  layer: 

111 

4  -.73  (201:-301:-211:311)  -202  302  212 

-312 

-411  401 

imp:p=4 

lines  skipped 


c  Sphere  surrounding  room 

19  3  -.0012  (207:-307:-217:317:416:-406)  -43 

20  0  43 

c  Detector 

30  10  -7.00  60  -61  71  -70  72  -73  74 

c  Psuedo  Gantry 
40  11  -2.69  -81  80  83  -82  (-71:70) 

c  Shield 

50  11-2.69  -9190-92  73 


imp:p=4096 

imp:p=0 

imp:p=4 

imp:p=4 

imp.‘p=4 
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Section  II:  Surfaces 


c  MERD  Phantom  head  in  center  of  beam: 

c 

1  sq  100  0  49  0  0  0  -4900  0  1  0  $face 

2  sq  7225  4900  3540.25  0  0  0  -354025  0  1  0  Stop  o'  head 

5  py  1  Stop  of  face  and  side  of  detector 

6  py  -14.5  Stop  of  torso 

7  py  -41  Sbottom  of  lungs 

8  py-84.5  Stop  of  legs 

10  sq  100  0  400  0  0  0 -40000  0  1  0  Stonso 

12  sq  32400  1406.25  14400  0  0  0  -810000  8.5  -41  0  S  lungs 

13  sq  32400  1406.25  14400  0  0  0  -810000  2.5  -41  0 

14  sq  32400  1406.25  14400  0  0  0  -810000  -8.5  -41  0 

15  sq  32400  1406.25  14400  0  0  0  -810000  -2.5  -41  0 

16  sq  16  1  64  0  0  0  -6400  0  -84.5  0  Slegs 
c  50'  Sphere  of  air  surrounding  room: 

43  so  1500 
c  Wall  to  -X  of  patient 

201  p -1.901  1  0  526.599 

202  p  -1.901  1  0  529.6068 

lines  skipped 


406  pz  -123.6 
c  Ceiling 

411  pz  287.02 

412  pz  290.02 

413  pz  293.02 

414  pz  296.02 

415  pz  299.02 

416  pz  302.02 
c  detector 

60  c/y-63  0110 

61  c/y-63  0115 

70  p -0.50953  01  32.1001 

71  p  0.50953  01 -32.1001 

72  py-2 

73  py2 

74  px  0 

c  Aluminum  block  to  represent  power  supply,  support,  etc  of  gantiy 

80  cy  65 

81  cy  68 

82  py  15 

83  py  -15 
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Section  III:  Materials 


mode  p 

c  Defmition  for  lung  tissue: 

ml  1000  -.1021 
6000  -.1001 
7000  -.0280 
8000  -.7596 
11000  -.0019 
15000  -.00081 
16000  -.0023 
17000  -.0027 
19000  -.0020 
26000  -.00037 

c  Standard  tissue  (non-lung): 
m2  1000  -.1047 
6000  -.2302 
7000  -.0234 
8000  -.6321 
11000 -.0013 
12000  -.00015 
15000  -.0024 
16000  -.0022 
17000  -.0014 
19000  -.0021 


lines  skipped 


c  Cadmium  Tungstite,  density=~7.00  g/cc 
mil  74000  1  48000  1  8000  4 


Section  IV:  Source 


c  CT  source  with  head  in  center 
c  spectrum  from  G£ 
c  120  kVp  with  bowtie  filter 
sdef  pos=dl  dir=d2  vec=Q)os=d3  eig=Q)os=d4 
sil  1  -63  0  0 

-63  0  0 

-63  0  0 


lines  skipped 


-63  0  0 
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-63  0  0 
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Section  V:  Termination 
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Appendix  C:  Bone-In  vs  No-Bone  Spectra 


Figure  C.l  Normalized  spectra  without  Roentgen  conversion  for  positive  x  tally  detector  for  bone-in 
and  no-bone  comparison. 


Figure  G.2  Normalized  spectra  without  Roentgen  conversion  for  negative  x  tally  detector  for  bone-in 
and  no-bone  comparison. 
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Appendix  C;  Bone-In  vs  No-Bone  Spectra 


Figure  C.4  Normalized  spectra  without  Roentgen  conversion  for  negative  y  tally  detector  for  bone-in 
and  no-bone  comparison. 


C-2 


Appendix  C:  Bone-ln  vs  No-Bone  Spectra 


Figure  C.6  Normalized  spectra  without  Roentgen  conversion  for  floor  tally  detector  for  bone-in  and  no¬ 


bone  comparison. 
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Appendix  C;  Bone-In  ys  No-Bone  Spectra 


Figure  C.7  Normalized  spectra  without  Roentgen  conversion  for  positive  gantry  tally  detector  for  bone- 
in  and  no-bone  comparison. 


Figure  C.8  Normalized  spectra  without  Roentgen  conversion  for  negative  gantry  tally  detector  for  bone- 
in  and  no-bone  comparison. 
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Figure  D.l  Positive  x  exposure  (R/s/particle),  CT  in  X-ray  room.  0  degree  direct 
beam  position. 


Figure  D.2  Negative  x  exposure  (R/s/particIe)  for  CT  In  X-ray  room.  180  degree 
direct  beam  position. 
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Appendix  D:  CT  in  X-Ray  Room 
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Figure  D.3  Ceilii^  exposure  (R/s/particle)  for  CT  in  X-ray  room.  90  degree  direct 
beam  position. 
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Figure  D.4  Floor  exposure(R/s/particle)  for  CT  in  X-ray  room.  270  degree  direct 
beam  position. 
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Appendix  D:  CT  in  X-Ray  Room 


Figure  D.5  Positive  y  exposure(R/s/particIe)  for  CT  in  X-ray  room.  No  direct  beam 
position. 


Figure  D.6  Negative  y  exposure  (R/s/particle)  for  CT  in  X-ray  room.  No  degree 
direct  beam  position. 
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Appendix  D:  CT  in  X-Ray  Room 


Figure  D.7  Source  exposure  (R/s/particle)  for  CT  in  X-ray  room.  This  is 
representative  of  source  detector  tallies  for  all  analyses. 
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Appendix  D:  CT  in  X-ray  Room 


Figure  D.8  Positive  and  negative  x  detector  normalized  average  spectra  without 
Roentgen  conversion  for  CT  in  X-ray  room. 


Figure  D.9  Positive  and  negative  y  detector  normalized  average  spectra  without 
Roentgen  conversion  for  CT  in  X-ray  room. 
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Appendix  D:  CT  in  X-ray  Room 


Figure  D.IO  Ceiling  and  floor  detector  normalized  average  spectra  without 
Roentgen  conversion  for  CT  in  X-ray  room. 


Appendix  E:  CT  Inside  WPMC,  Head  Scan 


Figure  E.1  Exposure  (R/s/particIe)  for  positive  x  source  axis  detector  inside  CT  suite  at 
WPMC,  bead  scan.  0  degree  direct  beam. 


Figure  E.2.  Exposure  (R/s/particle)  for  negative  x  source  axis  detector  inside  CT  suite  at 
WPMC,  bead  scan,  direct  beam  at  180  degrees. 
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Appendix  E:  CT  Inside  WPMC,  Head  Scan 


Figure  E.3  Exposure  (R/s/particle)  for  positive  x  room  axis  detector  inside  CT  suite  at 
WPMC,  bead  scan,  no  direct  beam. 
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Figure  E.4  Exposure  (R/s/particie)  for  negative  x  room  iixis  detector  inside  CT  suite  at 
WPMC,  bead  scan,  no  direct  beam. 


Appendix  E:  CT  Inside  WPMC,  Head  Scan 


2.5E-17 

2.3E-17 

2.1E-17 

1.9E-17 

1.7E-17 

1.5E-17 


0  45  90  135  180  225  270  315  360 


Source  Position,  Deg 

Figure  E.5  Exposure  (R/s/particle)  for  positive  j  body  axis  detector  inside  CT  suite  at 
WPMC,  head  scan,  no  direct  beam. 


Figure  E.6  Exposure  (R/s/particle)  for  negative  y  body  itxis  detector  inside  CT  suite  at 
WPMC,  head  scan,  no  direct  beam. 
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Appendix  E:  CT  Inside  WPMC,  Head  Scan 


Figure  E.7  Exposure  (R/s/particle)  for  positive  j  room  axis  detector  inside  CT  suite  at 


WPMC,  head  scan,  no  direct  beam. 
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Appendix  E:  CT  Inside  WPMC,  Head  Scan 


Figure  £.10.  Exposure  (R/s/particle)  for  floor  detector  inside  CT  suite  at  WPMC,  bead  scan. 
270  degree  direct  beam. 
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Appendix  E:  Inside  WPMC  CT  Suite,  Head  Scan 


Figure  £.11  Average  normalized  spectrum  and  spectra  for  beam  positions  0, 11  and  22  for 
positive  X  source  axis  detectors  inside  WPMC  CT  suite  for  head  slice 


Figure  £.12  Average  normalized  spectrum  and  spectra  for  beam  positions  157, 168  and 
180  for  positive  x  source  axis  detectors  inside  WPMC  CT  suite  for  head  slice 


E-6 


Appendix  E:  Inside  WPMC  CT  Suite,  Head  Scan 


inside  WPMC  CT  suite  for  head  slice. 


Figure  E.14  Average  normalized  spectrum  for  positive  and  negative  y  body  axis  detectors 
inside  WPMC  CT  suite  for  head  slice. 
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Appendix  E:  Inside  WPMC  CT  Suite,  Head  Scan 


Appendix  E:  Inside  WPMC  CT  Suite,  Head  Scan 


Figure  E.17  Average  normalized  spectrum  and  individual  spectra  for  beam  positions  247 
through  292  for  floor  detector  inside  WPMC  CT  suite  for  head  slice. 
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Appendix  F:  Inside  WPMC  CT  Suite,  Abdomen  Scan 


Source  Position  (deg) 

Figure  F.l  Exposure  (R/s/particle)  for  positive  x  source  axis  detector  inside  CT  suite  at  WPMC, 
abdomen  scan. 


Figure  F.2  Exposure  (R/s/particle)  for  negative  x  source  iuxis  detector  inside  CT  suite  at  WPMC, 
abdomen  scan. 
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Appendix  F:  Inside  WPMC  CT  Suite^  Abdomen  Scan 


Figure  F.3  Exposure  (R/s/particle)  for  positive  x  room  axis  detector  inside  CT  suite  at  WPMC,  abdomen 
scan. 


scan. 
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Appendix  F:  Inside  WPMC  CT  Suite,  Abdomen  Scan 


Figure  F.5  Exposure  (R/s/particle)  for  positive  y  body  axis  detector  inside  CT  suite  at  WPMC,  abdomen 
scan. 


Figure  F.6  Exposure  (R/s/particle)  for  negative  y  body  axis  detector  inside  CT  suite  at  WPMC,  abdomen 
scan. 
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Appendix  F:  Inside  WPMC  CT  Suite,  Abdomen  Scan 


Figure  F.7  Exposure  (R/s/purticle)  for  positive  y  room  axis  detector  inside  CT  suite  at  WPMC,  abdomen 
scan. 


Figure  F.8  Exposure  (R/s/particle)  position  for  negative  y  room  axis  detector  inside  CT  suite  at  WPMC, 
abdomen  scan. 
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Appendix  F:  Inside  WPMC  CT  Suite,  Abdomen  Scan 
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Figure  F.9  Exposure  (R/s/particle)  for  ceiling  detector  inside  CT  suite  at  WPMC,  abdomen  scan. 


Figure  F.IO  Elxposure  (R/s/particle)  for  floor  detector  inside  CT  suite  at  WPMC,  abdomen  scan. 
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Source,  +x 
Source,  -x 


Figure  F.ll  Normalized  average  spectrum  for  positive  and  negative  x  source  axis  detectors  inside  CT 
suite  at  WPMC)  abdomen  scan. 


Figure  F.ll  Normalized  average  spectrum  for  positive  and  negative  x  room  axis  detectors  inside  CT 
suite  at  WPMC,  abdomen  scan. 
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Appendix  F:  Inside  WPMC  CT  Suite,  Abdomen  Scan 


Figure  F,13  Normalized  average  spectrum  for  positive  and  negative  y  body  axis  detectors  inside  CT  suite 


at  WPMC,  abdomen  scan. 


Appendix  F:  Inside  WPMC  CT  Suite,  Abdomen  Scan 


Figure  F.15  Normalized  average  spectrum  for  ceiling  and  floor  detectors  inside  CT  suite  at  WPMC, 
abdomen  scan. 

oc 
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Appendix  G:  WPMC  CT  Suite.  1/16”  Lead  Shielding,  Head  Scan 
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Figure  G.1  Exposure  (R/s/parUde)  at  each  source  point  for  positive  x  source  axis  detector  with  CT  in 
WPMC,  1/16"  lead,  head  scan. 


Figure  G.2  Exposure  (R/s/parttde)  at  each  source  point  for  negative  x  source  axis  detector  with  CT  in 
WPMC  suite,  1/16"  lead  shielding,  head  scan. 
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Appendix  G:  WPMC  CT  Suite.  1/16”  Lead  Shielding,  Head  Scan 


Figure  GJ  Exposure  (R/s/particle)  at  each  source  point  for  -He  room  axis  detector  with  CT  in  WPMC 
cuite,  1/16"  lead»  head  scan. 


WPMC  suite,  1/16"  lead  shielding,  head  scan. 
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Appendix  G:  WPMC  CT  Suite.  1/16”  Lead  Shielding,  Head  Scan 


Figure  G.5  Exposure  (R/s/parttcle)  at  each  source  point  for  positive  y  body  axis  detector  with  CT  in 
WPMC  room,  1/16”  lead,  head  scan. 


Figure  G.6  Exposure  (R/s/particle)  at  each  source  point  for  negative  y  body  axis  detector  with  CT  in 
WPMC  suite,  1/16”  lead  shielding,  head  scan. 
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Appendix  G:  WPMC  CT  Suite.  1/16”  Lead  Shielding,  Head  Scan 


Figure  G.7  Exposure  (R/s/particle)  at  each  source  point  for  positive  y  room  axis  detector  with  CT  in 
WPMC  suite,  1/16**  lead,  head  scan. 


Figure  G.8  Exposure  (R/s/particle)  at  each  source  point  for  negative  y  room  axis  detector  with  CT  in 
WPMC  suite,  1/16**  lead  shielding,  head  scan. 
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Appendii  G:  WPMC  CT  Suite,  1/16”  Lead  Shielding,  Head  Scan 
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Figure  G.9  Exposure  (R/s/particle)  at  each  source  point  for  ceiling  detector  with  CT  in  WPMC.  1/16” 
lead,  head  scan. 
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Figure  G.11  Normalized,  weighted-average  spectra  for  positive  and  negative  x  source  axes  detectors 
with  CT  in  WPMC,  1/16"  lead  shielding,  head  scan. 


Figure  G.12  Normalized,  weighted-average  spectra  for  positive  and  negative  x  room  axes  detectors  with 
CT  in  WPMC,  1/16"  lead  shielding,  head  scan. 


Appendix  G:  WPMC  CT  Suite.  1/16"  Lead  Shielding,  Head  Scan 


Figure  G.13  Normalued,  weighted-average  spectra  for  positive  and  negative  y  body  axes  detectors  with 
CT  in  WPMC,  1/16"  lead  shielding,  head  scan. 


Figure  G.14  Normalized,  weighted-average  spectra  for  positive  and  negative  y  room  axes  detectors  with 
CT  in  WPMC,  1/16"  lead  shielding,  head  scan. 
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Appendix  G:  WPMC  CT  Suite.  1/16"  Lead  Shielding,  Head  Scan 


1/16"  lead  shielding,  head  scan. 
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Appendix  H:  WPMC  CT  Suite,  1/16”  Lead  Shielding,  Abdomen  Scan 


Figure  H.1  Exposure  (R/s/particle)  at  each  source  point  for  positive  x  source  axis  detector  with  CT  in 
WPMC»  1/16"  lead,  abdomen  scan. 


Figure  H.2  Exposure  (R/s/particle)  at  each  source  point  for  negative  x  source  axis  detector  with  CT  in 
WPMC  suite,  1/16"  lead  shielding,  abdomen  scan. 


H-1 


Appendii  H;  WPMC  CT  Suite,  1/16"  Lead  Shielding,  Abdomen  Scan 
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Figure  H.4  Exposure  (R/s/particle)  at  each  source  point  for  negative  x  room  axis  detector  with  CT  in 
WPMC  suite,  1/16"  lead  shielding,  abdomen  scan. 
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Appendix  H:  WPMC  CT  Suite,  1/16**  Lead  Shielding,  Abdomen  Scan 


Figure  H.6  Exposure  (R/s/purticle)  at  each  source  point  for  negative  y  body  axis  detector  with  CT  in 
WPMC  suite,  1/16''  lead  shielding,  abdomen  scan. 
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Appendis  H;  WPMC  CT  Suite,  1/16"  Lead  Shielding,  Abdomen  Scan 


Figure  H.7  Exposure  (R/s/partIcle)  at  each  source  point  for  positive  y  room  axis  detector  with  CT  In 
WPMC  suite,  1/16"  lead,  abdomen  scan. 


Figure  H.8  Exposure  (R/s/particle)  at  each  source  point  for  negative  y  room  axis  detector  with  CT  in 
WPMC  suite,  1/16"  lead  shielding,  abdomen  scan. 


Appendix  H:  WPMC  CT  Suite,  1/16”  Lead  Shielding,  Abdomen  Scan 


Figure  H.8  Exposure  (R/s/particle)  at  each  source  point  for  negative  y  room  axis  detector  vrith  CT  in 


Figure  H.9  Exposure  (R/s/particle)  at  each  source  point  for  ceiling  detector  with  CT  in  WPMC,  1/16" 
lead,  abdomen  scan. 


lead  shielding,  abdomen  scan. 
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Appendix  H:  WPMC  CT  Suite,  1/16”  Lead  Shielding,  Abdomen  Scan 


Source,  +x 
Source,  -x 


Figure  H.11  Normalized,  weighted-average  spectra  for  positive  and  negative  x  source  axes  detectors 
with  CT  in  WPMC,  1/16"  lead  shielding,  abdomen  scan. 


Figure  H.12  Normalized,  weighted-average  spectra  for  positive  and  negative  x  room  axes  detectors  with 
CT  in  WPMC,  1/16"  lead  shielding,  abdomen  scan. 


H.6 


Figure  H.14  Normaliied)  weighted-average  spectra  for  positive  and  negative  y  room  axes  detectors  with 
CT  in  WPMC,  1/16”  lead  shielding,  abdomen  scan. 
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Appendix  H:  WPMC  CT  Suite,  1/16’*  Lead  Shielding,  Abdomen  Scan 


Figure  H.15  Normalized,  weighted-average  spectra  for  ceiling  and  floor  detectors  with  Ct  in  WPMC, 
1/16"  lead  shielding,  abdomen  scan. 
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Appendix  I:  WPMC  CT  Suite,  1/8”  Lead,  Head  Scan 


Figure  Ll.  Exposure  (R/s/partlde)  at  each  c4MU«e  point  for  the  positive  z  source 
axis  detector  fw  CT  in  WPMC  with  1/8"  lead  shldding,  head  scan. 


Figure  1.2  Exposure  (R/s/particle)  at  each  source  point  for  negative  x  source  axis 
inWPMC  with  1/8”  lead  shfadding,  head  scan. 


Appendix  I:  WPMC  CT  Suite,  1/8"  Lead,  Head  Scan 


Figure  U  Exposure  (R/s/paiUde)  at  each  source  point  for  positive  x  room  axis 
detector  for  CT  detectin' for  CT  in  WPMC  wlthl/S”  lead  shidding,  head  scan. 


Figure  1.4  Exposure  (R/s/particle)  at  each  source  pointfor  negative  x  room  axis 
detector  for  CT  in  WPMC  with  1/8”  lead  shidding,  head  scan. 
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Appendix  I:  WPMC  CT  Suite,  1/8”  Lead,  Head  Scan 


Figure  Exposure  (R/s/paiUde)  at  each  source  point  for  potittvey  body  axis 

detector  for  CT  iuMTMC  with  1/8"  lead  shielding,  head  scan. 
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Figure  1.6  Exposure  (R/s/partide)  at  each  source  point  for  negative  y  body  axis 
CT  In  WPMC  with  1/8"  lead  shldding,  head  scan. 
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Appendix  I:  WPMC  CT  Suite,  1/8”  Lead,  Head  Scan 


Figure  1.7  Exposure  (R/s/paiUde)  at  each  source  point  for  positive  yromn  axis 
detectm*  for  drtectorfwCT  in  WPMC  wHh  1/8”  lead  sUdding,  head  scan. 


Figure  1.8  Exposure  (R/s/partkle)  at  each  source  point  for  negative  x  room  axis 
detector  for  CT  inWPMCivMi  1/8"  lead  shididing,  head  scan. 
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Appendix  I:  WPMC  CT  Suite,  1/8”  Lead,  Head  Scan 


Figure  L9  Exposure  (R/s/partkie)  at  eadi  source  point  for  ceiling  detector  for  CT 
inWPMC  wHli  1/8"  lead  shielding,  head  scan. 


Figure  1.10  Exposure  (R/s/particle)  at  each  source  point  for  floor  detector  for  CT 
inlVPMC  with  1/8”  lead  shielding,  head  scan. 
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Appendix  I:  WPMC  CT  Suite,  1/8"  Lead  Shielding,  Head  Scan 


Figure  Lll  Spectra  for  weighted  average  at  the  positive  and  negative  x 
source  axes  for  the  CT  at  liVPMC  with  1/8"  lead  shielding,  head  scan. 


Figure  L12  Spectra  for  weighted  average  at  the  positive  and  negative  x 
room  axes  for  the  CT  at  WPMC  with  1/8"  lead  shielding,  head  scan. 
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Appendix  I:  WPMC  CT  Suite,  1/8”  Lead  Shielding,  Head  Scan 
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Figure  L13  Spectra  for  weighted  average  at  the  positive  and  negative  y 
body  axes  for  the  CT  at  WPMC  with  1/8”  lead  shielding,  head  scan. 
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Figure  L14.  Spectra  for  weighted  average  at  the  positive  and  negative  y 
room  axes  for  the  CT  at  WPMC  with  1/8”  lead  shielding,  head  scan. 


Appendix  1:  WPMC  CT  Suite,  1/8"  Lead  Shielding,  Head  Scan 


Figure  L15  Spectra  for  weighted  average  at  the  ceiling  and  floor  for  the 
CT  at  WPMC  with  1/8"  lead  shielding,  head  scan. 


Appendix  J:  WPMC  CT  Suite,  1^*'  Lead  Shielding,  Abdomen  Scan 


FIgureJ.l  Exposure  (R/s/paiUcle)  at  each  source  point  fw+x  source  axis  detector 
for  CT  in 'WPMC  with  1/8"  lead  shidding,  abdomen  scan. 


Figure  J.2  Exposure  (R/s/particle)  at  each  source  pointfor  -x  source  axis  detector 
for  CT  in  WPMC  with  1/8"  lead  shielding,  abdomen  scan. 
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Appendix  J;  WPMC  CT  Suite,  1/8”  Lead  Shielding,  Abdomen  Scan 


FlgureJJ  Exposure  (R/t/parttde)  at  each  source  point  for +x  room  axis  detector 
forCT  inWPMC  with  1/8”  iead  shMding,  abdomm  scan. 


FigureJ.4  Eiqiosure(R/s/particle)  at  each  source  point  for -x  room  axis  detector 
forCT  in  WPMC  with  1/8"  lead  shielding,  abdomen  scan. 


J-2 


Appendix  J:  WPMC  CT  Suite,  1/8"  Lead  Shielding,  Abdomen  Scan 


Figure  Exposure  (R/s/partlde)  at  each  source  point  for  -fy  source  axis  detector 
fw  CT  In  WPMC  with  1/8”  lead  shielding,  abdmnoi  scan. 


FlgureJ.6  Exposure  (R/s/partkle)  at  eadi  source  p<dnt  for -;bo^  axis  detector 


for  CT  In  WPMC  ifplth  1/8”  lead  shielding,  abifamien  scan. 
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Appendix  J:  WPMC  CT  Suite,  1/8"  Lead  Shielding,  Abdomen  Scan 
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Figure  J.7  Exposure  (R/s/parOde)  at  each  source  pdnt  for  +y  room  axis  detector 
for  CT  in  WPMC  with  1/8”  lead  shielding,  abdomen  scan. 


Figure  J.8  Exposure  (R/s/parOcle)  at  each  source  point  for  -y  room  axis  detector 
for  CT  in  WPMC  with  1/8”  lead  shielding,  abdomen  scan. 
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Appendix  J:  WPMC  CT  Suite,  I/S'*  Lead  Shielding,  Abdomen  Scan 


Figure  J.9  Exposure  (R/s/parttde)  at  each  source  point  for  ceUng  detector 
for  CT  in  WPMC  with  1/8”  lead  shielding,  abdomen  scan. 


Figure  J.IO  Exposure  (R/s/pardcle)  at  each  source  point  for  floor  detector  for 
CT  in  WPMC  with  1/8”  lead  shielding,  abdomen  scan. 
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Appendix  J:  WPMC  CT  Suite,  1/8”  Lead  Shielding,  Abdomen  Scan 


Figure  J.ll  Normalized  spectra  of  weighted  average  at  the  positive  and  negative  x  source  axes  for 
the  CT  at  WPMC  with  1/8”  lead  shielding,  abdomen  scan. 


Figure  J.ll.  Normalized  spectra  of  wei|^ted  average  at  the  positive  and  negative  x  room  axes  for 
the  CT  at  WPMC  with  1/8”  lead  shielding,  abdomen  scan. 
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Appendix  J;  WPMC  CT  Suite,  1/8”  Lead  Shielding,  Abdomen  Scan 


Figure  J.13  Normalised  spectra  of  wdghted  avouge  atthe  positive  and  negatative  y  body  axes  for 
flieCT  at  WPMC  wifli  1/8”  toad  shielding,  abdomm  scan. 


Figure  J.14  Nmmalized  spectra  of  weighted  average  at  tte  positive  and  negatative  y  body  axes  for 
aie  CT  at  WPMC  with  1/8”  toad  shielding,  abdomen  scan. 
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Appendix  J:  WPMC  CT  Suite,  1/8”  Lead  Shielding,  Abdomen  Scan 


Figure  J.IS  Normalized  spectra  of  wei^ted  average  at  the  positive  and  negative  y  body  axes  for  tiie 
CT  at  WPMC  with  1/8"  lead  shielding,  abdomen  scan. 


Appendix  K:  GE  Minimiim  Recommended  Suite,  1/16”  Lead  Shielding,  Head  Scan 


Figure  K.2  Exposure  (R/s/particle)  at  each  source  point  for  negative  x  source  axis  detector  with  with 
CT  in  GE  minimum  recommended  room,  1/16**  lead,  head  scan. 
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Appendix  K:  GE  Minimum  Recommended  Suite,  1/16”  Lead  Shielding,  Head  Scan 


Figure  K.3  Exposure  (R/s/particle)  at  each  source  point  for  positive  x  room  axis  detector  with  with  CT 
in  GE  minimum  recommended  room,  \I16”  lead,  head  scan. 


in  GE  minimum  recommended  room,  1/16"  lead,  head  scan. 
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Appendix  K:  GE  Minimum  Recommended  Suite,  1/16”  Lead  Shielding,  Head  Scan 


Figure  K.5  Exposure  (R/s/purticle)  «t  each  source  point  for  positive  y  body  axis  detector  with  with  CT 
in  GE  minimum  recommended  room,  1/16"  lead,  head  scan. 
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Appendix  K:  GE  Minimum  Recommended  Suite,  1/16”  Lead  Shielding,  Head  Scan 


Figure  K.7  Exposure  (R/s/particle)  at  each  source  point  for  positive  y  room  axis  detector  with  CT  in  GE 
minimum  recommended  room  suite,  1/16"  lead,  head  scan. 
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Figure  K.8  Elxposure  (R/s/particle)  at  each  source  point  for  negative  y  room  axis  detector  with  with  CT 


in  GE  minimum  recommended  room,  1/16"  lead,  head  scan. 
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Appendix  K:  G£  Minimum  Recommended  Suite,  1/16”  Lead  Shielding,  Head  Scan 

6.0E-20 

5.0E-20 

4.0E-20 

3.0E-20 

2.0E-20 

l.OE-20 

O.OE+00 

0  45  90  135  180  225  270  315  360 

Source  Position  (deg) 

Figure  K.9  Exposure  (R/s/particle)  at  each  source  point  for  ceiling  detector  with  with  CT  in  GE 
minimuiii  recommended  room,  1/16*'  lead,  head  scan. 
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Figure  K.10  Exposure  (R/s/particle)  at  each  source  point  for  floor  detector  with  CT  in  GE  miniraum 
recommended  room,  1/16"  lead  shielding,  head  scan. 
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Figure  K.11  Normalized,  weighted-average  spectra  for  positive  and  negative  x  source  axes  detectors 
with  CT  in  GE  minimum  recommended  room,  1/16*'  lead,  head  scan. 


Figure  K.12  Normalized,  weighted-average  spectra  for  positive  and  negative  x  room  axes  detectors  with 
CT  in  G£  minimum  recommended  room,  1/16"  lead,  head  scan. 


Appendix  K:  GE  Minimum  Recommended  Suite,  1/16'*  Lead  Shielding,  Head  Scan 
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Figure  K.15  Normalized,  weighted-average  spectra  for  ceiling  and  floor  detectors  with  CT  in  GE 


minimum  recommended  room,  1/16"  lead,  head  scan. 
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